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ABSTRACT 


Studies  of  the  basic  properties  of  small  structures  that  could  be  of  major  importance 
in  the  fabrication  of  new  microdevices  have  been  carried  out.  Some  of  the  basic 
areas  of  major  contribution  include:  transport  over  and  through  barrier  structures, 
ohmic  contacts  and  fluctuations,  band  offsets,  and  the  search  for  lattice  matches  be¬ 
tween  silicides  and  silicon  for  application  in  microdevices..  Seven  graduate  students 
received  all  or  part  of  their  support  on  this  program.  Four  of  them  have  obtained 
Ph.  D.’s.  The  other  three  are  still  completing  their  programs. 


I.  Introduction 


In  recent  years,  it  has  become  quite  clear  that  the  current  trend  in  the  de¬ 
velopment  of  modern  electronics  cannot  continue.  The  strategy  of  continuing  to 
shrink  standard  device  structures  to  ever  smaller  dimensions  to  produce  circuits 
of  increasing  density  and  complexity  will  begin  to  fail  in  a  few  more  generations. 
Hence,  it  has  become  important  to  carry  out  major  research  programs  looking  for 
new  device  phenomena  and  exploring  the  physics  and  chemistry  of  the  ultimate 
limitations  in  modern  electronics.  The  research  program,  supported  under  this  con¬ 
tract,  is  one  such  program  trying  to  make  major  progress  in  understanding  some 
of  the  key  phenomena  that  could  be  important  both  in  new  device  structures  and 
also  as  an  ultimate  limiter  of  the  current  device  technology. 

Transport  is  a  major  issue  in  the  limitations  on  current  devices  and  possibili¬ 
ties  for  new  devices.  In  this  research  program,  we  have  explored  transport  through 
GaAs-GaAlAs  heterostructures.  Particular  emphasis  has  been  placed  on  tunnel 
barrier  structures  which  have  the  possibililty  for  being  important  in  new  device 
structures.  We  have  also  explored  the  role  of  the  discreteness  of  atoms  in  limiting 
the  ultimate  performance  of  such  technically  important  aspects  of  devices  as  ohmic 
contacts.  Transport  over  barriers  of  GaAlAs  have  been  explored  using  photore¬ 
sponse  techniques.  The  precise  nature  of  band  offsets  which  critically  determine 
the  behavior  of  heterojunctions  has  also  been  explored.  Finally,  we  have  investi¬ 
gated  the  possibililty  of  new  and  novel  lattice-matched  structures  involving  silicon 
and  silicides. 

Major  results  of  this  study  are  contained  in  the  next  section.  An  appendix 
contains  all  of  the  publications  that  were  supported  under  this  program.  One  of 
the  major  outputs  of  a  research  program  such  as  this  one  is  graduate  students.  A 
list  of  the  students  supported  and  their  current  positions  is  contained  in  Section  II. 
Overall  conclusions  are  contained  in  Section  IV. 
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II.  Major  Results 


A.  Tunneling 

The  primary  research  area  was  the  study  of  transport,  particularly  tunneling 
through  very  narrow  barrier  structures  with  characteristic  thicknesses  of  50-100A. 

One  of  the  primary  questions  is  the  role  of  inelastic  tunneling.  In  these  processes 
phonons  are  admitted  in  the  transport  across  the  barrier.  In  a  series  of  experiments 
carried  out  under  this  program,  we  made  some  of  the  first  investigations  of  inelastij 
tunneling  processes  in  these  simple  device  structures.  Major  results  of  this  study 
were  presented  in  papers  1  and  2  in  the  Appendix. 

We  have  had  major  success  in  studying  what  has  become  a  very  important 
device  structure,  the  so-called  double  barrier  resonant  structure.  We  were  the 
first  group  to  observe  resonant  tunneling  metal-organic-chemical-vapor  deposited 
GaAs-AlAs-GaAs-AlAs-GaAs  heterostructures.  We  were  the  third  group  to  ob¬ 
serve  this  phenomenon.  We  have  also  shown  that  the  uniformity  of  such  tunneling 
phenomenon  can  be  very  good,  hence  opening  up  the  possibilities  for  actually  fab¬ 
ricating  meaningful  devices  with  this  technique.  These  very  important  results  are 
contained  in  papers  3  and  4  in  the  Appendix. 

We  have  also  investigated  the  fundamental  transport  mechanisms  in  these  dou¬ 
ble  barrier  structures.  In  particular,  we  have  investigated  the  role  of  transport  in¬ 
volving  other  parts  of  the  band  structure  in  addition  to  the  T-point  as  well  as  the  j 

role  of  phonons  in  the  transport.  These  results  are  contained  in  paper  5  in  the  ! 

Appendix.  ! 

I 

I 

I 

3 

i 

i 


3 


While  very  exciting  two-terminal  devices  have  been  fabricated  based  on  the  res¬ 
onant  tunneling  through  a  double  barrier  structure,  for  these  new  device  structures 
to  be  very  useful  a  third-terminal  is  essential.  Over  the  years,  tunneling  based  tran¬ 
sistors  have  been  proposed  but  have  had  great  difficulity  in  attaining  performance 
levels  which  would  make  them  interesting.  Under  this  research  project,  we  have 
proposed  a  novel  three-terminal  device  based  on  tunneling.  This  device  structure 
called  the  Stark  transistor  has  the  possibility  of  overcoming  some  of  the  historical 
problems  of  other  tunneling  transistor  structures.  While  we  have  not  had  success 
in  fabricating  this  device  under  this  program,  we  believe  that  this  proposed  device 
structure  opens  up  the  possibility  of  making  high-density,  high-speed  electronics 
based  on  tunnel  structures.  We  report  on  these  new  device  structures  in  paper  6  in 
the  Appendix. 

B.  Transport  Over  Barriers 

In  addition  to  the  transport  through  the  barrier,  reported  in  the  above  studies, 
we  have  also  investigated  the  transport  over  barriers  in  GaAs-GaAlAs-GaAs  het¬ 
erostructure  devices.  In  order  to  carry  out  these  studies,  we  have  used  photons  to 
excite  the  carriers  over  the  barrier.  While  the  initial  plans  for  these  experiments 
were  to  explore  the  details  of  resonances  in  the  transfer  of  electrons  across  barri¬ 
ers,  theoretical  studies  as  well  as  the  experiments  indicated  that  these  could  not 
be  observed  due  to  the  distribution  in  k  of  the  electrons  excited  by  the  photons. 
However,  a  number  of  novel  effects  were  observed  and  explained  theoretically  during 
this  study.  The  results  are  contained  in  papers  7,  8  and  9  in  the  Appendix. 


C.  Ohmic  Contacts  and  Fluctuation  Phenomena 


Ohmic  contacts  are  one  of  the  most  important  components  of  modern  devices. 
As  device  dimensions  that  are  shrunk  to  increase  density  and  performance,  ohmic 
contacts  are  likely  to  be  a  major  limiter  on  device  performance.  Ohmic  contacts  are 
made  using  black  art  recipes  which  give  little  insight  into  exactly  how  they  work. 
However,  it  is  widely  believed  that  the  ohmic  contact  actually  consists  of  a  heavily 
doped  region  in  the  semiconductor  in  which  tunneling  dominates  transport  through 
a  normally  present  Schottky  barrier.  There  have  been  very  few  theoretical  studies  of 
ohmic  contacts.  The  most  popular  theory  is  one  that  makes  a  number  of  simplifying 
assumptions  which  made  comparision  between  specific  contact  resistance  and  theory 
impossible.  Under  this  program,  we  have  developed  the  first  serious  theoretical 
study  of  ohmic  contact  behavior.  This  study  was  specifically  aimed  at  GaAs  but 
is  readily  extendable  to  other  systems.  The  results  of  this  study  are  presented  in 
paper  10  in  the  appendix. 

One  of  the  major  issues  in  examining  the  properties  of  ohmic  contacts  is  the 
discreteness  of  the  dopant  atoms.  Most  of  the  theories,  including  the  one  described 
in  the  preceeding  paragraph,  assume  that  the  dopant  atoms  are  uniformily  dis¬ 
tributed.  However,  in  actual  fact  there  are  large  fluctuations  in  potential  due  to  the 
discreteness  or  the  ionized  dopants.  During  this  program,  we  have  initiated  a  study 
of  this  fundamental  phenomenon  to  study  the  fluctuations  that  could  occur  in  the 
ohmic  contact  behavior.  This  pioneering  work  is  beginning  to  illustrate  the  ultimate 
limitations  on  shrinking  devices  due  to  the  discreteness  of  the  actual  underlaying 
atomic  structure.  The  preliminary  results  of  this  study  are  contained  in  paper  11 
in  the  Appendix. 
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D.  Band  Offsets 


l 


One  of  the  major  unsolved  problems  in  the  study  of  heterojunctions  is  the  de¬ 
termination  of  the  band  offset.  The  band  offset  measures  the  discontinuity  between 
the  valence  bands  or  conduction  bands  of  the  two  semiconductors  on  either  side 
of  the  heterojunction.  Within  the  last  few  years,  it  has  been  one  of  the  primary 
sports  of  the  theorists  in  this  field  to  attempt  to  explain  the  values  of  band  offsets. 
Under  this  program,  we  made  an  important  contribution  to  this  field  indicating 
the  kinds  of  effects  that  could  result  because  of  defects  and  dipole  layers  at  these 
heterojunctions.  Results  of  our  study  are  included  in  the  Appendix  in  paper  12. 


E.  Lattice  Match  Structure  of  Silicides  and  Silicon 
Lattice  match  is  thought  to  play  a  major  role  in  growing  device  quality  in¬ 
terfaces.  One  of  the  major  areas  of  interest  is  the  fabrication  of  device  structures 
involving  silicides  and  silicon.  During  this  contract,  we  made  a  major  search  for 
lattice  matches  between  silicide  and  silicon.  Results  of  this  study  are  contained  in 
paper  13  in  the  Appendix. 


III.  Students  Supported  Under  this  Contract 

One  of  the  major  outputs  of  any  research  program  in  a  university  environment  is 
the  graduate  students.  During  this  program  the  research  of  seven  graduate  students 
has  been  supported.  A  table  of  these  graduate  students  and  their  current  positions 
is  given  below. 
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Name 

Degree 

Year 

Current  Position 

R.  T.  Collins 

Ph.D 

1984 

Member  Technical  Staff,  IBM  T.  J.  Watson  Research 
Center 

A.  Zur 

Ph.D 

1983 

Member  Technical  Staff,  Rafphael,  Israel 

A.  R.  Bonnefoi 

Ph.D 

1986 

Prof,  of  Physics,  Uruguay 

T.  E.  Schlesinger  Ph.D 

1985 

Assistant  Professor  of  Electrical  Engineering,  Carnegie 
Mellon  Univ. 

W.  Boudville 

Graduate  Student  at  Caltech 

T.  K.  Woodward 

Graduate  Student  at  Caltech 

D.  H.  Chow 

Graduate  Student  at  Caltech 

IV.  Conclusions 

During  this  research  program,  major  results  were  obtained  in  the  study  of 
tunneling  structures  and  the  fabrication  of  two-terminal  resonant  tunneling  devices. 
We  also  proposed  a  novel  three-terminal  tunneling  device  which  has  the  possibility 
of  giving  both  high-speed  and  high-density  electronics.  Additionally,  we  carried 
out  extensive  studies  of  the  photovoltaic  properties  of  GaAs-GaAlAs-GaAs  barrier 
structures.  These  represented  some  of  the  first  studies  of  the  transport  across  these 
barriers.  Ohmic  contacts,  which  are  likely  to  be  a  major  limiter  in  1 1 1  -  V  device 
technology,  have  been  explored  both  in  a  continuum  theory  and  one  including  the 
fluctuations  due  to  dopants.  The  role  of  defects  and  dipole  layers  in  producing  band 
offsets  has  been  studied  in  some  very  simple  theoretical  investigations.  Lattice 
match  between  silicon  and  silicide  structures  has  been  studied  in  an  attempt  to 
identify  new  and  potentially  very  interesting  device  structures.  Under  the  program, 
seven  graduates  students  have  received  support.  Four  of  these  graduate  students 
have  obtained  degrees  and  are  holding  down  positions  in  academic  institutions  or 
in  industrial  research  laboratories.  Three  of  the  students  are  continuing  to  work 
towards  degrees  at  Caltech. 
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Inelastic  tunneling  characteristics  of  AlAs/GaAs  heterojunctions 

R.  T  Collins,  J.  Lambe,  and  T.  C.  McGill 

California  Institute  of  Technology.  Pasadena.  California  9/125 

R.  D.  Burnham 

Xerox  Corporation.  Palto  Alto.  California  94304 

(Received  7  November  1983;  accepted  for  publication  12  December  1983) 

We  report  the  first  observation  of  inelastic  tunneling  in  electronic  transport  perpendicular  to  a 
thin  AlAs  layer  sandwiched  between  two  GaAs  layers.  Temperature  dependent  I-V,  first 
derivative  [dl  /dV),  and  second  derivative  \d2I  /dV2)  measurements  were  made  on  AlAs/GaAs 
double  heterojunctions  for  a  range  of  AlAs  layer  thicknesses  and  dopings.  For p-type  AlAs 
harriers  current  transport  at  4.2  K  was  due  to  tunneling,  and  reproducible  structure  was  seen  in 
the  second  derivative  spectrum.  This  structure  was  associated  with  the  inelastic  excitation  of 
AlAs  optical  phonons  and  with  a  density-of-states  effect  caused  by  optical  phonon-electron 
coupling  in  the  GaAs.  A  different  second  derivative  spectrum  which  also  exhibited  reproducible 
structure  was  obtained  for  rt-type  AlAs  layers.  Several  possible  explanations  for  these  differences 
are  proposed. 

PACS  numbers:  73  40. Gk,  73.40. Lq,  73.60. Fw 

Several  studies  of  the  electrical  properties  of  the  ducibie  structure  was  seen  in  the  second  derivative  spectrum 

Al,Ga,  _ ,  As/GaAs  double  heterojunctions  have  recently  at  the  AlAs  and  GaAs  optical  phonon  frequencies.  These 

been  reported.  In  these  papers  the  capacitance-voltage  (C-  structures  were  attributed  to  the  inelastic  excitation  of  AlAs 

Fjandcurrent-voltage  (I-V)  properties  of  thin  Al,  Ga,  ,  As  longitudinal  optical  phonons  and  to  a  density-of-states  effect 

barriers  sandwiched  between  two  GaAs  layers  were  in  vest  i-  caused  by  coupling  between  electrons  and  optical  phonons  in 

gated.  The  currents  obtained  in  the  I-V  measurements  were  the  GaAs.  In  contrast,  I-V  curves  and  second  derivative 

attributed  to  thermionic  emission  over  the  offset  barrier  and  spectra  for  samples  with  n-type,  Se-doped  AlAs  layers  were 
tunneling  through  the  barrier.  Chang  et  al 4  and  Solner  eta/.'  considerably  different.  Possible  explanations  for  these  differ- 
looked  at  the  tunneling  behavior  of  two  Al,Ga,  ,  As  bar-  ences  will  be  given. 

ners  separated  by  a  GaAs  quantum  well  and  Vojak  et  al.''  The  samples  used  ir.  this  study  were  grown  by  melalor- 

observed  tunneling  behavior  in  multiple  quantum  well  struc-  ganic  chemical  vapor  deposition  (MOCVD)  using  a  tech- 

tures,  but  to  our  knowledge,  no  one  has  looked  for  structure  mque  described  elsewhere.7  The  substrate  used  in  making 

in  the  tunneling  /-Fcurve  of  a  single  AI,Ga,  . ,  As  barrier.  the  structures  was /r-type  GaAs  doped  with  2-3  X  1018  era'1 

[n  this  letter  we  report  the  first  observation  of  inelastic  Si  An  approximately  3-pm-thick  epitaxial  layer  of  GaAs 

tunneling  in  electronic  transport  perpendicular  to  a  single  doped  with  Se  was  grown  on  the  substrate  followed  by  a  thin 

AlAs  barrier  Temperature-dependent  I-V.  first  derivative  AlAs  layer.  The  thicknesses  of  the  AlAs  layers  varied  from 

I  dl  /dV ),  and  second  derivative  [d2I  /dV2)  measurements  60  to  250  A  from  wafer  to  wafer.  A  second  Se-doped  GiAs 

were  made  on  AlAs/GaAs  double  heterojunctions  Both  the  layer  was  grown  on  top  of  the  AlAs  barrier.  The  epitaxial 

thickness  and  the  doping  of  the  AlAs  layers  were  varied.  GaAs  dopings  varied  but  were  typically  1-5  X  10IS  cm  ‘  \ 

When  the  AlAs  layer  was  doped p  type  with  Mg,  transport  at  The  AlAs  was  either  doped  n  type  at  approximately  1  X  I018 

4  2  K  occurred  mainly  through  electron  tunneling.  Repro-  cm  '  with  Se  or p  type  at  about  1  x  10IS  cm  "  '  with  Mg  The 
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widths  and  dopings  of  the  AlAs  given  above  are  extrapola¬ 
tions  from  measurements  made  on  low  aluminum  content 
Al,  Ga,  _  ,  As  and  could  be  in  error  since  they  have  not  actu¬ 
ally  been  calibrated  for  AlAs. 

Circular  mesas  80-360 /mi  in  diameter  and  more  than  4 
/j.m  in  height  were  defined  on  the  epitaxial  GaAs  surface 
using  photolithography  and  a  4: 1 : 1  H,S04:H:0::H;0  GaAs 
etch.  A  (Au.Ge)-Ni  ohmic  contact  was  fabricated  on  the 
mesa  surface  and  on  the  substrate.  The  samples  were  then 
mounted  on  standard  transistor  headers  and  wire  bonded  to 
provide  contacts. 

Temperature-dependent  I-V first  derivative,  and  sec¬ 
ond  derivative  spectra  were  taken  on  the  samples  from  room 
temperature  to  4.2  K.  The  I-V  measurements  were  made  by 
applying  a  dc  current  source  to  the  samples  and  measuring 
the  resulting  dc  voltage.  First  and  second  derivatives  were 
obtained  using  standard  modulation  techniques.  First  deri¬ 
vative  spectra  were  taken  by  applying  a  5-kHz  current 
source  to  the  samples  and  measuring  the  voltage  across  the 
sample  at  5  kHz  to  obtain  dV /dl.  Second  derivative  spectra 
were  taken  using  a  second  harmonic  detection  system  as  de¬ 
scribed  by  Lambe  and  Jaklevic.8  It  was  operated  at  a  funda¬ 
mental  frequency  of  50  kHz  with  typical  modulation  vol¬ 
tages  of  1-5  mV. 

First,  the  results  of  measurements  on  samples  with  Mg- 
doped,  //-type  AlAs  barriers  will  be  presented.  The  I-V 
curves  for  p-type  barriers  exhibited  increasing  zero-bias  re¬ 
sistance  as  temperature  was  decreased  from  room  tempera¬ 
ture  to  77  K.  From  77  to  4.2  K.  the  curves  maintained  the 
same  shape  and  resistance.  As  the  thickness  of  the  AlAs 
layer  was  increased  from  60  to  120  A,  the  zero-bias  resistiv¬ 
ity  increased  by  about  2.5  orders  of  magnitude.  The  depen¬ 
dence  of  sample  resistance  on  temperature  and  AlAs  layer 
thickness  strongly  suggests  that  current  flow  at  low  tempera¬ 
tures  was  mainly  due  to  electronic  tunneling  through  the 
AlAs  barrier.  If  this  were  true,  the  zero-bias  resistivity 
should  have  an  approximately  exponential  dependence  on 


-200  -100  0  100  200 


V  (mV) 

FIG  I  Second  derivative  spectrum  at  4  2  K  for  a  p-\ ype  AlAs  layer  sand¬ 
wiched  between  two  w-type  GaAs  layers  The  layer  thickness  was  approxi¬ 
mately  60  A  The  structure  labeled  (al  is  attributed  to  an  electronic  density- 
of-states  effect  caused  by  electron -optical  phonon  coupling  in  the  GaAs 
The  peak  labeled  |b|  is  caused  by  the  inelastic  excitation  nf  longitudinal  opti¬ 
cal  phonons  in  the  AlAs  barrier 
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barrier  width.  Additional  60  A  increases  in  barrier  width  did 
result  in  greater  zero-bias  resistivity  but  by  less  than  2.5  or¬ 
ders  of  magnitude.  It  is  likely  that  impurity  states,  edge  ef¬ 
fects,  and  other  forms  of  leakage  in  the  AlAs  contribute  to 
the  low-temperature  current.  For  small  thicknesses  the  cur¬ 
rent  would  be  dominated  by  tunneling  through  the  barrier, 
while  for  larger  thicknesses  impurity  conduction  and  other 
forms  of  leakage  would  be  more  important. 

A  second  derivative  spectrum  taken  at  4.2  K  for  a  60- A- 
thick  p-type  AlAs  barrier  is  given  in  Fig.  1.  Two  structures 
are  present  in  both  forward  and  reverse  bias.  The  first,  la¬ 
beled  (a),  is  at  approximately  the  optical  phonon  frequencies 
of  GaAs.  It  resembles  similar  structures  previously  observed 
in  tunneling  through  GaAs  Schottky  barriers  which  were 
attributed  to  a  density-of-states  effect  in  the  degenerate 
GaAs  layers.0  This  self-energy  effect  is  caused  by  coupling 
between  electrons  and  longitudinal  optical  phonons  in  the 
GaAs.  The  structure  labeled  (a)  is  believed  to  arise  from  the 
same  effect.  The  peak  labeled  |b)  occurs  at  50  meV  which  is 
the  longitudinal  optical  phonon  frequency  of  AlAs.10  This 
peak  is  attributed  to  the  inelastic  excitation  of  longitudinal 
optical  phonons  in  the  AlAs  by  the  electrons  as  they  tunnel 
through  the  barrier.  This  is  the  first  observation  of  these 
effects  in  the  AlAs/GaAs  system.  The  same  results  were 
obtained  in  two  different  wafers  with  approximately  60-A- 
thick  AlAs  barriers.  For  barriers  1 20  A  thick  or  greater,  this 
structure  was  not  observed  in  the  second  derivative  spec¬ 
trum.  This  is  probably  due  to  the  large  decrease  in  tunneling 
current  at  these  thicknesses  which  reduced  the  second  deri¬ 
vative  signal  below  the  detectable  limit  of  the  measurement 
system. 

When  the  measurements  described  above  were  made  on 
samples  with  Se-doped,  n- type  AlAs  layers,  the  results  were 
not  the  same.  I-V  curves  for  the  n-type  AlAs  layers  also 
displayed  increasing  resistance  with  decreasing  tempera¬ 
ture.  In  contrast  to  the  p-type  layers,  the  resistance  contin¬ 
ued  to  increase  until  approximately  30  K.  Furthermore.the 
zero-bias  resistivity  increased  by  less  than  an  order  of  magni¬ 
tude  as  the  AlAs  layer  thicknesses  increased  from  60  to  250 
A,  and  the  current  flow  was  nearly  an  order  of  magnitude 
larger  than  in  the  60-A  p-type  barrier  (Fig.  2).  Changes  in 
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FIG  2  /-land  dl  /«/r  curves  at  4  2  K  for  an  n-type  AlAs  barrier  sand¬ 
wiched  between  two  n-type  GaAs  layers  The  AlAs  layer  thickness  was  ap¬ 
proximately  80  A 
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FIG  3.  Second  derivative  spectrum  at  4  2  K  for  an  approximately  8U-A- 
thick,  n-type  AlAs  layer  sandwiched  between  two  n-type  GaAs  layers 

slope  were  visible  in  the  I-V  curves  taken  at  4.2  K.  These 
changes  in  slope  were  as  large  as  60%  as  can  be  seen  from  the 
dl /dV curve  given  in  Fig.  2.  To  study  these,  second  deriva¬ 
tive  spectra  were  taken. 

A  second  derivative  spectrum  at  4.2  K  for  an  n-type 
barrier  device  is  shown  in  Fig.  3.  Five  peaks  representing 
increases  in  conductivity  are  evident  in  the  positive  bias  di¬ 
rection.  The  two  most  prominent  peaks  occur  at  approxi¬ 
mately  50  and  70  meV.  Two  peaks  are  also  seervat  about  1 25 
and  145  meV.  In  addition,  there  is  a  23-meV  shoulder  pres¬ 
ent  on  a  strong  zero-bias  anomaly.  The  spectrum  is  fairly 
symmetric  with  respect  to  bias  direction.  Nearly  identical 
second  derivative  spectra  were  observed  for  all  AlAs  thick¬ 
nesses.  The  same  peaks  were  seen,  but  there  were  differences 
in  the  measured  voltages  for  the  peaks  and  asymmetries 
between  the  voltages  measured  in  both  bias  directions.  These 
asymmetries  and  voltage  differences  were  usually  around  5- 
10  meV. 

A  definitive  explanation  for  the  differences  between 
samples  with  n-type  and  p-t ype  barriers  has  not  been  found, 
but  some  possible  explanations  can  be  proposed.  If  transport 
in  the  n-type  AlAs  samples  is  the  result  of  tunneling,  a  model 
for  these  differences  must  describe  why  the  actual  tunneling 
barrier  seen  by  the  electrons  is  not  increasing  with  the  phys¬ 
ical  width  of  the  AlAs  layer  in  the  n-type  samples,  as  it  is  in 
the/>-type  samples.  One  possible  explanation  is  that  the  dop¬ 
ing  in  the  n-type  AlAs  layer  is  actually  larger  than  anticipat¬ 
ed,  and  band  bending  in  the  AlAs  results  in  the  formation  of 
two  back-to-back  AlAs/GaAs  tunneling  barriers  with  a  con¬ 
necting  layer  of  AlAs.  Increasing  the  AlAs  layer  thickness 
would  only  increase  the  width  of  the  connecting  AlAs  layer. 
Then,  the  peaks  in  the  second  derivative  spectrum  may  arise 
from  the  excitation  of  characteristic  phonons  in  the  AlAs/ 
GaAs  system  by  tunneling  electrons  (possibly  wave  vector- 
conserving  phonons  resulting  from  electrons  moving  from 


direct  to  indirect  conduction-band  minima).  A  voltage  drop 
would  occur  across  each  barrter,  causing  the  voltages  at 
which  the  phonon  associated  peaks  occur  to  differ  from  the 
phonon  energies. 

A  second  explanation  includes  the  possibility  that  more 
impurities  or  defects  are  incorporated  into  the  n-type  layers 
than  into  the  p-type  layers,  giving  rise  to  impurity  conduc¬ 
tion  in  the  barrier.  This  would  not  have  an  exponential  de¬ 
pendence  on  barrier  width.  The  peaks  in  the  second  deriva¬ 
tive  spectrum  could  indicate  the  onset  of  impurity-assisted 
tunneling  or  the  excitation  of  vibrational  modes  of  the  im¬ 
purities.  Work  is  currently  under  way  to  test  these  models  of 
current  transport  in  the  n-type  barriers  samples. 

To  summarize,  we  have  observed  current  transport  per¬ 
pendicular  to  a  thin  AlAs  layer  sandwiched  between  two 
degenerately  doped,  n-type  GaAs  layers.  Dependences  of 
the  associated  I-V  curves  on  temperature  and  thickness  of 
the  AlAs  layer  seem  to  indicate  that  for p-type  AlAs  barriers 
current  flow  at  4.2  K  is  largely  due  to  tunneling  through  the 
AlAs  barrier.  Reproducible  structure  was  observed  in  the 
second  derivative  spectrum  for p-type  AlAs  barriers  and  as¬ 
sociated  with  the  inelastic  excitation  of  longitudinal  optical 
phonons  in  the  AlAs  and  a  density-of-states  effect  in  GaAs 
caused  by  coupling  between  the  electrons  and  optical  phon¬ 
ons.  Samples  with  n- type  AlAs  barriers  did  not  have  the 
same  characteristics.  The  resistivity  of  the  samples  was  near¬ 
ly  independent  of  AlAs  layer  thickness,  and  different  peaks 
were  observed  in  the  second  derivative  spectrum.  These 
peaks  resulted  from  very  large  changes  in  the  slope  of  the  I-  V 
curves  and  were  very  reproducible.  A  definitive  explanation 
for  the  diffferences  between  n-type  and  p- type  barriers  has 
not  been  found,  but  some  possibilites  have  been  proposed. 
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We  report  the  observation,  at  low  temperatures  and  bias  vol¬ 
tages,  of  elastic  and  inelastic  tunneling  through  AlAs  bar¬ 
riers  sandwiched  between  two  GaAs  layers.  Previous  studies 
of  electronic  transport  perpendicular  to  Al„  Ga,  _  x  As  bar¬ 
riers  sandwiched  between  GaAs  layers  ,~4  have  attributed 
currents  passing  through  the  barriers  to  thermionic  emission 
at  temperatures  above  150  K  or  to  tunneling  through  the 
Al,  Gai  _  ,  As  barrier  at  low  temperatures,  but  for  bias  vol¬ 
tages  which  were  comparable  to  the  barrier  height.  Two  of 
these  studies  noted  that  at  low  temperatures  and  low  biases 
more  current  passed  through  the  barriers  than  could  be  ac¬ 
counted  for  by  either  of  these  processes,  indicating  that  some 
form  of  leakage  may  be  dominating  current  transport  at  low 
temperatures.1-2  In  all  of  these  studies  the  barriers  were  300 
A  thick  or  thicker.  The  measurements  reported  here  were 
made  on  samples  with  barriers  thinner  than  300  A  in  an 
effort  to  enhance  the  tunneling  current  at  low  bias  voltages 
relative  to  leakage.  A  more  detailed  account  of  most  of  these 
measurements  is  given  in  Ref.  5. 

The  samples  used  in  this  study  were  grown  by  metalor- 
ganic  chemical  vapor  deposition  (MOCVD)  and  by  molecu¬ 
lar  beam  epitaxy  (MBE).  In  both  cases  an  «-type  epitaxial 
layer  of  GaAs  was  grown  on  a  degenerately  doped  n-type 
GaAs  substrate.  A  thin  layer  of  AlAs  was  grown  on  top  of 
this  GaAs  layer,  followed  by  a  final  n-type  layer  of  GaAs.  In 
the  MOCVD  samples  the  epitaxial  GaAs  was  doped  with  Se 
at  1-5 X  10'®  cm-3,  and  the  AlAs  layer  was  either  doped  n 
type  with  Se  or  p  type  with  Mg  at  approximately  1  X  10'® 
cm'3.  The  AlAs  layers  were  from  50-300  A  thick.  The  MBE 
samples  came  from  two  sources.  In  both  cases  the  epitaxial 
GaAs  layers  were  doped  with  Si  at  1 X  10'*  cm-3,  and  the 
AlAs  was  unintentionally  doped.  The  estimated  AlAs  layer 
thicknesses  were  from  50-150  A.  Devices  were  made  on  the 
wafers  by  defining  mesas  and  fabricating  (Au,Ge)-Ni  con¬ 
tacts  on  the  mesas.  Temperature  dependent  current-voltage 
(/-K),  first  derivative  (dl /dV),  and  second  derivative  [d 2 / / 
dV2)  spectra  were  taken  on  the  devices  from  room  tempera¬ 
ture  to  4.2  K.  Zero-bias  resistances  were  obtained  from  the 
slope  of  the  7-  V  curves  at  ±  3  mV. 

Results  of  the  measurements  made  on  MOCVD  grown 
samples  will  be  discussed  first.  There  was  a  large  difference 
in  the  characteristics  of  the  samples,  depending  upon 
whether  the  AlAs  barrier  was  doped  p  type  with  Mg  or  n 
type  with  Se.  The  p- type  barriers  exhibited  increasing  resis¬ 
tance  as  temperature  was  decreased  to  77  K,  at  which  point 
the/-Fcurve  maintained  the  same  shape  down  to  4.2  K.  The 


zero-bias  resistances  of p-type  barrier  samples  increased  ra¬ 
pidly  with  barrier  width  as  is  seen  in  Fig.  1.  Figure  1  also 
includes  the  results  of  a  WKB  calculation  of  the  zero-bias 
resistances  expected  for  various  barrier  widths.  This  calcula¬ 
tion  was  based  on  the  technique  of  Kurtin  et  al.6  Only  contri¬ 
butions  to  the  tunneling  probability  from  the  AlAs  r  points 
were  included.  The  band  offsets  and  effective  masses  were 
taken  from  Casey  and  Panish.7  Although  this  model  is  only 
presented  to  give  an  idea  of  the  orders  of  magnitude  expected 
for  tunneling  resistances,  the  agreement  is  reasonable  for 
small  barrier  widths.  As  the  thickness  of  the  barrier  is  in¬ 
creased,  the  measured  zero-bias  resistance  does  not  increase 
at  the  rate  predicted  for  tunneling,  indicating  that  for  thick 
AlAs  layers,  most  of  the  low  temperature  current  is  due  to 
some  form  of  leakage.  At  smaller  thicknesses  the  tunneling 


Fig.  I .  Zero-bias  resistivity  at  4.2  K  for  MCXVD  samples  with/Mype  bar¬ 
riers  and  for  MBE  samples  which  had  uniform  resistivities  between  neigh¬ 
boring  devices.  The  solid  line  is  a  theoretical  calculation  of  zero-bias  resis¬ 
tivity,  using  a  WKB  approximation  to  the  tunneling  probability. 
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Fig.  2.  Second  derivative  spectrum  at  4.2  K  for  a  p- type  AlAs  barrier  sam¬ 
ple  grown  by  MOCVD.  The  layer  thickness  was  approximately  50  A.  The 
structure  labeled  (a|  is  attributed  to  an  electronic  density-of-states  effect 
caused  by  electron-optical  phonon  coupling  in  the  GaAs.  The  peak  labeled 
(b)  is  caused  by  the  inelastic  excitation  of  longitudinal  optical  phonons  in  the 
AlAs  barrier. 


current,  which  decreases  exponentially  with  width,  begins  to 
dominate  leakage. 

When  first  and  second  derivatives  were  taken  of  the  I~V 
curves  of  samples  with  50  A  thick,  p-type,  AlAs  layers,  re¬ 
producible  structure  was  evident  (Fig.  2).  Structure  [labeled 
(a)]  is  present  at  the  GaAs  longitudinal  optical  (LO)  phonon 
energy  (36  meV)  in  both  forward  and  reverse  bias  in  Fig.  2 
and  is  attributed  to  a  density-of-states  effect  caused  by  elec¬ 
tron-phonon  coupling  in  the  degenerate  GaAs  layers.  Simi¬ 
lar  observations  have  been  made  in  Schottky  barrier  tunnel¬ 
ing  in  GaAs.8  A  peak,  labeled  (bj,  is  also  visible  at 
approximately  the  AlAs  LO  phonon  energy  (50  meV)  in  the 
second  derivative  spectrum.9  This  peak  is  associated  with 
the  inelastic  excitation  of  longitudinal  optical  phonons  in  the 
AlAs  by  the  tunneling  electrons.  This  is  the  first  observation 
of  these  effects  in  the  AlAs/GaAs  system.  These  results  were 
reproducible  across  a  given  wafer  and  were  seen  in  two  wa¬ 
fers  with  barriers  of  about  the  same  thickness.  When  the 
AlAs  barrier  width  increased  to  100  A  or  greater,  no  repro¬ 
ducible  structure  was  observed. 

When  the  MOCVD  grown  AlAs  barriers  were  doped  n 
type,  the  results  were  different.  The  zero-bias  resistance 
changed  by  less  than  an  order  of  magnitude  as  the  AlAs  layer 
width  was  increased  from  60  to  250  A.  Five  peaks  were  pres¬ 
ent  in  the  second  derivative  spectra  in  both  forward  and  re¬ 
verse  bias  for  all  of  the  n-type  barrier  samples  regardless  of 
width.  These  peaks  did  not  match  those  seen  in  the  p- type 
barrier  MOCVD  samples.  A  definitive  explanation  for  the 


results  of  measurements  on  n-type  barriers  has  not  been 
found.  A  more  complete  discussion  of  the  differences 
between  n- and p-type  barrier  MOCVD  samples  and  the  pos¬ 
sible  explanations  for  these  differences  is  given  in  Ref.  5. 

The  zero-bias  resistances  for  devices  made  from  the  two 
MBE  grown  samples  with  the  thinnest  AlAs  layers  (one 
from  each  source)  are  given  in  Fig.  1.  These  resistances  are 
also  consistent  with  the  tunneling  calculations,  although 
from  Fig.  1,  it  appears  that  either  the  MBE  or  MOCVD 
barrier  widths  may  be  in  error  by  around  10  A.  The  I-V 
curves  for  devices  made  from  these  wafers  showed  approxi¬ 
mately  the  same  behavior  with  temperature  as  the  p-type 
MOCVD  barrier  samples.  Devices  made  on  the  two  wafers 
for  which  resistances  are  given  in  Fig.  1  showed  some  fluctu¬ 
ation  from  device  to  device  in  the  zero-bias  resistance.  The 
remaining  MBE  samples  studied  had  barrier  widths  of  100  A 
or  greater.  The  zero-bias  resistances  of  neighboring  devices 
made  on  these  wafers  fluctuated  by  many  orders  of  magni¬ 
tude.  The  largest  resistances  measured  for  the  MBE  samples 
were  two  or  more  orders  of  magnitude  larger  than  the  resis¬ 
tances  of  p-type  MOCVD  samples  with  the  same  barrier 
width.  Much  smaller  fluctuations  occurred  in  the  resistivi¬ 
ties  of  devices  on  the  MOCVD  grown  wafers.  The  first  and 
second  derivative  spectra  for  the  MBE  wafers  did  not  show 
any  reproducible  structure,  possibly  because  they  had  such 
large  resistances  that  tunneling  currents  were  beyond  the 
detectable  limits  of  our  system.  A  reason  for  the  large  flunc- 
tuations  in  device  resistance  for  the  thicker  barrier  MBE 
samples  has  not  been  found. 
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We  report  the  first  observations  of  resonant  tunneling  in  electronic  transport  perpendicular  to  two 
AlAs  layers  separated  by  a  GaAs  quantum  well  in  GaAs/AlAs  heterostructures  grown  by 
metalorganic  chemical  vapor  deposition.  Resonant  tunneling  can  be  observed  as  inflections  in  the 
I-  V curves  at  room  temperature.  These  inflections  become  more  pronounced  as  the  temperature  is 
reduced,  until  negative  differential  resistance  regions  become  visible  for  temperatures  below  260 
K.  At  low  temperatures,  the  I-  V curves  not  only  reveal  two  large  negative  resistance  regions 
corresponding  to  the  first  energy  level  in  the  GaAs  quantum  well  but  also  a  structure  which  shows 
evidence  of  resonant  tunneling  through  the  second  and  possibly  the  third  energy  states  in  the  well. 

Second  derivative  \d2I  /dV2)  measurements  confirm  the  existence  of  the  resonances  seen  in  the  /-  V 
curves. 

Semiconductor  heterostructures  with  small  character-  ber  of  studies  of  electronic  transport  perpendicular  to  the 

istic  dimensions  are  now  frequently  fabricated  from  GaAs  layers  have  been  reported. 1-6  Because  an  important  current 

and  Al,  Ga,  ,  As  using  techniques  such  as  molecular  beam  transport  mechanism  is  tunneling,  the  associated  tunneling 

epitaxy  (MBE)  and  metalorganic  chemical  vapor  deposition  effects  yield  novel  electronic  properties  with  very  important 

(MOCVDI.  The  conduction-band  offset  between  the  GaAs  theoretical  and  maybe  practical  applications.  An  interesting 

and  the  Al,  Ga,  _  ,  As  causes  Al,  Ga,  ,  As  layers  to  act  as  configuration  is  one  in  which  two  barriers  are  separated  by  a 

potential  energy  barriers  to  electrons  in  the  GaAs  layers.  quantum  well.  There  have  been  observations  of  resonant 

This  makes  it  possible  to  grow  quantum  well  and  quantum  tunneling  in  MBE  grown  samples  with  a  GaAs  well  adjacent 

barrier  structures,  which  are  the  object  of  considerable  theo-  to  two  Al,  Ga,_,As  barriers.7,8  Following  some  prelimi- 

retical  and  experimental  work.  For  these  structures,  a  num-  nary  work,‘>,'°  we  report  here  observations  of  resonant  tun- 
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neling  in  MOCVD  grown  structures  with  two  AlAs  harriers 
separated  by  a  single  GaAs  quantum  well.  For  the  first  time, 
resonant  tunneling  can  be  observed,  not  only  through  the 
lowest  energy  level  in  the  GaAs  well  at  room  temperature, 
but  also  through  higher  energy  states  at  lower  temperatures, 
once  the  thermal  current  has  been  reduced. 

The  structure  used  in  this  study  was  grown  by  an 
MOCVD  technique"  '  on  a  GaAs  substrate  doped  n  type 
with  Si  at  2  —  3  X  10'"  cm  The  first  layer  grown  on  the 
substrate  was  one  of  GaAs,  doped  n  type  with  Se  at  about 
2x  10IK  cm  1  and  with  a  thickness  of  2-3  fin\.  Then,  two 
thin  layers  of  AlAs,  doped  p  type  with  Mg"  14  at  approxi¬ 
mately  1  X  101"  cm  "  \  and  separated  by  a  thin  layer  of  nomi¬ 
nally  undoped  GaAs,  were  grown.  Finally,  a  1 ,8-/jm  top  lay¬ 
er  of  GaAs,  doped  n  type  with  Seat  about  2  x  10,K  cm  was 
grown.  The  thicknesses  of  the  AlAs  barriers  and  the  GaAs 
well  were  obtained  from  transmission  electron  microscopy 
measurements.  The  AlAs  layer  closer  to  the  substrate  was  52 
A  thick,  whereas  the  other  one  was  only  42  A  thick  The 
GaAs  well  between  the  two  AlAs  layers  was  62  A  wide. 
Although  the  barriers  were  doped  p  type,  they  were  fully 
depleted  and  still  acted  as  potential  energy  barriers  to  the 
electrons  in  the  GaAs  electrodes.  Devices  were  made  by  de¬ 
fining  mesas  on  the  epitaxial  sample  face  using  conventional 
photolithography  as  a  GaAs  etch  |4I.T, 
H;S04:H:0::H,0|.  Ohmic  contacts  were  made  on  the  sur¬ 
face  of  the  mesas  and  on  the  substrate  by  evaporating  Au-Ge 
and  Au,  and  annealing  at  410  °C.  The  samples  were  then 
either  mounted  on  transistor  headers  and  wire  bonded,  or 


FIG.  1  Current -voltage  curves  at  (at  300  K.  (bl  250  K.  and  ic)  77  K  for 
circular  devices  15  /im  in  diameter,  in  the  voltage  range  *  500  mV  The 
inflections  and  the  negative  resistance  regions  are  due  to  resonant  tunneltng 
thorugh  the  lowest  energy  level  in  the  GaAs  quantum  well 
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probed  with  a  whisker.  The  mesas  were  circular  and  1  5->r 
/im  in  diameter 

Current-voltage  curves  i/-F),  as  well  as  first  derivamo 
\dl  /dV  |  and  second  derivatives  I d  :I  /</F:tof  the  /- 1 'curses, 
were  measured  from  room  temperature  to4.2  K.  /-^'curves 
were  taken  using  a  Tektronix  577  curve  tracer  First  and 
second  derivatives  were  obtained  using  modulation  tech¬ 
niques,  with  typical  modulation  voltages  of  1-5  mV  peak  it! 
peak. 

Figures  1(a),  l(bi,and  l(clshowexperimental/-Fcur\es 
at  300  K,  250  K,  and  77  K,  respectively.  These  curves  were 
taken  on  1 5-/zm-diam  devices,  and  were  reproducible  for 
devices  of  the  same  area  across  the  whole  wafer.  Resonant 
tunneling  is  visible  in  the  /-  F curves  at  room  temperature  as 
inflections  [Fig.  1(a)].  These  become  more  pronounced  as  the 
temperature  is  reduced  and  yield  negative  differential  resis¬ 
tance  regions  at  temperatures  as  high  as  260  K  In  this  sam¬ 
ple,  the/-Fcurvesat  77  and  4.2  K  were  almost  identical.  The 
negative  resistance  regions  exhibit  features  [Fig  lie)]  which 
have  not  yet  been  identified  but  may  be  attributable  to  inelas¬ 
tic  excitations  of  electronic  states  of  impurities  by  the  tunnel¬ 
ing  electrons. "  Also  note  that  the  zero  bias  resistance  varies 
very  little  with  temperature,  while  the  current  at  voltages 
greater  than  the  resonance  peaks  decreases  rapidly  as  the 
temperature  is  lowered.  This  is  an  indication  that  most  of  the 
current  near  zero  bias  is  due  to  resonant  tunneling,  and  that, 
at  equilibrium,  the  Fermi  level  is  close  to  the  first  resonance 
in  the  GaAs  well. 

In  the  low-temperature  /-F curves  [Fig.  llci],  the  rath 
of  the  current  at  the  peak  to  the  current  at  the  valley \lp/I  1 11 
3.8  in  reverse  bias  and  3.4  in  forward  bias. 

Whereas  Fig.  I(a|— (cj  show  resonant  tunneling  effects 
associated  with  the  lowest  energy  state  in  the  GaAs  well,  an 
additional  structure  was  observed  at  larger  bias  voltages  m 
the  low-temperature  spectra.  Although  these  features  disap¬ 
peared  in  the  thermal  current  at  room  temperature,  thev 
could  be  identified  without  ambiguity  at  77  K.  Figure  2, 
which  is  a  4.2  K.  /-Fcurve.  exhibits  the  structure  at  approxi¬ 
mately  —  750  and  950  mV,  in  reverse  bias  and  forward  bias. 


v oil oqe  .mv  ' 


FIG  2  Current- voltage  curve  at  4  2  K  for  circular  devices  1 5 /4mm  dome 
ter.  in  the  voltage  range  ♦  2.5  V  Besides  the  negative  resistance  regions  due 
to  resonant  tunneling  through  Et\  ♦  200  mVi.  there  are  inflections  due  to 
resonant  tunneling  through  t\  the  second  energy  level  in  the  GaAs  *e|] 
i  7 so  mV  and  950  mVi,  and  maybe  also  through  £\.  the  third  energs  levei 
i  -  2  Vi 
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FIG  3  Second  derivative  spectrum  at  4.2  K  for  the  same  device  as  that  of 
Fig  2.  in  the  voltage  range  ±  l  5  V  The  very  large  peaks  correspond  to 
the  first  energy  level  in  the  GaAs  well,  and  the  small  peaks  to  £\,  the  second 
energy  state  in  the  well 

respectively  Although  the  amplitude  of  these  effects  with 
respect  to  the  background  current  is  much  smaller  than  for 
the  negative  resistances  discussed  previously,  they  do  corre¬ 
spond  to  tunneling  through  the  first  excited  state  in  the 
GaAs  well  These  effects  also  appear  in  the  derivative  spec¬ 
tra.  Figure  3  is  a  second  derivative)^  ZI  /dV  :)at  4.2  K  for  the 
same  device  as  that  shown  in  Fig.  2.  The  peaks  labeled  E: 
match  the  structure  identified  in  the  /-Fcurve  (Fig.  2|  as  the 
second  resonance. 

Calculations  based  on  the  Tsu  and  Esaki  model,1''  and 
taking  into  account  the  fact  that  part  of  the  voltage  drop 
occurs  in  the  GaAs,  give  reasonable  agreement  with  the  ex¬ 
perimental  data.  These  calculations  give  additional  evidence 
that  the  structure  observed  in  the  7-F curves  corresponds  to 
the  first  and  second  resonances  in  the  well.  Furthermore,  for 
a  given  GaAs  well,  the  tunneling  transmission  probability 
associated  with  a  resonance  decreases  rather  quickly  as  the 
AlAs  barriers  become  more  asymmetric.  This  could  explain 
the  weakness  of  the  second  resonance  with  respect  to  the  first 
one. 

It  should  be  mentioned  that,  although  the  structure  cor¬ 
responding  to  the  first  resonance  is  reproducible  across  the 
wafer  both  in  position  and  amplitude,  it  is  not  quite  so  for  the 
second  resonance.  This  might  be  attributable  to  fluctuations 
in  the  quality  of  the  contacts  at  the  high  currents  imposed  on 
the  devices  in  this  voltage  range,  and  to  gradual  fluctuations 
in  doping  and/or  layer  thicknesses  across  the  wafer 

Besides  the  two  resonances  discussed  above,  additional 
effects  appear  in  the  /-  F  curves  in  reverse  bias  at  about 
-  2  V  (Fig  2)  This  inflection  is  wide  and  its  amplitude  with 


respect  to  the  background  current  is  very  small.  Its  position 
would  be  consistent  with  the  location  of  the  third  resonance 
in  the  well.  Technical  limitations  have  prevented  a  more  de¬ 
tailed  study  of  this  feature. 

In  summary,  we  have  studied  current  transport  perpen¬ 
dicular  to  two  AlAs  barriers  separated  by  a  GaAs  quantum 
well  in  a  GaAs/ AlAs  heterostructure  grown  by  MOCVD. 
The  7-F  curves  and  their  derivatives  show  evidence  of  reso¬ 
nant  tunneling  through  the  structure.  Effects  corresponding 
to  the  first  resonance  in  the  GaAs  quantum  well  are  visible  in 
the  /-F curves  at  room  temperature,  and  negative  resistance 
regions  are  observed  at  temperature  as  high  as  260  K.  The 
low-temperature  7-  F curves  also  display  an  additional  struc¬ 
ture  corresponding  to  resonant  tunneling  through  the  sec¬ 
ond  and  perhaps  the  third  energy  levels  in  the  well. 
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Uniformity  of  current- voltage  characteristics  in  GaAs/AIAs  tunnel  structures  grown  by 
metalorganic  chemical  vapor  deposition  has  been  investigated  by  studying  electronic  transport 
perpendicular  to  GaAs  layers  separated  by  one  or  two  AlAs  barriers.  For  each  sample, 
measurements  of  the  current-voltage  curves  and  their  first  and  second  derivatives,  for  series  of  10 
to  80  identical  devices  taken  randomly  across  the  sample,  gave  reproducible  and  uniform  results. 

These  showed  evidence  that  the  doping  concentrations  and  layer  thicknesses  were  uniform  across 
the  wafers.  In  single  barrier  structures,  they  further  suggested  that  the  average  fluctuations  in  the 
thicknesses  of  ultrathin  layers  should  be  within  one  atomic  layer. 


Semiconductor  heterostructures  with  small  characteris¬ 
tic  dimensions  are  the  object  of  considerable  theoretical  and 
experimental  work.  They  are  frequently  fabricated  from 
GaAs  and  AI,  Ga,  _  ,  As  using  techniques  such  as  molecular 
beam  epitaxy  (MBE)  and  metalorganic  chemical  vapor  depo¬ 
sition  (MOCVD).  Since  the  conduction-band  offset  between 
the  GaAs  and  the  Al.Ga,  _  ,  As  causes  Al,  Ga,  ,  As  layers 
to  act  as  potential  energy  barriers  to  electrons  in  the  GaAs 
layers,  it  becomes  possible  to  grow  quantum  well  and  quan¬ 
tum  barrier  structures.  In  these  structures,  an  essential  cur¬ 
rent  transport  mechanism  is  tunneling.  The  tunneling  effects 
are  associated  with  very  important  theoretical  properties 
and  are  expected  to  yield  novel  electronic  features  of  consid¬ 
erable  interest.  In  the  simplest  configuration,  where  two 
GaAs  layers  are  separated  by  a  single  Al,  Ga,  _  „  As  barrier, 
elastic  and  inelastic  tunneling  take  place. 1-4  In  samples  con¬ 
taining  a  GaAs  quantum  well  adjacent  to  two  Al,  Ga,  ,  As 
barriers,  resonant  tunneling  occurs.'''10 

We  report  in  this  letter  a  systematic  study  of  the  repro¬ 
ducibility  and  the  uniformity  of  these  phenomena.  We  mea¬ 
sured  current- voltage  (/-  V )  curves  and  their  first  \dl  /dV)  and 
second  {d2I /dV2)  derivatives  for  MOCVD  grown  structures 
with  both  one  and  two  AlAs  barriers.  For  each  sample,  sev¬ 
eral  series  of  10  to  80  identical  devices  of  given  sizes,  taken  at 
random  across  the  sample,  were  systematically  tested  The 
results  of  the  measurements,  which  were  strongly  correlated 
to  the  parameters  of  the  structures,  showed  evidence  of  the 
growth  uniformity  of  the  wafers. 

The  samples  used  in  this  study  were  grown  by  an 
MOCVD  technique"  12  on  a  GaAs  substrate  doped  n  type 
with  Si  at  (2-3) X  10'*  cm'’.  The  first  layer  grown  on  the 
substrate  was  one  of  GaAs,  2-3  /im  thick,  doped  n  type  with 
Se  at  about  2  X  10IK  cm  “  \  Following  this,  either  a  single  thin 
layer  of  AlAs,  or  two  thin  layers  of  AlAs  separated  by  a  thin 
layer  of  nominally  undoped  GaAs,  were  grown.  Finally,  a 
top  layer  of  GaAs,  doped  n  type  with  Se,  was  grown  The 
AlAs  layers  were  doped p  type  with  Mg  (Refs  13  and  14)  at 
approximately  1  x  10IH  cm  ‘  ’  In  the  single  barrier  samples, 
the  AlAs  layer  was  typically  from  50  to  1 50  A  wide,  whereas 
the  top  GaAs  layer  was  approximately  ipm  thick  and  doped 
at  about  5x  10IK  cm  V  In  the  double  barrier  samples,  the 
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AlAs  layers  and  the  GaAs  well  were  on  the  order  of  50  A 
wide,  whereas  the  top  GaAs  layer  was  approximately  2 pm 
thick  and  doped  at  about  2  X  1018  cm-1.  Although  the  bar¬ 
riers  were  doped  p  type,  they  were  fully  depleted  and  still 
acted  as  potential  energy  barriers  to  the  electrons  in  the 
GaAs  electrodes.  Devices  were  prepared  on  the  samples  by 
defining  mesas  on  the  epitaxial  surface  using  con¬ 
ventional  photolithography  and  a  GaAs  etch  (4:1:1, 
HjS04:H202:H30).  Ohmic  contacts  were  made  on  the  sur¬ 
face  of  the  mesas  and  on  the  substrate  by  evaporating  Au-Ge 
and  Au,  and  annealing  at  400-420 ’C  for  15-20  s.  Finally, 
the  samples  were  either  mounted  on  transistor  headers  and 
wire  bonded,  or  probed  with  a  whisker.  The  mesas  were  cir¬ 
cular  and  12-700  pm  in  diameter.  For  the  purpose  of  this 
study,  special  care  was  taken  to  prepare  all  the  samples  and 
devices  under  similar  conditions.  In  particular,  since  many 
electrical  properties  depend  on  the  device  area,  the  actual 
device  diameters  were  carefully  measured.  Furthermore,  in 
order  to  reduce  the  error  introduced  by  undercutting  during 
the  etching  process,  the  top  GaAs  layer  was  preliminary 
etched  down  to  about  0.2  pm.  The  mesas  were  then  obtained 
by  etching  only  about  1  pm  down  from  the  surface.  This 
procedure  significantly  reduced  the  uncertainty  in  the  device 
areas  and  yielded  well  defined  mesas. 

Starting  from  wafers  of  about  2  cm2,  at  least  five  small 
pieces  of  approximately  4  mm2,  chosen  as  far  apart  from 
each  other  as  possible,  were  cleaved.  A  large  number  of  me¬ 
sas  were  then  defined  on  each  one  of  these  pieces  and,  de¬ 
pending  on  their  size,  about  10  to  80  devices  taken  at  random 
were  systematically  tested.  Although  this  procedure  was  re¬ 
peated  for  different  device  sizes  ranging  between  12  and  700 
£im,  a  more  thorough  study  was  performed  for  devices  close 
to  50  and  80^m  in  diameter.  Current-voltage  curves,  as  well 
as  their  first  and  second  derivatives,  were  measured  at  tem¬ 
peratures  ranging  from  4.2  to  300  K.  The  I-V  curves  gave 
information  about  current  transport  mechanisms.  They  re¬ 
vealed  mainly  elastic  tunneling  in  the  single  barrier  struc¬ 
tures  and  resonant  tunneling  in  the  double  barrier  samples. 
Structure  in  the  derivatives  of  the  l-V  curves  indicated  the 
presence  of  inelastic  and  resonant  tunneling  processes  and 
allowed  the  identification  of  fundamental  excitations  creat- 
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Device  Diameter  (fim) 

FIG.  1.  Lof ,«  (zero-bias  resistance)  vs  lof  10  (device  diameter)  for  a  single¬ 
harrier  sample. 

ed  by  the  tunneling  electrons.  The  I-V  curves  were  taken 
using  a  Tektronix  577  curve  tracer.  Point-by-point  measure¬ 
ments  using  a  dc  power  supply  were  also  made.  First  and 
second  derivatives  were  obtained  using  modulation  tech¬ 
niques,  with  typical  modulation  voltages  of  1-5  mV  peak  to 
peak.13 

Some  of  the  properties  which  were  found  to  be  of  inter¬ 
est  are  described  below.  For  the  one-barrier  samples,  in 
which  elastic  and  inelastic  tunneling  occur,  we  studied  (i)  the 
zero-bias  resistance  at  different  temperatures  and  (ii)  the  uni¬ 
formity  and  reproducibility  of  the  I-V  curves  as  well  as  their 
first  and  second  derivatives.  In  the  double-barrier  structures, 
resonant  tunneling  takes  place  and  the  I-V  curves  feature 
regions  of  negative  differential  resistance.  For  these  samples, 
we  performed  a  careful  study  of  (i)  the  zero-bias  resistance  at 
different  temperatures,  (ii)  the  voltages  and  currents  at  the 
peaks  and  valleys  of  the  negative  resistance  regions,  and  (iii) 
the  uniformity  and  reproducibility  of  the  features  displayed 
in  the  derivatives.  Only  a  few  of  these  properties  are  dis¬ 
cussed  in  the  present  letter. 

Figure  1  is  a  plot  of  the  logarithm  of  the  zero-bias  resis¬ 
tance  as  a  function  of  the  logarithm  of  the  device  diameter 
for  a  single-barrier  sample.  The  error  bars  correspond  to 
averages  of  more  than  20  measurements  made  on  different 
devices  of  a  given  size  and  taken  at  random.  The  crosses 
correspond  to  single  measurements.  The  curve  passing 
through  the  experimental  data  is  a  straight  line  of  slope  al¬ 
most  exactly  (  —  2),  indicating  that  the  zero-bias  resistance 
does  scale  with  area.  The  zero-bias  resistance,  which  de¬ 
pends  on  the  transmission  rate  through  the  AlAs  barrier  for 
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FIG.  2.  Second  derivative  spectra  at  4.2  K  for  seven  devices  7S/jm  in  diam¬ 
eter,  selected  randomly  among  100  of  the  same  size,  for  a  single  barrier 
structure. 

elastic  tunneling,  should  decrease  exponentially  as  the  width 
of  the  barrier  is  increased.  Therefore,  it  is  expected  to  depend 
strongly  on  the  thickness  of  the  AlAs  barrier  through  which 
the  electrons  actually  tunnel.  Since  the  tunneling  electrons 
will  always  tend  to  pass  through  the  thinnest  part  of  a  bar¬ 
rier,  variations  in  the  values  of  the  zero-bias  resistance  are 
directly  correlated  to  the  fluctuations  in  the  barrier  width. 
Furthermore,  keeping  all  the  other  parameters  constant,  a  7 
A  change  in  the  barrier  thickness  makes  the  zero-bias  resis¬ 
tance  vary  by  an  order  of  magnitude.13  Consequently,  tun¬ 
neling  measurements  provide  a  sensitive  test  of  barrier 
thickness.  In  addition,  the  zero-bias  resistance  depends  sig¬ 
nificantly  on  the  doping  concentrations,  which  define  the 
actual  position  of  the  Fermi  level  in  the  GaAs  electrodes.  It 
should  also  depend  on  the  contacts.  However,  these  had  re¬ 
sistances  which  were  negligible  compared  to  the  tunneling 
resistances  of  the  barriers.  Therefore,  they  did  not  alter  sig¬ 
nificantly  the  tunneling  currents.  In  Fig.  1,  it  should  be  noted 
that  the  error  bars,  which  represent  the  actual  standard  devi¬ 
ation  of  the  results,  are  very  small.  This  clearly  indicates  that 
the  results  of  the  measurements  were  reproducible  across  the 
whole  wafer.  This  further  shows  that  leakage  and  edge  cur¬ 
rents,  as  well  as  effects  associated  with  the  local  fluctuations 
in  the  electrostatic  potential  due  to  the  ionized  dopants,  were 
negligible  compared  to  the  tunneling  currents.  The  above 
discussion  about  the  dependence  of  the  zero-bias  resistance 
upon  barrier  thickness  and  doping  concentrations  allows  us 
to  conclude  that  both  the  doping  profiles  and  the  layer  thick¬ 
nesses  were  uniform  across  the  entire  sample.  It  further  sug¬ 
gests  that  the  average  fluctuations  in  the  thicknesses  of  ul- 
trathin  layers  in  single-barrier  structures  should  be  within 
one  atomic  layer. 

Figure  2  is  the  superposition  of  seven  second  derivative 
\d 2 1  /dV2)  curves  at  4.2  K  for  seven  different  devices,  78  fim 
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FIG.  3.  /-Kcurvesat  300  K  and  77  K  for  ten  devices,  78  /jm  in  diameter,  for 
a  double-barrier  sample. 

in  diameter,  taken  randomly  among  100  of  the  same  size. 
The  sample  is  the  same  one  as  in  Fig.  1.  The  fact  that  these 
spectra  overlap  exactly  was  a  common  feature  of  all  the  se¬ 
ries  of  devices  tested.  Second  derivative  curves  reveal  inelas¬ 
tic  tunneling  through  the  AIAs  barrier.  Tunneling  electrons 
can  create  fundamental  excitations,  provided  the  applied 
bias  is  increased  above  the  excitation  energy  (Au0).  Thus, 
there  exists  a  threshold  voltage  K0  =  ftcog/e,  which  corre¬ 
sponds  to  a  change  in  the  slope  of  the  I-V  curve.  Since  the 
structure  is  symmetric,  the  same  changes  should  occur  in 
forward  and  in  reverse  biases.  Therefore,  steps  which  are 
symmetric  with  respect  to  bias  appear  in  the  first  derivative 
spectra  at  voltages  ±  Vg,  whereas  the  second  derivative 
curves  show  peaks  at  the  same  voltages,  but  antisymmetric 
with  respect  to  bias.  In  Fig.  2,  a  reproducible  structure  is 
observed  at  the  GaAs  and  the  AIAs  longitudinal  optical  (LO) 
phonon  energies,  ±  36  meV  and  ±  50  meV,  respectively. 

For  double-barrier  samples,  series  of  more  than  60  de¬ 
vices  were  tested  for  devices  close  to  20,  50,  and  80  fxm  in 
diameter.  I-V  curves  as  well  as  first  and  second  derivatives 
were  measured.  In  a  given  series,  the  curves  would  all  coin¬ 
cide,  except  for  an  occasional  one,  which  could  be  due  to  the 
preparation  procedure,  a  bad  contact,  or  an  improper  bond. 
Figure  3  shows  experimental  I-  V  curves  at  300  and  77  K  for 
10  different  devices,  78  pm  in  diameter,  taken  randomly 
across  the  wafer  of  a  double-barrier  sample.  The  breaks  in 
the  low-temperature  I-V  curves  are  due  to  the  fact  that  the 
curve  tracer  used  for  these  measurements  could  not  reliably 
resolve  regions  of  negative  differential  resistance.  But  a  more 


detailed  discussion  of  the  negative  resistance  regions  is  not 
the  subject  of  this  letter  and  was  reported  elsewhere.10  Reso¬ 
nant  tunneling  currents  are  critically  dependent  upon  the 
layer  thicknesses  and  the  doping  concentrations.  Therefore, 
the  reproducibility  and  uniformity  of  the  experimental  data 
are  again  a  clear  indication  that  the  geometrical  parameters 
and  the  doping  profiles  were  uniform  throughout  the  entire 
sample. 

In  summary,  we  have  studied  the  reproducibility  and 
uniformity  of  I-V  characteristics  and  their  derivatives  in 
MOCVD  grown  GaAs/ AIAs  tunnel  structures.  We  have  in¬ 
vestigated  electronic  transport,  mainly  elastic,  inelastic, and 
resonant  tunneling,  perpendicular  to  GaAs  electrodes  sepa¬ 
rated  by  one  or  two  AIAs  layers.  For  each  sample,  systema¬ 
tic  measurements  were  made  on  a  large  number  of  devices 
taken  randomly  across  the  wafer.  The  measured  properties 
were  critically  dependent  on  the  layer  thicknesses  and  the 
doping  concentrations.  This  study  did  not  allow  us  to  draw 
conclusions  about  the  abruptness  of  the  interfaces  and  the 
presence  of  defects  and  impurities.  However,  the  reproduc¬ 
ibility  and  uniformity  of  the  experimental  results  showed 
evidence  that  the  doping  profiles  and  the  geometrical  param¬ 
eters  of  the  samples  were  uniform  across  the  wafers.  In  the 
single  barrier  samples,  they  further  suggested  that  the  aver¬ 
age  fluctuations  in  the  thicknesses  of  ultrathin  layers  should 
be  within  one  atomic  layer. 
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Elastic  and  inelastic  tunneling  processes  are  investigated  in  GaAs-AlAs-GaAs  double 
heterojunctions  grown  in  the  [100]  direction  by  metal-organic  chemical  vapor  deposition 
(MOCVD).  The  AlAs  quantum  barriers  in  the  heterostructures  studied  are  doped  p-type  with 
Mg.  Theoretical  calculations  of  tunneling  currents  are  performed  and  compared  with 
experimental  I-V  data.  It  is  found  that  for  structures  with  thin  AlAs  barriers,  the  dominant 
current  transport  mechanism  at  low  temperatures  is  tunneling  through  the  AlAs  band  gap  at  both 
the  T  and  X  points.  This  is  consistent  with  inelastic  processes  observable  in  first  (dl /dV)  and 
second  (d2I /dV *)  derivative  spectra  obtained  with  modulation  techniques.  A  simple  model, 
developed  for  calculating  impurity-assisted  tunneling  currents,  shows  that  the  role  of  barrier 
impurities  becomes  more  important  as  the  barrier  is  grown  thicker.  Implications  of  some  of  these 
results  for  resonant  tunneling  heterostructures  consisting  of  two  AlAs  quantum  barriers 
separated  by  a  GaAs  quantum  well  are  discussed.  Experimental  second  derivative  spectra 
showing  reproducible  features  are  also  presented  for  these  double  barrier  structures. 


I.  INTRODUCTION 

Because  of  their  numerous  potential  device  applications,  se¬ 
miconductor  quantum  well  and  quantum  barrier  hetero- 
structures  are  the  object  of  considerable  theoretical  and  ex¬ 
perimental  work.  These  structures  are  frequently  fabricated 
from  GaAs  and  Al*  Ga,  _  *  As,  using  techniques  such  as  mo¬ 
lecular  beam  epitaxy  (MBE)  and  metal-organic  chemical 
vapor  deposition  (MOCVD)  to  achieve  small  characteristic 
dimensions.  In  structures  with  thin  quantum  barriers,  an 
important  current  transport  mechanism  is  electron  tunnel¬ 
ing.  Several  studies  of  electronic  transport  in  single  and  dou¬ 
ble  barrier  tunnel  structures  have  been  reported. 110  The  as¬ 
sociated  tunneling  efects  yield  novel  electronic  properties, 
some  of  which  remain  to  be  fully  understood.  Most  of  the 
studies  to  date,  both  theoretical  and  experimental,  have  used 
AljGa,  _ ,  As  quantum  barriers  in  which  the  A1  composi¬ 
tion,  x,  was  small  enough  that  the  alloy  was  still  direct.5'10 
Furthermore,  most  of  the  experimental  investigations  have 
been  performed  on  MBE  grown  wafers.4-8  Only  a  few  stud¬ 
ies  have  been  reported  on  MOCVD  grown  samples, 1-3 
and/or  structures  in  which  the  barriers  consisted  of  indirect 
Al,Ga,  _„As,5  4  or  pure  AlAs. 1,2 

In  this  paper,  we  concentrate  primarily  on  current  trans¬ 
port  mechanisms  at  low  temperatures  in  GaAs-AlAs-GaAs 
double  heterojunctions  grown  by  MOCVD  in  the  [100]  di¬ 
rection.  In  Sec.  II,  we  briefly  present  the  sample  growth  and 
preparation  procedures.  In  Sec.  Ill,  we  report  both  experi¬ 
mental  and  theoretical  results  for  structures  with  thin  p-type 
AlAs  barriers.  We  first  calculate  current-voltage  (I-V) 
characteristics  at  low  temperatures  due  to  elastic  tunneling 
through  the  AlAs  band  gap  at  the  T  point,  using  the  actual 
shape  of  the  barrier  obtained  from  band  bending  calcula¬ 
tions.  Upon  comparing  the  calculated  and  experimental  I-  V 
curves,  it  is  found  that  tunneling  through  the  AlAs  f -point 


barrier  fails  to  explain  ihe  experimental  data.  This  indicates 
that  other  transport  mechanisms  must  contribute  to  the  total 
current.  Tunneling  through  the  AlAs  band  gap  at  the  X 
point  is  then  discussed.  Good  theoretical  fits  of  the  experi¬ 
mental  I-V  data  can  be  obtained  by  considering  tunneling 
through  both  the  T-  and  T-point  barriers.  In  Sec.  IV,  we 
introduce  a  new  element  to  the  calculation  of  the  theoretical 
I-V  curves:  the  presence  of  trap  levels  in  the  AlAs  band  gap, 
due  to  defects  or  impurities.  A  simple  model  is  developed  to 
calculate  the  impurity-assisted  tunneling  current.  This  cur¬ 
rent  becomes  more  important  as  the  quantum  barriers  be¬ 
come  thicker,  and  can  actually  become  the  dominant  current 
transport  mechanism  through  thick  enough  barriers.  In  Sec. 
V,  implications  of  some  of  these  results  for  resonant  tunnel¬ 
ing  heterostructures  made  of  two  AlAs  quantum  barriers 
separated  by  a  GaAs  quantum  well  are  discussed.  For  these 
structures,  experimental  derivative  spectra  showing  repro¬ 
ducible  and  periodic  features  are  also  reported.  Finally,  the 
results  of  the  present  study  are  summarized  in  Sec.  VI. 


II.  EXPERIMENTAL  PROCEDURES 


The  samples  used  in  this 
MOCVD  technique”12  on  a 


study  were  grown  by  an 
[  100] -oriented  GaAs  sub¬ 


strate  doped  n-type  with  Si  at  2-3  X  1018  cm-'.  A  first  epi¬ 
taxial  layer  of  degenerate  GaAs,  doped  n-type  with  Se,  was 
grown  2-3  gm  thick.  In  single  barrier  structures,  this  layer 
was  followed  first  by  an  AlAs  layer  doped  p- type  with 
Mg, 15,14  and  then  by  another  GaAs  layer,  degenerately 
doped  with  Se.  Samples  were  grown  with  GaAs  dopings 
ranging  between  lxlO18  andSxlO18  cm-5,  AlAs  barrier 
thicknesses  between  40  and  250  A,  and  barrier  dopings  esti¬ 
mated  to  vary  between  I  X  IO17  and  3  X  1018  cm-5.  In  dou- 
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ble  barrier  structures,  the  first  epitaxial  GaAs  layer  was  fol¬ 
lowed  by  two  thin  AlAs  layers  separated  by  a  quantum  well 
of  nominally  undoped  GaAs.  The  AlAs  barriers  were  doped 
p-type  with  Mg  at  approximately  1  ,<  10ls  cm  Finally,  a 
GaAs  top  layer  was  grown,  degenerately  doped  with  Se.  The 
thicknesses  of  the  AlAs  barriers  and  the  GaAs  wells  were 
determined  by  transmission  electron  microscopy  (TEM). 
Electrode  dopings  were  obtained  from  Polaron  doping  pro¬ 
files.  Barrier  dopings  were  estimated  from  the  flow  rates 
used  during  growth.  They  could  thus  be  in  error  due  to  the 
difficulty  of  accurately  calibrating  dopings  in  such  thin  lay¬ 
ers. 

Devices  were  made  by  defining  mesas  on  the  epitaxial 
sample  face  using  conventional  photolithography  and  a 
GaAs  etch  (4:1:1,  H;S04  :H- O-  :H,  O).  The  mesas  were  cir¬ 
cular  and  50-700  pm  in  diameter.  Ohmic  contacts  were 
made  on  the  surface  of  the  mesas  and  on  the  substrate  by 
evaporating  Au-Ge,  or  Au-Ge,  Ni,  Au,  and  then  annealing 
at  380-410  °C  for  20-30  s. 

Measurements  of  current-voltage  curves  (I-V),  as  well  as 
first'  (dl  /dV)  and  second  (d1 1  /dV1 )  derivativesof  the/-F 
curves  were  performed  at  temperatures  ranging  from  300  to 
4.2  K.  The/-F curves  were  measured  with  an  HP  4145  semi¬ 
conductor  parameter  analyzer.  The  first  and  second  deriva¬ 
tive  spectra  were  obtained  using  modulation  techniques  at  5 
and  50  kHz,  respectively. 

III.  RESULTS  ON  THIN  SINGLE  BARRIER 
STRUCTURES 

In  this  section,  we  present  both  experimental  and  theoreti¬ 
cal  results  for  single  barrier  heterostructures  having  a  thin 
AlAs  barrier  layer  totally  depleted  of  carriers.  A  number  of 
studies  have  reported  that  the  dominant  current  transport 
mechanism  at  low  temperatures  in  Al.Ga,  _  ,  As  quantum 
barrier  structures  is  elastic  tunneling  through  the 
Al^Ga,  As  direct  band  gap  at  the  T  point.3'5'6  In  these 
studies,  theoretical  I-V  curves  were  calculated  assuming  a 
trapezoidal  shape  for  the  quantum  barrier  under  applied 
bias,  and  compared  to  experimental  data  on  logarithmic 
scales.  Nevertheless,  discrepancies  were  observed3  which  in¬ 
dicated  that  other  current  transport  mechanisms  can  com¬ 
pete  with  elastic  T -point  tunneling.  Theoretical  studies  have 
also  reported  that  T -point  electrons  in  the  GaAs  electrodes 
should  tunnel  through  the  Al,  Ga,  _  ,  As  band  gap  at  the  T 
point  even  in  indirect  alloys.15'1*’  However,  we  show  here 
that  experimental  data  for  structures  with  pure  AlAs  bar¬ 
riers  cannot  be  explained  by  T-point  tunneling  alone.  We 
propose  that  there  is  a  small,  but  nonnegligible  probability 
that  the  incident  electrons  tunnel  through  the  AlAs  .T-point 
barrier. 

Current  densities  are  calculated  using  the  following  ap¬ 
proach  First,  the  actual  shape  of  the  conduction  band  edge 
in  the  electrodes  and  the  barrier  is  determined  from  band 
bending  calculations.  These  calculations  are  found  to  be  of 
critical  importance,  particularly  in  samples  with  heavily 
doped  barriers  They  are  performed  by  solving  Poisson's 
equation  self-consistently  throughout  the  entire  heterostruc¬ 
ture  at  each  applied  voltage  The  transmission  coefficient  for 


tunneling  electrons  is  then  calculated  using  the  WKB  ap¬ 
proximation.  The  attenuation  constant  in  the  barrier  is  de¬ 
termined  from  a  two-band  model,  k-p  theory  formula'7  for 
electrons  tunneling  through  the  AlAs  T  point,  a^d  a  simple 
"one-band  model”  formula  for  electrons  tunneling  through 
the  AlAs  X point.  Finally,  the  current  density,  J,  is  calculat¬ 
ed  as  a  function  of  applied  voltage,  V,  using  the  approach  of 
Tsu  and  Esaki.9 

Although  our  study  was  performed  on  a  wide  selection  of 
samples  and  at  temperatures  ranging  from  300  to  4.2  K,  we 
only  illustrate  the  most  important  results  by  discussing  data 
obtained  at  low  temperatures  from  a  single  barrier  structure 
having  a  thin  p- type  AlAs  barrier  layer.  A  more  detailed 
report,  including  results  for  heterostructures  with  /i-type 
AlAs  barriers,  will  be  presented  in  a  forthcoming  publica¬ 
tion. 

In  the  structure  of  interest,  the  GaAs  electrodes  are  doped 
n-type  with  Se  at  4—5  x  1018  cm-3.  The  doping  in  the p-type 
AlAs  barrier  layer  is  estimated  to  be  on  the  order  of  10'* 
cm"1.  A  barrier  thickness  of  48  A  was  obtained  from  TEM 
measurements,  accurate  to  within  one  monolayer.  Figure  1 
shows  the  calculated  conduction  band  edge  of  the  structure 
at  the  T  and  X  points  as  a  function  of  distance  in  the  direction 
perpendicular  to  the  plane  of  the  layers.  Figure  1(a)  depicts 
thermal  electronic  equilibrium  while  Fig.  1  (b)  corresponds 
to  an  applied  bias  of  200  mV.  In  the  band  bending  calcula¬ 
tions.  the  dopings  in  the  n-type  GaAs  electrodes  and  the  p- 
type  AlAs  barrier  are  taken  to  be  4x  10’8  and  3x  10,v 
cm'!,  respectively.  The  barrier  thickness  is  48  A.  The  caicu 
lations  are  performed  assuming  a  temperature  of  0  K.  The 
values  of  the  band  offsets  are  key  parameters  in  determining 
the  shape  of  the  barriers  seen  by  the  tunneling  electrons. 
These  values  are  still  a  matter  of  debate.  In  a  recent  experi¬ 
mental  study,  Batey  e/a/.18  obtained  a  valence  band  discon¬ 
tinuity  of  0.55  eV  in  GaAs-AlAs  heterostructures.  Using 
this  result,  values  close  to  190meV  and  1.04  eV  can  be  calcu¬ 
lated  for  the  conduction  band  offsets  at  the  X  and  T  points, 
respectively.  It  should  be  noted  that  when  such  values  are 
used,  the  A'-point  conduction  band  edge  lies  at  higher  ener¬ 
gies  in  GaAs  than  in  AlAs.  This  is  illustrated  in  Fig.  1,  with 
the  solid  line  corresponding  to  the  T -point  conduction  band 
edge  and  the  dashed  line  to  the  A'-point  conduction  band 
edge.  Figure  1(a)  shows  that  the  barrier  is  totally  depleted  of 
carriers  and  that  the  effect  of  band  bending  is  to  increase  the 
average  barrier  height  seen  by  the  tunneling  electrons.  This 
effect  increases  with  barrier  doping  and  thickness.  In  Fig 
1(b),  a  voltage  of  200  mV  is  applied  across  the  structure. 
This  voltage  alters  the  shape  of  the  conduction  band  edge  in 
the  cladding  layers  as  well  as  in  the  barrier.  This  is  another 
essential  factor  in  calculating  tunneling  currents.  In  Fig.  2. 
we  show  the  experimental  characteristic,  J„p  (  V  ),  at  4.2  K. 
in  the  voltage  range  0-200  mV  and,  on  the  same  linear  scale, 
the  calculated  current  density,  (  V  ),  for  elastic  T -point 
tunneling.  Since  the  structure  is  symmetric,  the  J-  V  curves 
are  symmetric  with  respect  to  the  origin  and  we  only  show 
the  forward  bias  direction  The  two  curves  are  in  good  agree¬ 
ment  at  very  low  biases  (  s 20  mV ),  but  start  to  deviate  sig. 
nificantly  in  magnitude  and  shape  at  higher  voltages  While 
J ,r,  variesalmost  linearly  with  voltage,  theexpenment.il  cur- 
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Flo  1  Calculated  conduction  band  edge  at  the  T  point  (solid  line)  and  at 
the  X  point  (dashed  line),  as  a  function  of  distance  in  the  direction  perpen¬ 
dicular  to  the  plane  of  the  layers,  for  a  GaAs-AlAs-GaAs  single  barrier 
heterostructure  The  GaAs  electrodes  are  doped  e-type  at  4  -  JO'*  cm 
The  AlAs  barrier  is  48  A  thick  and  doped  p-type  at  3  ■  10"  cm  The 
conduction  band  discontinuities  at  the.T and  T  points  are  IWmcV  and  I  04 
eV,  respectively  (a)  Corresponds  to  thermal  electronic  equilibrium,  and 
( b)  to  an  applied  voltage  of  200  mV  The  solid-dashed  lines  correspond  to 
the  Fermi  level.  E, .  in  the  electrodes  The  temperature  is  assumed  to  beO  K 


rent  density,  displays  a  nearly  exponential  voltage  de¬ 
pendence.  Similar  discrepancies  between  experimental  and 
theoretical  data  are  observed  in  all  of  our  samples.  This  is  a 
clear  indication  that  elastic  tunneling  through  the  AlAs 
band  gap  at  the  f  point  does  not  account  for  all  of  the  experi¬ 
mental  current  and  that  other  transport  mechanisms  are 
present  These  observations  are  further  supported  b>  deriv  a- 
tne  spectra  which  reveal  threshold  voltages  for  impuritv- 
assisted  and/or  inelastic  tunneling  processes  This  is  illus¬ 
trated  in  the  second  derivative  curve.  1  d  /  </(  jil 
shown  in  Fig  .’  This  spectrum  displays  two  prominent  and 
broad  features,  labeled  i  a  )  and  (  b  ),  peaked  at  approximate¬ 
ly  20  and  80  mV  It  also  reveals  two  narrower  peaks  at  about 


J-V  Curves 


Hi,  2  Experimental  and  theoretical  current- voltage  characteristics  for  the 
single  barrier  structure  shown  in  Fig  I  The  solid  line  corresponds  to  the 
experimental  current  density.  7,, r  The  dashed  line  is  the  calculated  T -point 
elastic  current  density.  J[, 


37  and  50  mV  These  have  been  attributed  previously  to  elec¬ 
tron-optical  phonon  coupling  in  the  GaAs  electrodes,  and 
the  excitation  of  longitudinal  optical  (LO)  phonons  in  the 
AlAs  barnei,  respectively.'  Leakage  has  been  eliminated  as 
a  possible  source  of  current  by  verifying  that  the  experimen¬ 
tal  J-  V  curves  were  reproducible  and  independent  of  device 
area.  The  behavior  of  the  T -point  elastic  current  can  easily 
be  understood  in  terms  of  the  complex  band  structure  of 
AlAs:  In  the  band  gap,  the  imaginary  part  of  the  complex 
wave  vector  at  the  T  point  vanes  very  slowly  with  energy, 
except  near  the  conduction  and  valence  band  edges.  This  is 
due  to  the  small  T -point  effective  mass.  Therefore,  because 
of  the  large  conduction  band  offset  (  ~  1  eV),  the  T-point 
current  remains  slowly  varying  with  voltage  up  to  about  O  S 
V. 

In  order  to  tentatively  explain  the  discrepancies  between 
experimental  J-  V  data  and  calculated  T -point  elastic  tunnel¬ 
ing  currents,  a  number  of  transport  mechanisms  can  be  sug¬ 
gested.  The  simplest  is  the  emission  of  phonons  or  the  excita¬ 
tion  of  collective  modes  by  electrons  tunneling  through  the 
AlAs  T -point  barrier.  Although  some  of  these  inelastic  pro¬ 
cesses  may  be  resolved  in  derivative  spectra,  they  yield  cur¬ 
rents  which  are  typically  at  least  two  or  three  orders  of  mag¬ 
nitude  smaller  than  elastic  currents,  and  therefore  too  small 
to  account  for  the  observed  discrepancies.  Another  possible 
current  transport  mechanism  is  tunneling  via  energy  levels 
in  the  AlAs  hand  gap  or  at  the  GjAs/AIAs  interfaces  A 
simple  model  for  impurity-assisted  tunneling  is  developed  in 
Sec  IV  Assuming  that  no  impurity  bands  or  clusters  are 
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Fig  3  Experimemal  second  derivative  spectrum  id'I  /'JV:  )t  1  ),  at  4.2  K 
for  the  single  barrier  heterostructure  illustrated  in  Fig  I  The  two  most 
prominent  peaks  labeled  (a)  and  (b)  may  correspond  to  inelastic  processes 
by  which  GaAs  F -point  electrons  scatter  into  Al>\s  A'-point  states  The  two 
narrower  peaks  at  31  and  50  mV  have  been  attributed  previously  to  elec¬ 
tron-optical  phonon  coupling  in  the  GaAs  electrodes,  and  the  excitation  of 
longitudinal  optical  ( LO)  phonons  in  the  AlAs  barrier,  respectively  ( Ref. 
1) 


formed,  this  mechanism  gives  negligible  current  densities  for 
the  48  A  barrier  structure  considered  here.  Another  possibil¬ 
ity  is  to  take  into  account  tunneling  through  the  AlAs  band 
gap  at  symmetry  points  other  than  the  T  point.  In  the  case  of 
epitaxial  growth  in  the  [  100]  direction,  an  important  contri¬ 
bution  could  come  from  the  AlAs  A'-point  valleys.  For  elec¬ 
trons  tunneling  from  left  to  right,  new  channels  open  when 
the  Fermi  level  in  the  left  electrode  and  the  A-point  conduc¬ 
tion  band  edge  in  the  right  electrode  align.  However,  such 
processes  can  only  occur  at  voltages  >350  meV  in  the  struc¬ 
ture  discussed  here.  This  is  again  in  disagreement  with  ex¬ 
perimental  data.  Other  mechanisms  can  allow  electrons  oc¬ 
cupying  T -point  states  in  both  GaAs  electrodes  to  tunnel 
through  the  AlAs  A'-point  barrier.  Two  such  processes  are 
discussed  below. 

First,  T -point  electrons  near  the  GaAs-AlAs  interface 
may  scatter  inelastically  into  virtual  states  beneath  the  four 
AlAs  X  valleys  lying  along  the  kr  and  k2  directions,  parallel 
to  the  plane  of  the  layers.  They  may  then  tunnel  through  the 
AlAs  band  gap  at  the  X  point,  and  scatter  again,  near  the 
AlAs-GaAs  interface,  into  the  T  valley  in  the  right  elec¬ 
trode  Since  electrons  travel  along  the  x  axis,  perpendicular 
to  the  plane  of  the  layers,  the  effective  mass  in  the  barrier  is 
the  small  transverse  effective  mass  m*  =  0  10m  ,  where  m, 
is  the  free  electron  mass  Such  inelastic  processes,  which  re¬ 
quire  two  scattering  events,  can  occur  through  a  number  of 
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mechanisms  such  as  phonon  scattering,  for  example.  Calcu¬ 
lations  of  scattering  matrix  elements  and  current  densities 
are  being  performed  to  determine  which  of  these  processes 
may  explain  the  experimental  data.  The  matrix  elements  can 
be  complicated  functions  of  the  scattering  mechanism,  ap¬ 
plied  voltage,  structure  characteristics,  electron  energy,  and 
temperature.  In  the  present  study,  we  simply  assume  that 
there  is  a  small  probability  that  incident  electrons  tunnel 
through  the  AlAs  band  gap  at  the  X  point  and  that  the  scat¬ 
tering  matrix  element  for  a  given  process  is  constant.  The 
calculated  current  density  for  any  such  mechanism  thus  con¬ 
tains  an  empirical  factor,  Bm .  Threshold  voltages  for  these 
inelastic  processes  are  chosen  to  coincide  with  some  of  the 
most  prominent  peaks  observ  ed  in  experimental  second  deri¬ 
vative  spectra  (e/’Z/c/k0  )  f  /'  )  The  total  theoretical  current 
density,  /,h  (  V  ),  is  then  calculated  by  summing  the  contri¬ 
butions  from  T -point  elastic  tunneling  and  A'-point  inelastic 
tunneling.  Figure  4  is  a  plot  of  J  th  as  a  function  of  applied 
voltage  for  the  structure  discussed  here.  Only  two  inelastic 
processes  are  used  to  calculate  the  A'-point  inelastic  tunnel¬ 
ing  current.  These  two  processes  are  assigned  threshold  vol¬ 
tages  of  20  and  80  mV,  and  empirical  Z?,n’sof0.6x  10-4  and 
1.3X  10~4,  respectively.  The  threshold  voltages  correspond 
to  the  two  most  prominent  peaks  labeled  (a)  and  (b)  in  the 
second  derivative  curve  illustrated  in  Fig.  3.  Since  two  scat¬ 
tering  events  are  required  to  allow  electrons  to  tunnel 
through  the  AlAs  A'-point  barrier,  peaks  should  be  observed 


Voltage  (mV) 

Fig  4  Total  theoretical  current  density,  /lh,  obtained  by  summing  the 
contributions  of  T -point  elastic  tunneling  and  two  ,T-pomt  inelastic  tunnel¬ 
ing  mechanisms,  for  the  same  single  barrier  structure  as  in  Fig  2.  For  the 
inelastic  processes,  threshold  voltages  of  20  and  80  mV  are  taken  tn  agree¬ 
ment  with  experimental  second  derivative  (d2 1  /dVl  )  data  The  corre¬ 
sponding  empirical  B  s  are  0.6  x  10' 4  and  1.3x  10  ~4,  respectively 
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in  the  (d2 1  /dV2  )  (  V  )  spectra  at  voltages  corresponding  to 
the  sum  of  the  two  scattering  energies.  Peak  ( a )  would  thus 
be  consistent  with  the  inelastic  excitation  of  two  A'- point 
transverse  acoustic  (TA)  phonons.  On  the  other  hand,  the 
broad  peak  (b)  may  result  from  a  number  of  scattering 
mechanisms  which  are  not  individually  resolved  in  the  spec¬ 
trum  shown  in  Fig.  3.  Such  inelastic  processes  could  be  the 
excitation  of  Af-point  longitudinal  acoustic  ( LA  ) ,  transverse 
optical  (TO),  and  longitudinal  optical  (LO)  phonons,  for 
example.  The  experimental  and  theoretical  J-V curves  pre¬ 
sented  in  Figs.  2  and  4  are  in  good  agreement. 

Another  possible  way  for  electrons  to  tunnel  through  the 
AlAs  band  gap  at  the  X  point  arises  from  the  coupling  of 
AlAs -T-point  states  to  GaAs  F -point  states  due  to  the  break¬ 
ing  of  translational  symmetry  in  the  direction  perpendicular 
to  the  layers.  This  makes  it  possible  for  T -point  electrons  in 
one  GaAs  electrode  to  tunnel  elastically  through  the  AlAs 
T-point  barrier  into  T-point  states  in  the  other  GaAs  elec¬ 
trode.  In  this  process,  the  virtual  T-point  states  to  consider 
are  those  below  the  two  X  valleys  lying  along  the  kx  direc¬ 
tion,  perpendicular  to  the  plane  of  the  layers.  Since  electrons 
travel  along  the  x  axis,  the  effective  mass  of  importance  in 
the  barrier  is  the  large  longitudinal  effective  mass  mf 
=  1.1m,.  Preliminary  calculations  indicate  that  the  cou¬ 
pling  of  states  having  different  symmetries  in  GaAs  and 
AlAs  may  be  significant  in  the  heterostuctures  studied  and 
could  yield  current  densities  large  enough  to  account  for  the 
discrepancies  observed  between  experimerttal  and  T-point 
elastic  tunneling  currents. 

We  have  proposed  two  separate  mechanisms  by  which  T- 
point  electrons  in  the  GaAs  electrodes  may  tunnel  through 
the  AlAs  Appoint  barrier.  In  fact,  both  processes  may  occur 
simultaneously,  and  the  empirical  factors  used  to  calculate 
theoretical  current  densities  can  be  modified  accordingly. 
More  detailed  calculations  of  scattering  matrix  elements  and 
the  mixing  of  states  of  different  symmetry  are  under  way  to 
obtain  a  quantitative  picture  of  electron  transport  in  these 
structures.  Nevertheless,  the  preliminary  results  presented 
here  suggest  that  (i)  elastic  tunneling  through  the  AlAs  T- 
point  barrier  is  insufficient  to  account  for  experimental  cur¬ 
rent  densities;  ( ii )  elastic  and  inelastic  tunneling  through  the 
AlAs  band  gap  at  the  X  point  could  account  for  the  observed 
discrepancies;  and  (iii)  theoretical  current  densities  in  good 
agreement  with  experimental  data  can  be  obtained  provided 
both  T-  and  T-point  tunneling  contributions  are  taken  into 
account. 

Although  these  concepts  have  only  been  illustrated  in 
terms  of  the  experimental  data  related  to  one  structure,  it 
was  found  that  samples  having  different  doping  concentra¬ 
tions  and  barrier  thicknesses  display  the  same  general  behav¬ 
ior.  Nevertheless,  further  investigations  are  necessary  to 
gain  a  better  understanding  of  the  specific  mechanisms  oc¬ 
curring  in  each  sample  and  to  build  a  unifying  model  for  all 
of  them  This  work  will  be  discussed  more  thoroughly  in  a 
forthcoming  publication. 

IV.  IMPURITY-ASSISTED  TUNNELING 

In  this  section,  a  new  element  is  added  to  the  current  den¬ 
sity  calculations,  namely  the  effect  of  impurity  encrgv  levels 


in  the  AlAs  band  gap.  Impurity  atoms  and  crystal  defects 
can  create  localized  energy  levels  (traps)  in  the  AlAs  bar¬ 
rier.  I"-2"  If  these  "impurity  states”  lie  in  the  AlAs  band  gap, 
they  can  affect  tunneling  currents  through  the  structure.  In 
this  paper,  we  view  traps  as  intermediate  states  through 
which  two-step  tunneling  processes  may  occur-  an  electron 
first  tunnels  from  one  GaAs  electrode  to  an  intermediate 
state,  and  then  from  the  intermediate  state  to  the  other  elec¬ 
trode.  This  approach  is  similar  to  that  used  by  Parker  and 
Mead  in  their  treatment  of  traps  in  Schottky  barriers.21 

A  simple  expression  for  the  two-step  tunneling  current 
can  be  derived  by  assuming  that  there  are  N  localized  trap 
states  per  unit  volume  of  AlAs,  with  each  trap  energy  level, 
E;.  being  a  fixed  energy,  E„,  below  the  AlAs  conduction 
band  edge.  Ec .  For  an  arbitrary  voltage  applied  to  the  struc¬ 
ture,  the  AlAs  conduction  band  edge  and  thus  the  trap  ener¬ 
gy  level  vary  with  position,  .t: 

E,  (x)  =  £,(*)  -  E0  .  ( 1 ) 

In  Fig.  5,  this  concept  is  illustrated  by  plotting  the  T-point 
conduction  band  edge  and  a  trap  energy  level  for  a  single p- 
type  barrier  heterostructure  under  applied  bias.  This  figure 
is  the  result  of  a  band  bending  calculation,  in  which  the  elec¬ 
trode  and  barrier  dopings  are  taken  to  be  3xl018  and 
1  X  I018  cm-’,  respectively.  The  barrier  thickness  is  100  A 
and  the  applied  voltage  is  150  mV.  The  impurity  level  is 


Distance  ! 

Fk.  5  Calculated  T-point  conduction  hand  edge  (solid  line)  as  a  function 
iff  distance  m  the  direction  perpendicular  to  the  layers  for  a  GaAs-AlAs 
CiuAs  single  harrier  structure  having  an  impurity  level  m  the  AlAs  harrier 
layer  The  impurity  state  (dashed  line)  lies  S'O  mcV  below  the  T  point 
conduction  band  edge  in  the  barrier  1  he  GaAs  electrodes  are  doped  n-tvpe 
at*  •  I0H  cm  1  The  AlAs  barrier  Liver  is  l(K)  thick  and  doped /7ts  pe  ji 
1  •  H)"1  cm  1  A  voltage  of  150  mV  is  applied  to  the  structure 
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taken  to  be  870  meV  below  the  AlAs  T -point  conduction 
band  edge 

The  elastic  tunneling  current  density.  J  ,  .  between  .1  trap 
state  at  position  a  and  the  left  GaAs  electrode  is  of  the  lot  m 

7,  oc  r,  [  E  (  a  .  )  ..v  ]  {  /.  [  E  <  a  )]-/•'[  E  t  v  i  ]  }  .  1  2 ) 

where  T ,  is  the  WK.B  transmission  probability  for  an  elec¬ 
tron  with  energy  E ,  (a„  )  tunneling  from  the  left  electrode  to 
x0  ,f  is  the  occupation  probability  of  the  left  electrode  state 
with  energy  E,  (x0 ),  and  Fis  the  occupation  probability  of 
the  trap  state.  Similarly,  the  elastic  tunneling  current  den¬ 
sity,  J2 ,  between  the  trap  state  at  position  xB  and  the  right 
electrode  is  given  by: 

y;<x  7\[£,  (a„  ),a„  ]{F  [F  (k„  )]—/,[£,  (a  , )  ] }  .  (3) 

where  7%  is  the  transmission  probability  for  an  electron  with 
energy  E,  (x0  )  tunneling  from  x0  to  the  right  electrode,  and 
fr,  is  the  occupation  probability  of  the  right  electrode  state 
with  energy  E,  (xB  ). 

Under  steady-state  conditions,  the  trap  level  occupation 
probability,  F,  will  adjust  itself  to  make  the  two  current  den¬ 
sities  equal.  It  follows  that  the  impurity-assisted  tunneling 
current,  y,mp,  takes  the  form: 

/  I  +  /  2 

A  comparison  of  Eq.  (4)  with  the  usual  expression  for  the 
one-step  elastic  tunneling  current,"  J ^ ,  rfiveals  that: 


J  vs.  Barrier  Thickness 
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where  /  (  E  )  is  the  I ransmission  coefficient  for  electrons 
with  energy  E  tunneling  from  the  left  electrode  to  the  right 
electrode,  and  C  is  a  constant  w  hich  is  proportional  to  the 
number  of  trap  levels  per  unit  volume.  .V.  and  the  effective 
cross-sectional  area  of  the  traps,  o 

Figure  b  is  a  calculated  plot  of  y„np  as  a  function  of  A I  As 
barrier  thickness.  The  calculations  are  performed  using  a  p- 
type  barrier  doping  of  1  -  10"'  cm-',  and  an  n-type  elec¬ 
trode  doping  of  3  x  10"*  cm  -  The  applied  voltage  is  100 
mV  and  the  impurity  energy  level  is  taken  to  be  950  meV 
below  the  AlAs  T -point  conduction  band  edge.  .V  is  chosen 
to  he  1  v  1017  cm  -  and  a  is  assigned  a  value  of  ten  times  the 
cross-sectional  area  of  the  AlAs  primitive  cell.  Changes  in  .V 
and  (j  alter  only  the  multiplicative  constant,  C.  in  Eq.  f  5  j, 
and  therefore,  do  not  change  the  qualitative  behavior  ofJ,m„ 
y„  is  also  calculated  and  plotted  in  Fig  6  for  comparison 
with  ymp .  This  reveals  that  the  contribution  to  the  total  cur¬ 
rent  density  of  one-step  elastic  tunneling  is  larger  than  that 
of  impurity-assisted  tunneling  in  heterostructures  with  thin 
barriers.  The  opposite  is  true  in  structures  with  thicker  bar¬ 
riers.  For  the  parameters  used  in  Fig.  6,  the  impurity-assist¬ 
ed  tunneling  current  becomes  larger  than  the  one-step  elastic- 
current  when  the  AlAs  burner  is  thicker  than  about  80  A. 


V.  IMPLICATIONS  FOR  RESONANT  TUNNELING 
STRUCTURES 

In  this  section,  we  briefly  discuss  implications  of  some  of 
the  previous  results  for  resonant  tunneling  in  GaAs/AlAs 
double  barrier  heterostructures.  Resonant  tunneling  struc¬ 
tures  made  of  two  Al,  Ga,  .  ,  As  quantum  barriers  separat¬ 
ed  by  a  GaAs  quantum  well  display  negative  differential  re¬ 
sistances  in  their  I-V  characteristics.  This  effect  is 
commonly  attributed  to  the  tunneling  of  electrons  via  quasi- 
stationary  states  in  the  well  Furthermore,  it  is  usually  as¬ 
sumed  that  the  effective  barrier  height  for  tunneling  is  deter¬ 
mined  by  the  T -point  conduction  band  in  the  barriers,  even 
when  the  Al,  Ga,  ,  As  is  indirect.  Detailed  calculations  of 
band  bending,  energy  level  positions,  and  voltage  drops  are 
required  to  determine  without  ambiguity,  the  correlation 
between  the  resonances  in  the  well  and  the  voltages  at  which 
negative  differential  resistances  are  observed.  Simple  theo¬ 
retical  models  yield  current  densities  in  disagreement  with 
experimental  data.  In  structures  for  which  the  best  peak-to- 
vallev  current  ratios  realized  to  date  are  10:1  at  77  K.::  cal¬ 
culations  predict  values  as  large  as  several  orders  of  magni¬ 
tude  1  This  indicates  that  the  transport  properties  of  these 
stuctures  are  more  complex  than  simple  models  assess.  In 
particular,  the  obtainable  values  of  peak-to-v alley  current 
ratios  are  limited  by  other  current  transport  mechanisms 
competing  with  ("  point  resonant  tunneling  "'  Derivative 
spectra  are  important  tools  for  identifying  and  investigating 
the  nature  of  such  processes  Figure  1  illustrates  second  deri¬ 
vative  curves  u/  '  /  <ii  )  i  I  ).  obtained  at  4  2  K  for  a  reso¬ 
nant  tunneling  heterostructure  The  GaAs  electrodes  in  the 
sample  are  doped  n-tvpe  with  Se  at  2  ■  1 0 1 "  cm  '.  The  AlAs 
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Fig.  7  Second  derivative  spectra  (d;  I  /dV  '~  )  (  V  ),  at  4  2  K  for  a  GaAs/ 
AlAs  resonant  tunneling  double  barrier  heterostructure.  The  GaAs  elec¬ 
trodes  are  doped  retype  at  2x  10'“  cm  The  AlAs  quantum  barriers, 
approximately  50  A  thick,  are  doped  p-type.  The  nominally  undoped  GaAs 
well  is  about  50  A  wide  The  curves  (a),  (b),  and  (c)  correspond  to  three 
different  devices,  having  diameters  of  145,  1 15,  and  75  pm. 

barriers,  about  50  A  thick,  are  doped  p-type  with  Mg.  The 
nominally  undoped  GaAs  well  is  about  50  A  wide.  Three 
curves  are  shown  in  Fig.  7,  corresponding  to  different  size 
devices.  They  display  reproducible  features  which  have  a 
nearly  periodic  spacing  of  15-16  mV.  Similar  spectra  are 
obtained  from  other  samples.  A  definitive  explanation  of  the 
observed  features  has  not  yet  been  found,  but  a  possible  in¬ 
terpretation  may  be  proposed.  The  periodic  structure  could 
arise  from  the  discrete  nature  of  fluctuations  in  the  quantum 
well  width.  Such  fluctuations  in  barrier  and  well  thicknesses 
would  contribute  in  reducing  and  broadening  the  negative 
differential  resistances  observed  in  the  I-V characteristics. 

In  the  light  of  the  results  obtained  for  single  barrier  struc¬ 
tures,  it  is  possible  to  propose  other  qualitative  explanations 
for  the  small  peak-to-valley  current  ratios  observed  in  sever¬ 
al  of  our  samples.2  Whenever  electrons  can  scatter  to  X- 
point  states  and  tunnel  through  the  AlAs  barriers  at  the  X 
point,  a  large  background  current  can  be  expected,  even  for 
relatively  thick  barriers  and  wells.  Depending  upon  the  rela¬ 
tive  amplitudes  of  the  resonant  P-point  current  and  the  X- 
point  current,  the  negative  differential  resistances  may  be 
partially  or  totally  suppressed.  Determining  whether  or  not 
this  will  happen  is,  in  fact,  a  nontrivial  problem  which  de¬ 
pends  upon  the  structure  parameters  and  scattering  mecha¬ 
nisms.  In  single  barrier  samples,  electrons  may  tunnel 
through  the  AlAs  band  gap  at  the .V  point  either  inelasticallv 
by  scattering  near  each  interface,  or  elastically  through  the 
coupling  of  AlAs  A'-point  states  and  GaAs  T -point  states.  In 
double  barrier  structures,  the  problem  is  different  due  to  the 
presenceof  the  GaAs  well  between  the  two  AlAs  barriers.  At 
voltages  which  are  low  enough  that  the  A'-point  conduction 
band  edge  in  the  GaAs  well  is  higher  in  energy  than  the 
Fermi  level  in  the  left  electrode.  A-point  tunneling  may  oc¬ 
cur  through  several  mechanisms  First,  electrons  may  scat¬ 
ter  once  near  each  of  the  four  GaAs/AIAs  interfaces,  thus 
tunneling  through  the  AlAs  barriers  at  the  X  point,  but  occu- 
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pying  states  of  T  symmetry  in  the  GaAs  electrodes  and  well. 
Secondly,  electrons  may  scatter  only  near  the  first  and  last 
GuAs/AlAs  interfaces,  thus  tunneling  in  both  barriers  and 
the  well  through  the  A'-point  gaps  of  AlAs  and  GaAs.  Final¬ 
ly,  electrons  may  tunnel  elastically  due  to  the  coupling  of 
wave  functions  having  different  symmetries.  Simple  esti¬ 
mates  hav£  shown  that  the  relative  importance  of  these 
mechanisms  depends  on  the  well  width,  barrier  thicknesses, 
and  applied  voltage.  These  processes  either  involve  higher 
order  mechanisms  than  those  taking  place  in  single  barrier 
structures,  or  tunneling  through  thick  A'-point  barriers  con¬ 
sisting  of  the  two  AlAs  barriers  and  the  GaAs  well.  There¬ 
fore,  their  contribution  can  certainly  be  minimized  by  care¬ 
fully  optimizing  the  structure  parameters,  improving  the 
quality  of  the  interfaces,  and  lowering  the  background  dop¬ 
ing  in  nominally  undoped  layers. 


The  objective  of  this  study  was  to  investigate  elastic  and 
inelastic  tunneling  processes  in  GaAs-AlAs-GaAs  double 
heterojunctions  grown  by  MOCVD  in  the  [  100]  direction. 
The  AlAs  quantum  barriers  in  the  heterostructures  studied 
were  doped  p-type  with  Mg.  Measurements  of  current-vol¬ 
tage  characteristics  (I-V  ),as  well  as  first  (dl  /dV )  and  sec¬ 
ond  (d2 I /dV2  )  derivatives  of  the  I-V  curves,  were  per¬ 
formed  at  temperatures  ranging  from  300  to  4.2  K. 
Experimental  and  theoretical  results  indicated  that  low  tem¬ 
perature  electron  transport  in  structures  with  thin  barriers  is 
achieved  through  two  mechanisms:  (i)elastic  electron  tun¬ 
neling  through  the  AlAs  T-point  barrier;  and  (ii)  inelastic 
and/or  elastic  tunneling  through  the  AlAs  A’-point  barrier. 
Inelastic  A'-point  tunneling  may  occur  via  virtual  states  be¬ 
neath  the  four  A'-point  valleys  lying  parallel  to  the  plane  of 
the  layers,  and  requires  the  scattering  of  electrons  near  each 
interface.  Elastic  A'-point  tunneling,  due  to  the  coupling  of 
AlAs  A'-point  states  and  GaAs  T -point  states,  can  take  place 
via  virtual  states  below  the  two  AlAs  A'-point  valleys  lying  in 
the  direction  perpendicular  to  the  plane  of  the  layers.  A  theo¬ 
retical  model,  which  treats  trap  levels  in  the  AlAs  barrier  as 
intermediate  states  for  two-step  tunneling  processes,  was  de¬ 
veloped.  Calculations  indicated  that  this  impurity-assisted 
tunneling  mechanism  can  become  important  when  the  AlAs 
barrier  is  thick  enough.  Implications  of  some  of  these  results 
for  resonant  tunneling  heterostructures  consisting  of  two 
AlAs  quantum  barriers  separated  by  a  GaAs  quantum  well 
were  discussed.  In  particular,  carefully  choosing  the  param¬ 
eters  of  these  structures  to  maximize  the  contribution  of  the 
T -point  resonant  tunneling  current  is  critical  for  increasing 
the  peak-to-valley  current  ratios  of  the  negative  differential 
resistances.  For  these  double  barrier  structures,  we  also  pre¬ 
sented  experimental  second  derivative  spectra  (d2 1  /dV:  ). 
in  which  periodic  features  may  be  indicative  of  discrete  fluc¬ 
tuations  in  well  thickness. 
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Two  novel  three-terminal  devices  based  on  tunneling  in  quantum  well  and  quantum  barrier 
heterostructures  are  proposed  and  analyzed  theoretically.  In  both  devices,  the  relative  positions  of 
the  base  and  collector  are  interchanged  from  conventional  emitter-base-collector  sequence.  This 
provides  a  means  for  obtaining  negligible  base  currents  and  large  current  transfer  ratios.  In  both 
cases,  a  base  voltage  controls  the  emitter-collector  tunneling  current  by  shifting  the  resonances  in 
a  quantum  well.  Calculations  indicate  that  significant  variations  in  the  emitter-collector  current- 
voltage  characteristics  can  be  obtained  for  reasonable  base-emitter  voltages.  We  call  the  two 
devices  a  Stark  effect  transistor  and  a  negative  resistance  Stark  effect  transistor,  respectively. 


Electronic  tunnel  structures  are  a  source  of  increasing 
interest.  Two-terminal  devices  based  on  tunneling  in  single 
and  double  barrier  heterostructures  have  been  studied. ,_1> 
While  it  is  true  that  these  two-terminal  structures  have  a 
number  of  potential  applications,  three-terminal  devices 
would  be  preferable  in  many  cases. 

The  idea  of  making  three-terminal  devices  based  on  tun¬ 
neling  was  pioneered  in  1960  by  Mead. 10  He  proposed  a  met- 
al-insulator-metal-insulator-metal  hot-electron  transistor. 
A  number  of  devices  based  on  this  concept  have  been  pro¬ 
posed  and  investigated  experimentally. 11-13  All  of  these 
structures  suffer  from  very  small  emitter-collector  current 
transfer  ratios  due  to  large  base  currents. 

Recently,  Jogai  and  Wang'4  calculated  the  tunneling 
current  for  a  conceptual  three-terminal,  double  barrier  de¬ 
vice  consisting  of  alternating  layers  of  GaAs  and 
AIj.Ga,  _  x  As,  forming  the  emitter,  base,  and  collector. 
They  proposed  a  configuration  in  which  the  base  contact 
would  be  made  to  the  GaAs  quantum  well  and  the  two  bar¬ 
riers  independently  biased.  A  highly  conductive  base  was 
required  so  that  a  potential  could  be  applied  to  it,  but  no  base 
current  was  allowed  to  flow.  This  is  a  somewhat  unphysical 
and  unrealistic  assumption.  In  addition,  varying  the  collec¬ 
tor-base  voltage  did  not  produce  very  significant  changes  in 
the  device  current- voltage  characteristics 

In  this  letter,  we  propose  and  analyze  two  three-termi¬ 
nal  devices  which  we  expect  to  have  reduced  base  currents 
and  improved  device  performances.  Although  other  semi¬ 
conductors  could  be  used,  the  devices  are  presented  here  in 
the  context  of  GaAs/Al,  Ga,  ... ,  As  heterojunction  technol¬ 
ogy.  Since  tunneling  is  the  main  current  transport  mecha¬ 
nism,  these  devices  should  feature  the  high-speed  character¬ 
istics  associated  with  tunnel  structures.  They  should  also 
have  large  emitter-collector  current  transfer  ratios.  The  key 
step  in  achieving  this  goal  is  to  interchange  the  relative  posi¬ 
tions  of  the  base  and  collector,  thus  locating  the  latter  in  the 
region  where  current  is  most  likely  to  flow.  Figure  1  shows 
schematically  the  first  proposed  structure,  together  with  its 
energy-band  diagram  at  equilibrium.  The  emitter  is  an  n- 
type  GaAs  layer.  Doping  concentrations  on  the  order  of 
5  x  10'6-5  X  10 1 7  cm  '  1  should  provide  sufficient  tunneling 
currents  and,  at  the  same  time,  allow  resonant  tunneling  ef¬ 
fects  to  be  observed  at  room  temperature.4  The  emitter  is 
followed  by  a  thin  undoped  Al,  Ga,  ,  As  tunneling  barrier 


and  a  lightly  doped  n-type  quantum  well.  The  collector  con¬ 
tact  is  made  to  the  well.  Alloys  which  give  shallow  and 
abrupt  ohmic  contacts,  such  as  Au/Ge-Ag-Au  or  Au/Ge- 
Au-WSi  for  example,  should  be  used.  The  next  layer  is  light¬ 
ly  doped  A1  Ga,  _  As  barrier,  sufficiently  thick  to  prevent 
electrons  from  tunneling  through  it.  In  addition,  an  alloy 
with  large  A1  composition  is  desirable  to  minimize  ther¬ 
mionic  emission  over  the  barrier.  Finally,  a  heavily  doped  rt- 
type  GaAs  substrate  serves  as  the  base  electrode. 

The  principle  of  operation  of  the  proposed  device  is  de¬ 
scribed  below.  First,  let  us  consider  the  case  in  which  no  base 
voltage  is  applied.  When  the  collector  is  biased  positively 
with  respect  to  the  emitter,  electrons  near  the  Fermi  level  in 
the  emitter  tunnel  through  the  thin  AlXiGa,_,  As  barrier 
into  the  collector.  As  long  as  the  emitter  Fermi  level  remains 
below  the  first  subband  in  the  well,  a  negligible  tunneling 
current  is  expected.  When  the  bias  voltage  is  such  that  the 
emitter  Fermi  level  reaches  the  first  subband,  the  current  is 
significantly  increased.  The  current-voltage  (/-F)  character¬ 
istics  should  thus  feature  enhancements  corresponding  to 
the  alignment  of  the  emitter  Fermi  level  with  each  resonance 
in  the  well.  If  a  potential  difference  is  now  applied  between 


(hi 


FIG  I  Schematic  diagrams  (not  lo  scale!  of  la)  a  cross  section  of  the  pro¬ 
posed  Stark  effect  transistor  (SET),  (b)  the  conduction-band  edge  at  equilib¬ 
rium  as  a  function  of  position  in  the  x  direction  (perpendicular  to  the  laversl 
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FIG.  2.  Calculated  current-voltage  charactensttcs  for  the  Stark  effect  tran¬ 
sistor  (SET|.  The  barriers  are  pure  AlAs,  50  and  1000  A  thick,  respectively 
The  well  is  a  50-A-thick  GaAs  layer  The  conduction-band  offset  is  taken  to 
be  0.96  eV. 


the  base  and  emitter,  an  electric  field  perpendicular  to  the 
layers  is  created.  The  field  will  modify  the  positions  of  the 
subbands  in  the  well  with  respect  to  the  emitter  Fermi  level 
and  thus  modulate  the  tunneling  current.  The  field  will  pene¬ 
trate  into  the  quantum  well  region  for  the  following  reasons: 
(i)  the  barriers  and  well  are  lightly  doped;  (ii)  the 
Alx_  Ga,  _  Xi  As  barrier  is  not  thick  enough  to  drop  all  of  the 
base-emitter  voltage;  (iii)  the  device  geometry  is  such  that  the 
collector  contact  does  not  completely  shield  the  emitter  from 
the  base. 

Theoretical  /-  V  characteristics  for  the  device  are  shown 
in  Fig.  2.  These  curves  were  calculated  using  Bardeen’s 
many-particle  tunneling  formalism.15  In  this  approach, 
which  starts  from  Fermi’s  Golden  Rule,  the  density  of  finals 
states,  pf,  appears  explicitly  in  the  calculation.  In  the  present 
case,  pf  is  simply  a  sum  of  step  functions  since  each  reso¬ 
nance  in  the  well  is  the  bottom  of  a  two-dimensional  energy 
band.  The  matrix  element  was  calculated  by  using  WKB 
wave  functions  in  the  barrier  region.  A  two-band  model ,  k-p 
theory  calculation  was  used  to  obtain  the  complex  band 
structure  in  the  barrier.  Because  the  barriers  and  well  are 
lightly  doped,  the  base-emitter  voltage  V BE  was  assumed  to 
drop  linearly  in  those  regions.  The  calculated  curves  display 
the  substantial  variations  of  the  emitter-collector  I-V  char¬ 
acteristics  which  may  be  obtained  by  modulating  V  8E .  Base- 
emitter  voltages  much  higher  than  those  needed  to  produce 
significant  transistor  action  can  be  applied  without  produc¬ 
ing  avalanche  breakdowns.  If  V 8E  is  negative,  the  levels  in 
the  well  are  shifted  upwards  and  a  small  negative  current 
might  flow  to  the  emitter  when  FCE  is  small.  This  is  an  addi¬ 
tional  incentive  to  make  the  well  and  the  Alx  Ga,  .  X)  As  bar¬ 
rier  lightly  doped. 

The  main  advantage  of  this  configuration  is  a  negligible 
base  current  and  thus  a  large  current  transfer  ratio.  Because 
the  quasistationary  states  in  the  well  are  modulated  by  an 
electric  field  to  produce  transistor  I-V  characteristics,  the 
proposed  device  could  be  called  a  Stark  effect  transistor 
(SET). 

The  transistor  1-V  characteristics  described  above  can 
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FIG.  3.  Schematic  diagrams  (not  to  scale)  of  |a|  a  cross  section  of  the  pro¬ 
posed  negative  resistance  Stark  effect  transistor  (NERSET);  (b)  the  conduc¬ 
tion-band  edge  at  equilibrium  as  a  function  of  position  in  the  x  direction 
(perpendicular  to  the  layers). 


be  modified  and  enhanced  by  adding  a  potential  step  in  the 
GaAs  well  constituting  the  collector.  The  step  can  be  a  thin 
Al.Ga,  _xAs  layer  with  an  A1  composition  smaller  than 
that  of  the  barriers.  Its  main  effect  is  to  modify  the  relative 
positions  of  the  resonant  states  in  the  well.  The  structure  can 
now  be  designed  to  obtain  a  certain  energy  spectrum  or  to 
produce  given  shifts  of  some  of  the  subbands.  This  makes  it 
possible  to  optimize  device  performance  according  to  the 
requirements  or  applications. 

To  further  illustrate  the  concept  of  subband  modulation 
by  means  of  an  electric  field  applied  from  a  controlling  elec¬ 
trode,  another  device  configuration  is  proposed.  It  is  charac¬ 
terized  by  a  resonant  tunneling  double  barrier  heterostruc¬ 
ture.  Such  a  device  structure  is  shown  schematically  in  Fig. 
3,  together  with  its  energy-band  diagram  at  equilibrium. 
This  device  operates  on  the  same  principle  as  the  Stark  effect 
transistor.  The  values  of  the  fields  and  the  shifts  of  the  levels 
resulting  from  the  base-emitter  modulating  voltage  were  es¬ 
timated  to  be  of  the  same  order  of  magnitude  as  before.  In 
this  case,  however,  no  contact  needs  to  be  made  to  the  quan¬ 
tum  well.  The  essential  feature  of  this  device  configuration  is 
the  presence,  in  the  emitter-collector  I-V  characteristics,  of 
negative  differential  resistances  controlled  by  the  base-emit¬ 
ter  voltage. 

Theoretical  1-V  characteristics  for  the  device  are  shown 
in  Fig.  4.  The  approach  of  Vassell  et  a/.16  was  used  to  calcu¬ 
late  the  tunneling  current  through  the  double  barrier  hetero- 
structure  separating  the  collector  from  the  emitter.  The 
base-emitter  modulating  voltage  was  once  again  assumed  to 
drop  linearly  in  each  region.  This  device,  which  could  be 
called  a  negative  resistance  Stark  effect  transistor  (NER¬ 
SET),  offers  several  advantages  over  any  double  barrier  con¬ 
figuration  in  which  the  base  is  located  in  the  quantum  well 
(l)  easier  fabrication;  (ii)  no  base  current,  and  thus  larger  cur- 
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FIG  4.  Calculated  current-voltage  charactenstics  for  the  negative  resis¬ 
tance  Stark  effect  transistor  NERSET)  The  tunnel  bamers  are  20-A-thick 
Alo.Gao,  As  layers  and  the  well  a  50-A-thick  GaAs  layer  The  conduction- 
band  offset  is  taken  to  be  0. 5  e V . 

rent  transfer  ratios;  (iii)  undoped  well;  jiv)  no  transverse  elec¬ 
tric  fields  in  the  well  which  destroy  the  coherence  of  wave 
functions  across  the  entire  double  barrier  structure.  These 
last  two  conditions  are  essential  for  optimum  resonant  tun¬ 
neling. 

In  summary,  we  have  proposed  and  analyzed  two  novel 
three-terminal  devices  based  on  tunneling  in  quantum  well 
and  quantum  barrier  heterostructures.  Their  main  charac¬ 
teristic  is  that  the  relative  positions  of  the  collector  and  base 
electrodes  have  been  interchanged  with  respect  to  the  con¬ 
ventional  emitter-base-collector  sequence.  This  makes  it 
possible  to  obtain  negligible  base  currents  and  large  current 
transfer  ratios.  Electric  fields  produced  by  applying  a  vol¬ 
tage  to  the  base  modulate  the  positions  of  the  subband  levels 
in  the  quantum  well  and  thus  control  the  emitter-collector 
tunneling  current.  Calculations  showed  that  significant  var¬ 
iations  in  the  emitter-collector  /-  V charactenstics  can  be  ob¬ 
tained  by  modulating  the  base-emitter  voltage.  Although 


these  devices  were  analyzed  in  the  context  of  GaAs/ 
Al.Ga,  , As  heterojunction  technology,  other  semicon¬ 
ductors  could  be  used.  Experiments  are  under  way  to  mea¬ 
sure  the  properties  of  these  structures  and  to  explore  their 
applications  in  high-frequency  analog  and  digital  circuits. 

John  Lambe  played  a  key  role  in  calling  our  attention  to 
the  possible  use  of  electnc  fields  in  changing  the  positions  of 
subbands.  A.  Zur  developed  the  computer  programs  used  to 
calculate  the  I-V  charactenstics  of  the  NERSET.  M.  B 
Johnson,  T.  E.  Schlesinger,  T.  K.  Woodward,  and  G.  Y.  Wu 
have  all  helped  us  in  the  clarification  of  our  thinking  on  these 
device  concepts.  This  work  was  supported  in  part  by  the 
Defense  Advanced  Research  Projects  Agency  under  con¬ 
tract  No.  N00014-84-C-0083  and  the  Office  of  Naval  Re¬ 
search  under  contract  No.  N00014-82-K-0556. 


'D  Delagebeaudeuf,  P  Delescluse,  P  Etienne,  J.  Missies.  M.  Laviron.  J 
Chaplart,  and  N.  T.  Linh,  Electron.  Lett.  18,  85  1 1982). 

-R  T  Collins,  J  Lambe.  T.  C  McGill,  and  R  D  Burnham.  Appl  Phys. 
Lett.  44,  532(1984). 

’I  Hase.  H.  ICawai,  K.  ICineko,  and  N.  Watanabc,  Electron.  Lett.  20,  491 

11984). 

*R  Tsu  and  L.  Esaki,  Appl.  Phys.  Lett.  22.  562  (1973). 

'T  C  L.  G.  Sollner,  W  D.  Goodhue,  P  E.  Tannenwald.  C.  D  Parker,  and 
D  D  Peck,  Appl.  Phys.  Lett.  43.  588  11983). 

'•R  T  Collins,  A.  R.  Bonnefoi,  J  Lambe,  T  C  McGill,  and  R  D  Burn¬ 
ham,  Proceedings  of  the  /  7th  International  Conference  on  the  Physics  of 
Semiconductors.  San  Francisco,  1984  (Sponger,  New  York.  1984i,  p  43T 
’T  C  L.  G  Sollner.  P  E.  Tannenwald.  D  D  Peck,  and  W  D  Goodhue. 
Appl  Phys.  Lett  45.  1319(1984! 

"A  R  Bonnefoi.  R  T  Collins,  T  C  McGill,  R  D  Burnham,  and  F  A 
Ponce,  Appl  Phys  Lett.  46,  285  (1985). 

T  J  Shewchuk.  P  C  Chapin.  P  D  Coleman,  W  Kopp,  R  Fischer,  and 
H  Morkos.  Appl.  Phys.  Lett.  46.  508  (19851 
|0C.  A.  Mead,  Proc.  IRE  48.  359  (I960) 

"D  v  Geppert.  Proc.  IRE  50,  1527(1962) 
l!M  Heiblum.  Solid  State  Electron.  24.  343  (1981) 

"N.  Yokoyama,  K.  Imamura,  T  Ohshima,  H.  Nishi,  S  Muto,  K  Rondo, 
and  S  Hiyamizu,  Ipn.  J  Appl  Phys.  23.  L31 1  1 1 984) 

,4B  Jogai  and  K.  L.  Wang,  Appl  Phys.  Lett.  46.  167  119851 
’'1  Bardeen,  Phys.  Rev  Lett  6.  57  (19611 

"’M.O  V asset!,  Johnson  Lee,  and  H.  F  Lockwood.)  Appl  Phys  54.520 

(19831. 


890 


Appl  Phys  Lett  Voi  47  No  8  '  5  October  '988 


Bonnefoi  Chow,  ana  McGill 


890 


Paper  7 


Photovoltaic  investigations  of  GaAs/AIAs  heterostructures 

T.  E.  Schiesinger,  R.  T.  Collins,  and  T.  C.  McGill 

California  Institute  of  Technology.  Pasadena.  California  W 125 

R  D.  Burnham 

Xerox  Corporation.  Palo  Alio.  California  1 4}()4 

(Received  2'  May  1984;  accepted  for  publication  3  July  1984) 

We  present  the  results  of  the  application  of  photoresponse  techniques  to  the  study  of  the  transport 
of  electrons  past  an  energy  barrier.  In  this  study,  the  barrier  was  provided  by  a  thin  layer  of  AlAs 
sandwiched  between  GaAs  layers.  The  experiment  measured  the  voltage  resulting  from  the 
migration  of  optically  excited  electrons  from  one  side  of  the  barrier  to  the  other.  The  voltage  is 
measured  as  a  function  of  the  wavelength  of  the  incident  light.  We  also  present  the  results  of 
calculations  which  explain  the  nature  of  the  observed  spectra  and  how  they  change  when  the 
thickness  of  the  top  layer  (the  illuminated  side)  of  GaAs  is  changed. 


In  this  letter  we  discuss  experimental  results  of  the  ap¬ 
plication  of  photoresponse  techniques  to  the  study  of  GaAs/ 
AlAs  heterostructures.  We  have  applied  this  technique  to 
the  study  of  the  transport  of  electrons  past  an  energy  barrier. 
In  particular,  we  have  investigated  these  electronic  transport 
properties  in  structures  consisting  of  a  thin  layer  of  AlAs 
sandwiched  between  layers  of  GaAs.  The  experimental  pro¬ 
cedure  consisted  of  illuminating  samples  with  light  in  order 
to  excite  electrons  in  the  conduction  band  of  the  GaAs  which 
could  tfien  migrate  through  the  AlAs  barrier.  We  then  mea¬ 
sured  the  resulting  open  circuit  voltage  produced  across  the 
AlAs  barrier  as  a  function  of  the  wavelength  of  the  incident 
light. 

The  basic  results  presented  in  this  letter  can  be  summar¬ 
ized  as  follows.  The  sign  of  the  voltage  that  we  measured  was 
consistent  with  electrons  migrating  from  the  illuminated 
side  of  the  AlAs  barrier  to  the  back  side  of  the  AlAs  barrier. 
The  magnitude  of  the  voltage  was  proportional  to  the  inten¬ 
sity  of  the  illuminating  light.  The  voltage  spectrum,  that  is, 
the  voltage  as  a  function  of  the  wavelength  of  the  incident 
light,  had  two  basic  shapes  which  were  common  to  all  the 
AlAs  barrier  samples  so  far  studied.  The  two  shapes  seemed 
mainly  to  depend  on  the  thickness  of  the  top  layer  of  GaAs. 
When  the  top  layer  was  on  the  order  of  a  few  microns  thick, 
the  resulting  spectrum  was  a  broad  peak  whose  exact  posi¬ 
tion  depended  on  doping  and  temperature  When  the  top 
layer  was  thinner  (less  than  one  micron  thick)  the  voltage 
spectrum  rose  slowly  from  shorter  wavelengths  to  longer 
wavelengths  followed  by  a  sharp  turnoff  at  some  particular 
wavelength,  depending  on  the  sample  doping  and  the  tem¬ 
perature. 

The  samples  used  in  this  study  were  grown  by  a  metal- 
orgamc  chemical  vapor  depostion  technique  '  The  AlAs  lay¬ 
er  thicknesses  varied  from  about  150  to  about  250  A,  and 
these  layers  were  dopeJ  p  type  with  Mg  at  about  1  X  10'" 
cm  ’  The  GaAs  top  layer  (the  illuminated  side)  varied  from 
about  0  5  to  4  4  thick,  depending  on  the  sample,  while 
the  bottom  layer  was  substantially  thicker  The  GaAs  was 
degenerately  doped  n  type  with  Se  from  1  to  5  x  10'*  cm  \ 
For  these  doping  levels  the  AlAs  barrier  is  fully  depleted 
Ohmic  contacts  were  made  to  the  samples  by  evaporation  of 
Au/Ge  followed  by  a  20  s  anneal  at  420  *C  Photolithogra¬ 
phic  techniques  were  employed  to  make  ring-shaped  con¬ 


tacts  to  the  top  surface  in  order  to  allow  light  to  penetrate 
into  the  GaAs  through  the  center  of  the  rings.  The  individual 
contacts  were  then  mesa  isolated  from  each  other.  The  sam¬ 
ples  were  mounted  on  transistor  headers,  and  wire  bonds 
we/e  made  to  the  Au/Ge  contacts.  The  light  source  used  to 
illuminate  the  samples  was  a  tungsten  filament  lamp.  The 
light  was  directed  through  a  spectrometer  (SPEX  model 
1 269)  and  then  focused  onto  a  particular  device  to  be  studied. 
The  measurement  of  the  voltage  resulting  from  the  illumina¬ 
tion  was  made  directly,  using  a  voltmeter  when  it  was  de¬ 
sired  to  know  the  sign  of  the  voltage.  In  general,  however, 
the  light  was  chopped  at  about  270  Hz  and  the  resulting 
voltage  was  measured  using  a  lock-in  amplifier.  The  signal 
was  then  digitized,  and  the  spectrum  stored  on  a  signal  aver¬ 
ager  (NIC  1 1 10).  In  this  letter,  we  discuss  the  results  ob¬ 
tained  for  only  one  particular  sample  at  room  temperature 
(300  K).  Further  results  describing  the  variations  in  the  spec¬ 
tra  from  sample  to  sample  as  a  function  of  doping  and  tem¬ 
perature  will  be  discussed  in  a  forthcoming  publication. 

To  model  the  absorption  of  the  light  in  these  structures, 
we  earned  out  the  following  calculations.  We  consider  a 
model  where  we  have  a  layer  of  GaAs,  a  layer  of  AlAs,  and  a 
semi-infinite  layer  of  GaAs.  Light  impinges  on  the  first  layer 
of  GaAs  from  a  vacuum  at  normal  incidence.  We  assume 
that  no  attenuation  of  the  light  takes  place  in  the  AlAs.  In 
the  GaAs,  however,  we  model  the  attenuation  of  the  light, 
using  a  coefficient  of  absorption,  a,  which  we  obtained  by 
fitting  and  extrapolating  the  curves  publised  by  Casey  et  al } 
for  heavily  doped,  n-type  GaAs  at  room  temperature.  We 
used,  however,  constant  values  for  the  index  of  refraction  of 
GaAs  and  AlAs  of  3.4  and  3.1,  respectively.3  The  light  inten¬ 
sity  was  then  calculated  throughout  the  model  structure  and 
was  then  integrated  in  the  region  of  GaAs  in  front  of  the 
AlAs  barrier  and  in  an  equal  region  behind  the  AlAs  barrier 
to  obtain  the  total  light  intensity  in  each  of  the  two  regions. 
The  total  light  intensity  behind  the  AlAs  barrier  was  then 
subtracted  from  the  total  in  front  of  the  AlAs,  and  this  differ¬ 
ence  was  plotted  as  a  function  of  wavelength.  When  there  is 
more  light  intensity  in  front  of  the  AlAs  than  behind,  there 
are  correspondingly  more  electrons  excited  there.  It  is  this 
difference  in  the  concentration  of  electrons  in  the  conduc¬ 
tion  band  which  drives  the  flow  of  electrons  from  the  front 
side  of  the  AlAs  barrier  to  the  back  side.  As  will  be  discussed 
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FIG.  3  Calculated  difference  in  the  integrated  light  intensities  in  front  of 
and  behind  the  AIAs  barrier  as  a  function  of  wavelength  The  coefficient  of 
absorption,  a,  used  in  calculating  these  curves,  was  for  a  GaAs  rt-type  dop¬ 
ing  of  3.3  X  10'*  cm  1  at  room  temperature.  The  results  are  presented  for 
top  layer  thicknesses  of  GaAs  of  4.4  /im  and  0.9  /im  as  indicated  in  the 
figure. 

intensity  of  the  light  is  much  greater  at  the  shorter  wave¬ 
length  tail  of  the  spectrum.  Since  the  number  of  electrons 
promoted  by  process  A  is  proportional  to  the  light  intensity, 
we  would  expect  the  voltage  signal  to  be  most  enhanced  at 
the  short  wavelength  tail  of  the  spectrum.  Thus,  a  spectrum 
for  a  sample  with  a  thick  top  layer  will  change  qualitatively 
in  the  way  it  does  in  Fig.  1  when  the  top  layer  is  made  thin¬ 
ner.  The  r -point  barrier  at  room  temperature  is  at  most5 
1.35  eV'  (9180  A).  This  barrier  is  decreased  by  the  fact  that 
the  Fermi  level  is  in  the  conduction  band  but  is  increased  by 
the  effects  of  band-gap  shrinkage  in  the  GaAs  due  to  the 
heavy  doping''  and  is  also  affected  by  band  bending  in  these 
structures.  This  number,  however,  gives  us  an  estimate  as  to 
the  longest  wavelength  at  which  we  could  expect  to  measure 
a  voltage  if  there  is  also  a  difference  in  the  light  intensity  in 
front  of  and  behind  the  AIAs  barrier.  At  wavelengths  longer 
than  9 1 80  A  even  electrons  promoted  via  process  A  may  not 
be  high  enough  in  the  conduction  band  to  go  over  the  barrier. 

In  Fig.  3  we  present  results  of  the  predicted  spectrum 
for  a  sample  with  a  doping  level  of  3.3 X  10IK  cm"’  in  the 
GaAs,  which  has  a  240- A  layer  of  AIAs  and  which  has  a  top 
layer  thickness  of  4.4  /um  or  0.9  /im  as  labeled.  We  see  that 
since  the  absorption  coefficient  data  contain  in  it  the  effects 
of  band  to  band  versus  free-carrier  absorption,  they  correctly 
predict  the  qualitative  behavior  of  the  spectra  we  measured. 
We  note  that  the  peak  position  |in  the  case  of  the  thick  top 
layer)  and  the  turnoff  wavelength  (in  the  case  of  the  thin  top 
layer)  are  not  exactly  predicted  by  the  calculations  though 
they  are  in  fairly  good  agreement.  These  positions  depend 
very  much  on  the  doping  of  the  sample  used  in  the  experi¬ 
ment.  In  the  calculations  these  positions  depend  on  the  be¬ 


havior  of  the  coefficient  of  absorption,  a,  which  in  turn  de¬ 
pends  on  the  doping  and  which  differs  from  the  actual 
doping  of  the  samples.  In  our  calculations,  we  have  assumed 
a  constant  doping  profile  while  in  the  actual  samples  the 
doping  in  the  GaAs  is  a  function  of  position.  The  difference 
in  peak  positions  for  samples  of  different  dopings  is  correctly 
predicted  by  the  calculations  and  will  be  presented  in  a  forth¬ 
coming  publication.  We  also  note  that  the  undulations  seen 
in  the  experimental  peak  are  not  as  marked  in  the  calculated 
peak.  There  are  two  possible  origins  for  these  undulations. 
One  explanation  is  that  they  are  the  result  of  the  optical 
properties  of  the  structures,  in  particular,  the  detailed  behav¬ 
ior  ofa.  Another  possible  explanation  is  that  they  area  result 
of  resonant  effects  in  the  transmission  of  electrons  through 
the  AIAs  barrier.4 

In  summary,  we  have  presented  results  of  an  experi¬ 
mental  technique  that  we  have  employed  to  study  the  trans¬ 
port  of  electrons  past  an  energy  barrier.  We  have  presented 
experimental  results  which  imply  that  the  voltage  we  are 
measuring  is  a  result  of  electrons  which  have  been  promoted 
in  the  conduction  band  of  the  GaAs  by  free-carrier  absorp¬ 
tion  to  energies  which  allow  them  to  migrate  over  the  barrier 
presented  by  the  AIAs.  We  have  argued  that  the  driving 
force  for  the  migration  of  the  electrons  through  the  AIAs 
barrier  is  the  difference  in  the  concentration  of  electrons  on 
either  side  of  the  barrier  which  is  a  result  of  the  difference  in 
light  intensity  on  either  side  of  the  barrier.  This  difference  in 
light  intensity  was  calculated  using  a  model  for  the  structure. 
These  calculations  explain  the  nature  of  the  voltage  spectra 
as  a  result  of  the  optical  properties  of  the  heterostructures,  in 
particular,  the  layer  thicknesses  and  doping  levels. 
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We  present  the  results  of  the  application  of  photoresponse  techniques  to  the  study 
of  the  transport  of  electrons  past  an  energy  barrier.  In  this  study,  the  barrier  was 
provided  by  a  thin  layer  of  ALAs  sandwiched  between  GaAs  layers.  The  experiment 
measures  the  voltage  resulting  from  the  flow  of  optically  excited  electrons  from 
one  side  of  the  barrier  to  the  other.  The  voltage  is  measured  as  a  function  of  the 
wavelength  of  the  incident  light.  We  also  present  the  results  of  calculations  which 
explain  the  nature  of  the  observed  spectra  and  how  they  change  when  the  thickness 
of  the  top  layer  (the  illuminated  side)  of  GaAs  is  changed  and  when  the  doping  in 
the  GaAs  is  varied. 


1.  Introduction 

In  this  paper  we  discuss  experimental  results  of  the 
application  of  photoresponse  techniques  to  the  study 
of  GaAs/AlAs  heterostructures.  We  have  applied  this 
technique  to  the  study  of  the  transport  of  electrons  past 
an  energy  barrier.  In  particular,  we  have  investigated 
these  electronic  transport  properties  in  structures  con¬ 
sisting  of  a  thin  layer  of  ALAs  sandwiched  between  layers 
of  GaAs.  The  experimental  procedure  consisted  of  illu¬ 
minating  samples  with  light  in  order  to  excite  electrons 
in  the  conduction  band  of  the  GaAs  which  could  then 
flow  past  the  ALAs  barrier.  We  then  measured  the  re¬ 
sulting  open  circuit  voltage  produced  across  the  AlAs 
barrier  as  a  function  of  the  wavelength  of  the  incident 
light. 

The  basic  results  presented  in  this  paper  can  be 
summarized  as  follows.  The  sign  of  the  voltage  that 
we  measured  was  consistent  with  electrons  flowing  from 
the  illuminated  side  of  the  AlAs  barrier  to  the  back  side 
of  the  AlAs  barrier.  The  magnitude  of  the  voltage  was 
proportional  to  the  intensity  of  the  illuminating  light. 
The  voltage  spectrum,  that  is,  the  voltage  as  a  func¬ 
tion  of  the  wavelength  of  the  incident  light,  had  two  ba¬ 
sic  shapes  which  were  common  to  all  the  AlAs  barrier 
samples  so  far  studied.  The  two  shapes  seemed  mainly 
to  depend  on  the  thickness  of  the  top  layer  of  GaAs. 
When  the  top  layer  was  on  the  order  of  a  few  microns 
thick,  the  resulting  spectrum  was  a  broad  peak  whose 
exact  position  depended  on  temperature.  When  the  top 
layer  was  thinner  (less  than  one  micron  thick)  the  volt¬ 
age  spectrum  rose  slowly  from  shorter  wavelengths  to 
longer  wavelengths  followed  by  a  sharp  turnoff  at  some 
particular  wavelength,  depending  on  the  temperature. 
The  peak  position  also  depended  on  the  doping  ievel 


in  the  GaAs  layers.  For  samples  with  a  higher  doping 
level,  the  peak  was  shifted  to  shorter  wavelengths.  The 
nature  of  the  spectra  is  explained  by  our  calculations  as 
will  be  shown  below. 

2.  Samples  and  Experimental  Methods 

The  samples  used  in  this  study  were  grown  by  a 
metalorganic  chemical  vapor  deposition  technique1,2 
with  all  the  AlAs  barriers  doped  p-type  with  Mg5,4 
while  the  GaAs  was  doped  n-type  with  Se.  The  GaAs 
substrates  on  which  these  structures  were  grown  were 
doped  n-type  at  about  3  x  1018  cm’5  with  Si.  We 
will  discuss  results  obtained  from  four  different  samples 
whose  characteristics  are  summarized  in  table  1.  The 


TABLE  1  Sample  Characteristics 


Sample 

Barrier 

Thickness 

GaAs 

Doping  (cm-5) 

Barrier 

Doping  (cm-5) 

1 

110  A 

Se  -  5  x  1018 

Mg  -  1  x  1018 

2 

160  A 

Se  -  5  x  1018 

Mg  -  1  x  1018 

3a 

160  A 

Se  -  3  x  1018 

Mg  -  lx  1018 

46 

240  A 

Se  -  3  x  1018 

Mg-  1  x  1018 

a  GaAs  doping  1.2  x  10I81  pm  either  side  of  the  AL4s 
barrier. 


b  GaAs  doping  gradually  dropped  to  2  x  1018  in  front 
of  the  AlAs  barrier. 

thicknesses  of  the  top  layer  (illuminated  side)  of  Gals 
varied  from  about  0.5  pm  thick  to  4.4  pm  thick.  Ohmic 
contacts  were  made  to  the  samples  by  evaporation  of  a 
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Au/Ge  alloy  followed  by  a  20  second  anneal  at  420°C.S 
Photolithographic  techniques  were  employed  to  make 
ring  shaped  contacts  to  the  top  surface  to  allow  light  to 
penetrate  into  the  GaAs  through  the  center  of  the  rings. 
The  individual  contacts  were  then  mesa  isolated  from 
each  other.  The  mesas  were  created  by  etching  with 
a  GaAs  etch  (4:1:1,  H2SC>4:H202  :H20)  to  a  depth  of 
about  10  pm.  The  samples  were  mounted  on  transistor 
headers,  and  wire  bonds  were  made  to  the  Au / Ge  con¬ 
tacts.  The  headers  could  then  be  mounted  in  a  Janis 
variable  temperature  dewar  (Janis  model  DT  11),  and 
thus  could  be  cooled  to  any  temperature  desired  from 
room  temperature  to  4.2K.  The  light  source  used  to  illu¬ 
minate  the  samples  was  a  tungsten  filament  lamp.  The 
light  was  directed  through  a  spectrometer  (SPEX  model 
1269)  and  then  focussed  onto  a  particular  device  to  be 
studied.  The  measurement  of  the  voltage  resulting  from 
the  illumination  was  made  directly,  using  a  voltmeter, 
when  it  was  desired  to  know  the  sign  of  the  voltage.  In 
genera],  however,  the  light  was  chopped  at  about  270 
Hz,  and  the  resulting  voltage  was  measured,  using  a 
lock-in  amplifier.  The  signal  was  then  digitized,  and 
the  spectrum  stored  on  a  signal  averager  (NIC  1170). 

3.  Theoretical  Model 

To  model  the  absorption  of  the  light  in  these  struc¬ 
tures,  we  carry  out  the  following  calculations.  We  con¬ 
sider  a  model  where  we  have  a  layer  of  GaAs,  a  layer 
of  A lAs,  and  a  semi-infinite  layer  of  GaAs.  Light  im¬ 
pinges  on  the  first  layer  of  GaAs  from  a  vacuum  at  nor¬ 
mal  incidence.  We  assume  that  no  attenuation  of  the 
light  takes  place  in  the  AlAs.  In  the  GaAs,  however,  we 
model  the  attenuation  of  the  light,  using  a  coefficient  of 
absorption,  o,  which  we  obtained  by  fitting  and  extrap¬ 
olating  the  curves  published  by  H.C.  Casey  et  al.6  for 
heavily  doped  n-type  GaAs  at  room  temperature.  We 
used  constant  values  for  the  index  of  refraction  of  GaAs 
and  AlAs  of  3.4  and  3.1,  respectively7.  The  light  inten¬ 
sity  was  then  calculated  throughout  the  model  structure 
and  was  integrated  in  the  region  of  GaAs  in  front  of  the 
AlAs  barrier  and  in  an  equal  region  behind  the  AlAs 
barrier  to  obtain  the  total  light  intensity  in  each  of  the 
two  regions.  The  total  light  intensity  behind  the  AlAs 
barrier  was  then  subtracted  from  the  total  in  front  of 
the  ALAs,  and  this  difference  was  plotted  as  a  function 
of  wavelength.  When  there  is  more  light  intensity  in 
front  of  the  AlAs  than  behind,  there  are  correspond¬ 
ingly  more  electrons  excited  there.  It  is  this  difference 
in  the  concentration  of  electrons  in  the  conduction  band 
which  drives  the  flow  of  electrons  from  the  front  side  of 
the  AlAs  barrier  to  the  back  side.  As  will  be  discussed 
below,  it  is  those  electrons  which  are  excited  in  the  con¬ 
duction  band  of  the  GaAs  by  free  carrier  absorption 
which  have  enough  energy  to  go  over  the  AlAs  barrier. 
The  number  of  electrons  excited  in  the  conduction  band 
due  to  free  carrier  absorption  is  proportional  to  the  light 
intensity.  Thus,  we  would  expect  the  measured  voltage 
spectrum  to  have  the  same  shape  as  the  calculated  dif¬ 
ference  in  the  light  intensity  in  front  of  and  behind  the 
AlAs  barrier.  Recent  calculations  by  C.  Mailhot  et  al.® 


6500  7000  7500  8000  8500  9000 

Wavelength  of  Incident  Light(Angstrom) 


Fig.  1.  Voltage  spectrum  taken  at  room  temperature 
for  a  sample  with  a  240  A  AlAs  barrier.  The  GaAs 
layers  were  doped  n-type  at  about  3  x  1018  cm~s.  The 
thickness  of  the  top  layer  of  GaAs  was  4.4  pm  and  0.9  pm, 
as  indicated  in  the  figure. 


have  shown  that  it  is  the  T-point  of  the  AlAs  which  de¬ 
termines  the  energy  barrier  presented  to  an  electron  in 
the  GaAs,  and  our  observations  so  far  have  been  consis¬ 
tent  with  this  conclusion  in  that  we  have  not  observed 
a  voltage  at  wavelengths  corresponding  to  energies  sub¬ 
stantially  below  this  value. 

4.  Results 

In  Fig.  1,  we  present  a  voltage  spectrum  for  sample 
4  at  room  temperature.  This  sample  had  a  top  layer  of 
GaAs  which  was  4.4  pm  thick  in  one  case  and  which  had 
been  etched  down  to  a  thickness  of  0.9  pm  in  the  other 
case.  We  note,  first,  that  in  both  cases  we  have  a  sub¬ 
stantial  amount  of  signal  at  energies  below  the  GaAs 
band  gap  (1.424  eV  or  8700  A  at  room  temperature) 
At  wavelengths  below  the  GaAs  band  gap,  no  holes  are 
created  in  the  structure  and,  hence,  the  photo-voltage 
is  not  a  result  of  hole  movement  in  the  GaAs.  At  wave¬ 
lengths  above  the  GaAs  band  gap  where  we  do  create 
electron  hole  pairs,  the  movement  of  electrons  and  holes 
in  the  GaAs  due  to  the  electric  fields  at  the  GaAs-AlAs 
interface  would  result  in  a  voltage  opposite  in  sign  to  the 
one  that  we  measure.  These  observations  strongly  sug¬ 
gest  that  it  is  electrons  promoted  in  the  GaAs  by  free 
carrier  absorption  to  energies  greater  than  the  GaAs- 
AIAs  offset  which  can  flow  past  this  barrier  and  produce 
our  observed  voltages.  To  explain  the  general  shape  of 
these  spectra,  we  refer  to  Fig.  2  which  schematically 
shows  the  band  structure  and  absorption  processes  in 
the  sample.  At  a  doping  level  of  3  x  1018  cm-5,  the 
fermi  level  can  be  as  much  as  110  meV  above  the  con¬ 
duction  band  minimum.  The  two  absorption  processes 
are  free  carrier  absorption  (A)  for  which  the  coefficient 
of  absorption  is  qa  and  band-to-band  absorption  (B) 
for  which  the  coefficient  of  absorption  is  05.  Since  free 
carrier  absorption  involves  a  phonon  and  band-to-band 
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Fig.  2.  Schematic  diagram  of  the  conduction  band  {£c) 
and  valence  band  (£„)  edges  of  the  GaAs/AlAs  het¬ 
erostructure  showing  the  fermi  level  (f/)  in  the  conduc¬ 
tion  band.  Indicated  in  the  figure  are  free  carrier  (A) 
and  band-to-band  (B)  absorption  of  light  by  electrons. 


Fig.  3.  Calculated  difference  in  the  integrated  light 
intensities  in  front  of  and  behind  the  AlAs  barrier  as  a 
function  of  wavelength.  The  coefficient  of  absorption,  a. 
used  in  calculating  these  curves  was  for  a  GaAs  n-type 
doping  of  3.3  x  1018  cm"5  at  room  temperature.  The 
results  are  presented  for  a  top  layer  thickness  of  GaAs 
of  4.4  p m  and  0.9  jim,  as  indicated  in  the  figure. 


absorption  does  not  aA  <  aB-  At  long  wavelengths 
where  the  energy  of  a  photon  is  insufficient  to  promote 
an  electron  from  the  valence  band  over  the  fermi  level, 
the  much  weaker  free  carrier  absorption  dominates,  and 
there  is  not  a  substantial  difference  between  the  light  in¬ 
tensity  in  front  of  and  behind  the  AlAs  barrier.  There 
is,  thus,  a  correspondingly  small  difference  in  the  num¬ 
ber  of  electrons  promoted  in  the  conduction  band  of 
the  GaAs,  and  we  expect  to  measure  a  small  signal. 
At  short  wavelengths  where  band-to-band  absorption 
is  very  strong,  the  light  intensity  is  attenuated  very 
quickly  near  the  surface  of  the  sample,  and  the  difference 
in  light  intensity  in  front  of  and  behind  the  ALAs  barrier 
is  very  small,  and,  again,  we  expect  to  measure  a  small 
signal.  Between  these  two  extremes,  we  have  some  wave¬ 
length  at  which  the  difference  in  light  intensity  (and, 
hence,  the  concentration  of  electrons  promoted  in  the 
conduction  band  of  the  GaAs)  is  a  maximum,  and  we 
expect  to  see  a  maximum  signal.  We,  thus,  have  a  spec¬ 
trum  with  a  peak  whose  position  is  mainly  determined 
by  the  doping  of  the  sample.  To  explain  the  difference 
in  the  spectra  when  we  have  a  thick  top  layer  (5  p m) 
and  a  thin  top  layer  (l  pm),  consider  the  following.  For 
a  sample  with  a  top  layer  thickness  of  5  pm  at  wave¬ 
lengths  where  process  B  cannot  occur  (the  energy  of 
the  light  is  insufficient  'o  promote  carriers  from  the  va¬ 
lence  band  over  the  fermi  level),  the  light  intensity  is 
proportional  to  e~iaA.  At  shorter  wavelengths  where 
process  B  can  occur,  the  intensity  of  the  light  is  pro¬ 
portional  to  e~iaB.  Consider  the  same  sample  with  a 
thin  top  layer  of  1  pm.  The  enhancement  in  the  light 
intensity  at  wavelengths  where  process  A  dominates  is 
proportional  to 

e~la* /e~ia*  =  e*°A, 


but  at  wavelengths  where  process  B  dominates  (on  the 
short  wavelength  tail  of  the  spectrum),  the  enhancement 
in  the  intensity  of  light  is  proportional  to 

e-iaB/e-So,B  _  e*oB 

but  since  aA  C  aB  then  e*aB  >  e*aA  and  the  enhance¬ 
ment  in  the  intensity  of  the  light  is  much  greater  at  the 
shorter  wavelength  tail  of  the  spectrum.  Since  the  num¬ 
ber  of  electrons  promoted  by  process  A  is  proportional 
to  the  light  intensity,  we  would  expect  the  voltage  signal 
to  be  most  enhanced  at  the  short  wavelength  tail  of  the 
Bpectrum.  Thus,  a  spectrum  for  a  sample  with  a  thick 
top  layer  will  change  qualitatively  in  the  way  it  does 
in  Fig.  1  when  the  top  layer  is  made  thinner.  The  T- 
point  barrier  at  room  temperature  is  at  most  1.35  eV9 
(9180  A).  This  barrier  is  decreased  by  the  fact  that  the 
fermi  level  is  in  the  conduction  band  but  is  increased 
by  the  effects  of  band  gap  shrinkage  in  the  GaAs  due  to 
the  heavy  doping10  and  is  also  affected  by  band  bend¬ 
ing  in  these  structures.  This  number,  however,  gives  us 
an  estimate  as  to  the  longest  wavelength  at  which  we 
could  expect  to  measure  a  voltage  if  there  is  also  a  dif¬ 
ference  in  the  light  intensity  in  front  of  and  behind  the 
ALAs  barrier.  At  wavelengths  longer  than  9180  A,  even 
electrons  promoted  via  process  A  may  not  be  energetic 
enough  to  go  over  the  barrier. 

In  Fig.  3,  we  present  results  of  the  predicted  spec¬ 
trum  for  a  sample  with  a  doping  level  of  3.3  x  1018  cm-3 
in  the  GaAs,  which  has  a  240  A  layer  of  AlAs  and  which 
has  a  top  layer  thickness  o'  4.4  pm  or  0.9  pm  as  la¬ 
beled.  We  see  that  since  the  absorption  coefficient  data 
contains  in  it  the  effects  of  band-to-band  versus  free 
carrier  absorption,  it  correctly  predicts  the  qualitative 
behaviour  of  the  spectra  we  measured.  We  note  that  the 
peak  position  (in  the  case  of  the  thick  top  layer)  and  the 
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Fig.  4.  Voltage  spectrum  for  samples  1  through  4  at 
80  K  illustrating  the  relative  peak  positions.  The  band 
gap  of  GaAs  at  80  K  is  about  1.51  eV  (corresponding 
to  about  8230  A);  samples  3  and  4  show  signal  at  wave¬ 
lengths  corresponding  to  energies  below  this  value. 


turnoff  wavelength  (in  the  case  of  the  thin  top  layer)  are 
not  exactly  predicted  by  the  calculations  though  they 
are  in  fairly  good  agreement.  These  positions  depend 
very  much  on  the  doping  of  the  sample  used  in  the  ex¬ 
periment.  In  the  calculations,  these  positions  depend  on 
the  behaviour  of  the  coefficient  of  absorption  a,  which, 
in  turn,  depends  on  the  doping  and  which  differs  from 
the  actual  doping  of  the  samples.  In  our  calculations, 
we  have  assumed  a  constant  doping  profile,  while  in  the 
actual  samples,  the  doping  in  the  GaAs  is  a  function  of 
position. 

In  Fig.  4,  we  present  the  voltage  spectrum  for  four 
different  samples  at  80  K.  We  note  that  there  are  dif¬ 
ferences  in  the  spectrum  from  sample  to  sample,  we  can 
nonetheless  identify  and  explain  the  relative  peak  po¬ 
sitions.  We  see  that  in  Fig.  4  the  peak  positions  for 
samples  1  and  2  are  the  same  and  that  the  peak  posi¬ 
tions  for  samples  3  and  4  are  the  same  though  shifted 
to  longer  wavelengths  than  for  1  and  2.  In  Fig.  5,  we 
present  the  results  of  the  calculations  for  the  light  in¬ 
tensity  difference  in  front  of  and  behind  the  AlAs  layer 
for  various  dopings  in  the  GaAs  regions  as  indicated. 
These  curves  would  correspond  to  spectra  obtained  at 
room  temperature  since  the  behaviour  of  a  was  taken 
from  data  obtained  at  room  temperature.  In  the  calcu¬ 
lations  we  have  kept  the  AlAs  layer  thickness  and  the 
GaAs  top  layer  thickness  constant.  This  is  done  in  order 
to  isolate  the  effects  of  the  doping  level.  We  note  that 
the  peak  positions  shift  to  shorter  wavelengths  as  the 
doping  in  the  GaAs  is  increased.  The  difference  in  the 
peak  position  is  about  200  A  when  the  doping  level  is 
changed  from  3.3  x  I01S  cm-*  to  6.7  x  1018  cm-5  in  the 
GaAs.  Since  these  calculations  are  performed  at  room 
temperature  and  the  experiment  at  80  K,  we  would  ex¬ 
pect  the  absolute  peak  positions  in  the  experiment  to 
be  at  shorter  wavelengths  than  in  the  calculations,  but 
we  would  expect  the  relative  peak  positions  for  similar 
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Fig.  5.  Calculated  difference  in  the  integrated  light 
intensities  in  front  of  and  behind  the  AlAs  barrier  as 
a  function  of  wavelength.  The  model  in  all  three  cases 
was  for  a  top  layer  of  GaAs  which  was  4.4  /nn  thick  and 
for  an  AlAs  layer  thickness  of  240  A.  The  coefficient  of 
absorption,  a,  used  in  calculating  these  curves  was  for  a 
GaAs  n-type  doping  of  2.0xl018  cm-*,  3.3xI018  cm-*, 
and  6.7  x  1018  cm-*,  as  indicated  in  the  figure. 


differences  in  doping  levels  to  be  the  same.  The  abso¬ 
lute  peak  position  is  a  function  of  temperature,  since 
for  lower  temperatures,  the  band  gap  of  the  GaAs  in¬ 
creases.  This,  in  turn,  means  that  at  lower  tempera¬ 
tures  the  absorption  edge  of  the  GaAs  shifts  to  shorter 
wavelengths.  We  see  that  the  experimental  peaks  are 
at  shorter  wavelengths  than  the  calculated  peaks  due  to 
the  temperature  difference  as  expected.  We  note  that 
the  peaks  in  the  spectra  for  samples  1  and  2  which  are 
both  doped  at  5  x  1018  cm-*  are  at  the  same  position 
and  are  at  shorter  wavelengths  than  the  peaks  in  the 
voltage  spectrum  for  samples  3  and  4.  This  is  consis¬ 
tent  with  the  fact  that  1  and  2  are  more  heavily  doped 
in  the  GaAs  region  than  samples  3  and  4.  The  peaks 
for  samples  3  and  4,  are  also  at  the  same  position  and 
this  too  is  consistent  with  the  fact  that  they  have  sim¬ 
ilar  doping  levels  in  the  GaAs  regions.  The  difference 
in  the  peak  positions  is  about  150  A,  and  this  is  con¬ 
sistent  with  the  predicted  difference  in  peak  positions 
when  one  considers  the  fact  that  the  difference  in  doping 
levels  used  in  the  calculations  is  larger  than  the  differ¬ 
ence  in  doping  levels  in  the  actual  samples.  Changing 
the  AlAs  layer  thickness  over  the  range  found  in  these 
samples  while  maintaining  a  constant  doping  level  does 
not  result  in  such  a  shift  in  the  peak  position.  In  the 
experimental  spectra  presented,  we  see  that  there  are 
undulations  which  are  not  as  marked  but  nonetheless 
visible  in  the  calculated  spectra  There  are  two  possible 
origins  for  these  undulations  in  the  experimental  data. 
One  explantion  is  that  they  are  the  result  of  the  optical 
properties  of  the  structures,  in  particular,  the  detailed 
behaviour  of  a.  In  our  calculations,  for  example,  we 
have  assumed  a  constant  doping  profile,  while  in  the 
actual  samples,  the  doping  in  the  Ga4s  is  a  function 
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of  position.  Another  possible  explanation  is  that  they 
are  the  result  of  resonant  effects  in  the  transmission  of 
electrons  over  the  AlAs  barrier.9 


5.  Summary 

In  summary,  we  have  presented  results  of  an  ex¬ 
perimental  technique  that  we  have  employed  to  study 
the  transport  of  electrons  past  an  energy  barrier.  We 
have  presented  experimental  results  which  imply  that 
the  voltage  we  are  measuring  is  a  result  of  electrons 
which  have  been  promoted  in  the  conduction  band  of 
GaAs  by  free  carrier  absorption  to  energies  which  allow 
them  to  flow  over  the  barrier  presented  by  the  AlAs. 
We  have  argued  that  the  driving  force  for  the  flow  of 
the  electrons  through  the  AlAs  barrier  is  the  difference 
in  the  concentration  of  excited  electrons  on  either  side 
of  the  barrier  which  is  a  result  of  the  difference  in  light 
intensity  on  either  side  of  the  barrier.  Thus,  two  condi¬ 
tions  must  be  satisfied  in  order  for  there  to  be  a  mea¬ 
sured  voltage.  The  wavelength  of  the  incident  light  must 
correspond  to  an  energy  larger  than  the  barrier  provided 
by  the  AlAs,  and  there  must  be  a  difference  in  the  light 
intensity  in  front  of  and  behind  the  AlAs  layer.  The  dif¬ 
ference  in  light  intensity  was  calculated  using  a  model 
for  the  structure.  We  have  also  presented  data  showing 
the  change  in  the  spectra  as  a  function  of  doping  level 
in  the  GaAs  layers.  We  presented  the  results  of  calcula¬ 
tions  which  explain  the  difference  in  peak  positions  for 
the  spectra  resulting  from  samples  with  different  doping 
levels  in  the  GaAs  and  different  AlAs  layer  thicknesses 
as  being  a  result  of  the  change  in  the  optical  absorp¬ 
tion  properties  of  the  structures  w’hich  result  from  the 
different  doping  levels. 
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We  present  the  results  of  the  application  of  photoresponse  techniques  to  the  study  of  the  transport  f 
electrons  past  an  energy  barrier.  In  this  study,  the  barrier  was  provided  by  a  thin  layer  of  AlAs 
sandwiched  between  GaAs  layers.  The  experiment  measures  the  voltage  resulting  from  the  flow  of 
optically  excited  electrons  from  one  side  of  the  barrier  to  the  other.  The  voltage  is  measured  as  a  function 
of  the  wavelength  of  the  incident  light.  We  present  results  obtained  for  samples  with  different  doping 
levels,  AlAs  layer  widths,  and  at  various  temperatures.  We  find  that  for  increasing  temperatures  the 
voltage  spectrum  shifts  to  longer  wavelengths,  and  that  for  higher  doping  levels  in  the  GaAs  the  voltage 
spectrum  shifts  to  shorter  wavelengths.  We  also  present  the  results  of  calculations  which  explain  why  the 
observed  spectra  change  as  a  function  of  temperature  and  as  the  doping  level  in  the  GaAs  is  varied. 


INTRODUCTION 

The  transport  properties  of  small  structures  is  currently 
an  area  of  very  active  research1'3.  The  precise  character  of 
the  barrier  and  its  influence  on  the  transport  is  a  key  aspect 
to  understanding  the  phenomena  governing  the  applications 
discussed  in  Refs.  1-3.  In  this  paper  we  discuss  basic  results 
obtained  by  photoresponse  techniques  applied  to  study 
GaAs/AIAs  single  barrier  heterostructures.  Photoresponse 
was  employed  to  study  the  transport  of  electrons  past  an 
energy  barrier  which  was  provided  by  a  thin  layer  of  AlAs 
sandwiched  between  layers  of  GaAs.  The  experimental  pro¬ 
cedure  consisted  of  illuminating  samples  with  light  to  excite 
electrons  in  the  conduction  band  of  the  GaAs  which  could 
then  flow  past  the  AlAs  barrier.  We  then  measured  the  re¬ 
sulting  open-circuit  voltage  produced  across  the  AlAs  bar¬ 
rier  as  a  function  of  the  wavelength  of  the  incident  light.  We 
present  results,  both  of  experiment  and  calculation  which 
describe  the  voltage  spectrum  obtained  as  a  function  of  tem¬ 
perature  and  for  samples  with  different  doping  levels  in  the 
GaAs. 

The  basic  results  presented  in  this  paper  can  be  sum¬ 
marized  as  follows.  The  sign  of  the  voltage  that  we  measured 
was  consistent  with  electrons  flowing  from  the  illuminated 
side  of  the  AlAs  barrier  to  the  back  side  of  the  AlAs  barrier. 
The  magnitude  of  the  voltage  was  proportional  to  the  inten¬ 
sity  of  the  illuminating  light  and  was  generally  greater  in  the 
samples  with  larger  zero-bias  impedances.  The  voltage  spec¬ 
tra,  that  is,  the  voltage  as  a  function  of  the  wavelength  of  the 
incident  light,  for  the  samples  discussed  in  this  paper  are 
broad  peaks.  The  reasons  for  this  have  been  discussed  pre¬ 
viously.4  In  this  paper  we  show  that  the  position  of  the  peak 
depends  on  the  sample  temperature  and  on  the  doping  in  the 
GaAs  cladding  layers.  It  was  found  that  for  lower  tempera¬ 
tures  or  higher  doping  levels  the  spectrum  shifted  to  shorter 
wavelengths.  Based  on  our  calculations  we  show  that  the 
shifts  in  the  spectra  are  explained  as  resulting  from  the 
change  in  the  optical  absorption  properties  of  these  struc¬ 
tures  which,  in  turn,  depend  on  the  temperature  and  doping, 
and  that  these  shifts  are  consistent  with  the  mechanism  that 
we  have  proposed  produces  the  observed  voltages. 


SAMPLES  AND  EXPERIMENTAL  METHODS 

The  samples  used  in  this  study  were  grown  by  a  meta- 
lorganic  chemical  vapor  deposition  technique,5  6  with  all  the 
AlAs  barriers  doped  p  type  with  Mg7  8  while  the  GaAs  was 
doped  n  type  with  Se.  The  substrates  on  which  these  struc¬ 
tures  were  grown  were  ( 100)-oriented  GaAs  wafers  doped  at 
about  3  x  1018  cm  “ 3  with  Si.  In  this  paper  we  will  discuss  the 
results  obtained  from  four  different  samples  whose  charac¬ 
teristics  are  summarized  in  Table  I.  For  these  doping  levels, 
the  AlAs  barrier  is  fully  depleted. 

In  Fig.  1  we  present  a  schematic  cross-sectional  view  of 
the  sample  geometry  which  is  discussed  below.  The  thick¬ 
nesses  of  the  top  layer  (illuminated  side)  of  GaAs  varied  from 
about  3.3  pm  thick  to  4.4  pm,  depending  on  the  sample, 
while  the  bottom  layer  was  substantially  thicker.  Ohmic 
contacts  were  made  to  the  samples  by  evaporation  of  a 
Au/Ge  alloy  followed  by  a  20-sec  anneal  at  420  ’C.9  Photo¬ 
lithographic  techniques  were  employed  to  make  ring-shaped 
contacts  to  the  top  surface  to  allow  light  to  penetrate  into  the 
GaAs  through  the  center  of  the  rings.  The  individual  con¬ 
tacts  were  then  mesa  isolated  from  each  other.  The  rings  had 
an  outer  diameter  of  about  300 pm  and  an  inner  diameter  of 
about  1 20  ^m  while  the  mesas  were  about  3  50  pm  in  diame¬ 
ter.  The  mesas  were  created  by  etching  with  a  GaAs  etch 
(4:1:1,  H2S04  :H202  :H20)  to  a  depth  of  about  10  pm.  The 
samples  were  mounted  on  transistor  headers,  and  wire  bonds 
were  made  to  the  Au/Ge  contacts.  The  headers  could  then 


TABLE  I  Sample  characteristics. 


Sample 

Barrier 
thickness  (A) 

GaAs  doping 

Se  ( x  10"  cm”’) 

Barrier  doping 

Mg  |  x  10"  cm"’) 

1 

110 

5 

I 

2 

160 

5 

1 

3* 

160 

3 

1 

4h 

240 

3 

1 

’GaAs  doping  1.2  x  10"  cm  '  I  pm  either  side  of  the  AlAs  barrier 
"GaAs  doping  gradually  dropped  to  2  X  10"  cm”’  in  front  of  the  AlAs 
barrier 
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FIG  1  Schematic  diagram  of  the  structure  studied.  The  mesas  are  approxi¬ 
mately  350  jim  in  diameter  and  about  10  /tm  high.  The  ring-shaped  top 
contacts  have  an  inner  diameter  of  about  1 20  |im.  The  experiment  consists 
of  measuring  the  voltage  between  the  top  and  bottom  contacts  as  a  function 
of  the  wavelength  of  the  incident  light. 

be  mounted  in  a  Janis  variable  temperature  dewar  (Janis 
model  DT  1 1 ),  and,  thus,  the  samples  could  be  cooled  to  any 
temperature  desired  from  room  temperature  to  4.2  K.  The 
light  source  used  to  illuminate  the  samples  was  a  tungsten 
filament  lamp.  The  light  was  directed  through  a  spectrom¬ 
eter  (SPEX  model  1269)  and  then  focused  onto  a  particular 
device  to  be  studied.  The  measurement  of  the  voltage  result¬ 
ing  from  the  illumination  was  made  directly,  using  a  volt 
meter,  when  it  was  desired  to  know  the  sign  of  the  voltage.  In 
general,  however,  the  light  was  chopped  at  about  270  Hz, 
and  the  resulting  voltage  was  measured  using  a  lock-in  am¬ 
plifier.  The  signal  was  then  digitized,  and  the  spectrum 
stored  on  a  signal  averager  (NIC  1170).  In  either  case,  the 
measurement  was  an  open  circuit  measurement  of  the  vol¬ 
tage  across  the  sample. 


THEORETICAL  MODEL 

To  model  the  absorption  of  the  light  in  these  structures, 
we  carried  out  the  following  calculations.  We  consider  a  one 
dimensional  model  where  we  have  a  layer  of  GaAs,  a  layer  of 
AlAs,  and  a  semi-infinite  layer  of  GaAs.  Light  impinges  on 
the  first  layer  of  GaAs  from  a  vacuum  at  normal  incidence. 
We  assume  that  no  attenuation  of  the  light  takes  place  in  the 
AlAs.  In  the  GaAs,  however,  we  model  the  attenuation  of 
the  light,  using  a  coefficient  of  absorption  a,  which  we  ob¬ 
tained  by  fitting  and  extrapolating  the  curves  published  by 
H.C.  Casey  et  a/.10  for  heavily  doped  n-type  GaAs  at  room 
temperature.  We  used,  however,  constant  values  for  the  in¬ 
dex  of  refraction  of  GaAs  and  AlAs  of  3.4  and  3.1,  respec¬ 
tively."  For  an  electromagnetic  wave  incident  onto  the  in¬ 
terface  between  two  dielectric  media  the  boundary  con¬ 
ditions  at  the  interface  are  the  continuity  of  the  tangential 
components  of  the  electric  and  magnetic  fields  E  and  H, 
respectively  If  we  have  a  wave  propogating  in  the  z  direction 
normal  to  the  interface  then  the  electric  field  is  given  by, 

E  =  E0  exp [i(eot  —  kz)]x  , 

where 


2 

and  where  n  is  the  real  part  of  the  index  of  refraction  and  a  is 
the  absorption  coefficient.  From  Maxwell’s  equation 

VxE+^i  —  =  0 
dt 

we  obtain  that  for  nonmagnetic  materials  the  continuity  of 
the  tangential  component  of  H  is  equivalent  to  the  contin¬ 
uity  of  dE  /dz.  In  matching  these  boundary  conditions  at  the 
interface  we  take  into  account  the  total  electric  field  due  to 
the  incident,  reflected,  and  transmitted  waves  in  any  region. 
The  computer  programs  employed  a  transfer  matrix  tech¬ 
nique  to  solve  for  the  electric  field  throughout  the  model 
structure.  This  allows  us,  in  principle,  to  extend  the  calcula¬ 
tion  to  more  layers.  The  light  intensity  was  then  calculated 
as  the  square  of  the  electric  field  throughout  the  model  struc¬ 
ture  and  was  then  integrated  in  the  region  of  GaAs  in  front  of 
the  AlAs  barrier  and  in  an  equal  region  behind  the  AlAs 
barrier  to  obtain  the  total  light  intensity  in  each  of  the  two 
regions.  The  quantitative  results  obtained  from  these  calcu¬ 
lations  (i.e.,  exact  peak  positions  and  symmetries)  depend 
somewhat  on  the  size  of  the  regions  on  either  side  of  the  AlAs 
barrier  one  chooses  to  integrate  over,  however,  this  does  not 
affect  the  qualitative  nature  of  the  results.  The  integration 
also  included  an  integration  over  the  phase  of  the  electric 
field  in  order  to  account  for  the  incoherence  of  the  light.  The 
total  light  intensity  behind  the  AlAs  barrier  was  then  sub¬ 
tracted  from  the  total  in  front  of  the  AlAs  and  this  difference 
was  plotted  as  a  function  of  wavelength.  At  present  the  cal¬ 
culations  assume  that  the  absorption  properties  of  the  GaAs 
are  those  of  bulk  GaAs  at  the  same  doping  levels.  This  is  not 
completely  true  at  the  GaAs/ AlAs  interface  where  the  band 
bending  can  drastically  alter  the  free-carrier  concentration, 
and  thus,  the  absorption  properties  of  the  GaAs  layers. 

RESULTS 

In  Fig.  2,  we  present  the  voltage  spectrum  for  sample  2 
with  a  4.0-/um  thick  top  layer  of  GaAs  for  various  tempera¬ 
tures.  We  note  that  there  is  measurable  signal  at  energies 


Wovelengih  of  Incident  Light  (Angstrom) 


FIG.  2.  Voltage  spectrum  for  sample  2  at  the  temperatures  indicated.  The 
GaAs  layers  were  doped  n-type  at  5x  10"  cm  “ 5  with  the  top  layer  of  GaAs 
4  0  jim  thick  The  AlAs  layer  was  doped  p  type  at  1  X  10"  cm  “  ’  and  was 
about  160  A  thick. 
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below  the  GaAs  band  gap  (1.52  eV  or  8160  A  at  6  K).  The 
heavy  doping  causes  the  band  gap  to  shrink.1*  but  it  also 
pushes  the  Fermi  level  into  the  conduction  band.  For  an  n- 
type  doping  of  5 x  10IS  cm"'  in  the  GaAs  we  estimate  that 
the  Fermi  level  is  160  meV  above  the  conduction  band  mini¬ 
mum.  The  net  effect  is  to  increase  the  energy  required  to  take 
an  electron  from  the  valence  band  into  the  conduction  band 
and  thus  increase  the  energy  required  to  produce  holes  in 
this  structure.  As  will  be  presented  below,  samples  3  and  4 
show  signal  at  even  longer  wavelengths  than  samples  1  and  2 
at  a  given  temperature.  At  wavelengths  below  the  GaAs 
band  gap,  no  holes  are  created  in  the  structure,  and  hence, 
the  photovoltage  is  not  the  result  of  hole  movement  in  the 
GaAs.  At  wavelengths  which  correspond  to  energies  that 
could  promote  an  electron  from  the  valence  band  into  the 
conduction  band  and,  thus,  also  create  a  hole,  band  bending 
at  the  GaAs/AlAs  interface  would  cause  electron  and  hole 
movement  which  would  create  a  voltage  opposite  in  sign  to 
the  one  that  we  measure.  We  also  note  that  the  sign  of  the 
voltage  is  consistent  with  electrons  flowing  from  the  illumi¬ 
nated  side  of  the  AlAs  barrier  to  the  back  of  the  AlAs  bar¬ 
rier.  On  the  basis  of  these  observations  we  co>  elude  that 
electrons  promoted  from  the  occupied  states  in  the  conduc¬ 
tion  band  of  the  GaAs  by  free-carrier  absorption  to  energies 
greater  than  the  GaAs/AlAs  conduction  band  offset  flow 
past  this  barrier  and  produce  the  observed  voltages.  The  dif¬ 
ference  in  light  intensity  in  front  of  and  behind  the  AlAs 
barrier  creates  a  difference  in  the  concentration  of  optically 
excited  electrons  on  either  side  of  the  AlAs  barrier.  It  is  this 
difference  in  the  concentration  of  optically  excited  electrons 
which  drives  the  electron  flow  from  the  front  to  the  back  of 
the  AlAs  barrier  and  which,  in  turn,  produces  the  voltage 
that  we  measure.  Since  the  number  of  electrons  excited  in  the 
conduction  band  due  to  free-carrier  absorption  is  propor¬ 
tional  to  the  light  intensity,  we  would  expect  the  measured 
voltage  spectrum  to  have  the  same  shape  as  the  calculated 
difference  in  the  light  intensity  in  front  of  and  behind  the 
AlAs  barrier.  Recent  calculations  by  C.  Mailhotef  a/.'3  have 
shown  that  it  is  the  /'-point  offset  between  the  AlAs  and  the 
GaAs  which  determines  the  energy  barrier  presented  to  an 
electron  in  the  GaAs.  The  energy  required  to  take  an  elec¬ 
tron  over  this  barrier  is  less  than  the  conduction  band  offset 
due  to  the  fact  that  the  Fermi  level  is  above  the  conduction 
band  minimum  but  is  increased  by  the  effects  of  band-gap 
shrinkage  in  the  GaAs  due  to  the  heavy  doping,  and  is  also 
affected  by  band  bending  in  these  structures.  The  energy 
required  to  promote  an  electron  over  this  barrier  also  de¬ 
pends  on  the  relative  change  in  the  band  gaps  of  GaAs  and 
AlAs  as  the  temperature  is  varied.  There  is  a  controversy 
regarding  the  value  of  the  conduction-band  energy-gap  dis¬ 
continuity  between  the  GaAs  and  the  AlAs. 14  At  room  tem¬ 
perature  if  the  offset  is  85%  of  the  total  energy-gap  discon¬ 
tinuity,  the  barrier  presented  to  an  electron  would  be  about 
1.35  eV.  For  a  smaller  value  of  the  offset  the  barrier  would  be 
substantially  less  than  this.  One  may  then  ask  why  we  do  not 
measure  an  appreciable  signal  at  long  wavelengths  that  still 
correspond  to  energies  greater  than  this  barrier0  We  do  not 
measure  any  signal  when  there  is  no  difference  in  the  number 
of  optically  excited  electrons  in  front  of  and  behind  the  AlAs 
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barrier.  Thus  we  measure  a  signal  as  long  as  there  is  a  differ¬ 
ence  in  the  light  intensity  in  front  of  and  behind  the  AlAs 
barrier.  The  wavelength  at  which  there  ceases  to  be  such  a 
difference  may  correspond  to  an  energy  that  is  still  substan¬ 
tially  greater  that  the  barrier  presented  to  an  electron  due  to 
the  conduction  band  offset  between  the  AlAs  and  the  GaAs. 
Thus,  the  cutoff  in  our  spectra  at  both  long  and  short  wave¬ 
lengths  is  determined  by  the  optical  properties  of  the  struc¬ 
ture  and  not  the  conduction-band  offset. 

We  can  explain  why  in  Fig.  2  the  peak  position  shifts  to 
longer  wavelengths  with  increasing  temperature.  It  is  band- 
to-band  and  free-carrier  absorption  which  determine  the  na¬ 
ture  of  the  optical  absorption  in  the  GaAs  and  which,  in 
turn,  determine  the  peak  position.4  Thus,  as  temperature 
increases,  the  shift  of  the  peak  position  to  longer  wave¬ 
lengths  is  a  reflection  of  the  shift  to  lower  energies  of  the 
band-to-band  absorption  edge  of  the  GaAs  for  increasing 
temperature.  The  band  gap  of  GaAs,  as  a  function  of  tem¬ 
perature,  can  be  approximated  by15 

Eg  =  1.519  -  5.405  X  10-JrV(204+  T)eV . 

To  see  whether  the  shift  in  the  peak  position  follows  the 
change  in  band  gap,  we  proceed  as  follows.  We  choose  the 
half  maximum  position  on  the  shorter  wavelength  edge  of 
the  peak  as  an  indication  of  the  peak  position.  We  then  plot 
the  change  in  the  peak  position  referenced  to  its  position  at  6 
K  versus  the  change  in  the  band  gap  referenced  to  its  value  at 
6  K.  If  then  the  shift  in  the  peak  position  is  indeed  a  reflec¬ 
tion  of  the  change  in  the  band  gap,  we  expect  to  have  a 
straight  line  with  slope  of  one  passing  through  the  origin  of 
our  plot.  In  Fig.  3,  we  present  this  plot.  We  see  that  we  do 
have  the  behavior  we  expect  though  the  band  gap  tends  to 
initially  change  slightly  faster  than  the  peak  position.  This 
difference  between  the  movement  in  the  peak  position  and 
the  change  in  band  gap  results  from  two  sources.  The  first  is 
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FIG  3  Change  in  the  band  gap  vs  ihe  shift  in  the  peak  position  in  the  vol¬ 
tage  spectrum  of  sample  2  as  a  function  of  temperature  The  straight  line 
indicates  the  values  for  the  change  in  the  band  gap  equal  to  the  shift  in  the 
peak  position. 
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the  inaccuracy  in  assigning  a  value  to  the  peak  position  at 
different  temperatures.  The  second  is  the  fact  that  we  have 
not  included  the  effect  on  the  absorption  edge  produced  by 
the  change  in  the  distribution  of  the  electrons  in  the  conduc¬ 
tion  band  of  the  GaAs  as  a  function  of  temperature. 

In  Fig.  4,  we  present  the  voltage  spectrum  for  four  dif¬ 
ferent  samples  at  80  K.  We  note  that  while  there  are  differ¬ 
ences  in  the  spectrum  from  sample  to  sample,  we  can  none¬ 
theless  identify  and  explain  the  relative  peak  positions.  We 
see  that  in  Fig.  4  the  peak  positions  for  samples  1  and  2  are 
the  same  and  that  the  peak  positions  for  samples  3  and  4  are 
the  same  though  shifted  to  longer  wavelengths  than  for  1  and 
2.  In  Fig.  5,  we  present  the  results  of  the  calculations  for  the 
light  intensity  difference  in  front  of  and  behind  the  AlAs 
layer  corresponding  to  various  dopings  in  the  GaAs  regions 
as  indicated.  These  curves  would  correspond  to  spectra  ob¬ 
tained  at  room  temperature  since  the  behavior  of  a  was  tak¬ 
en  from  data  obtained  at  room  temperature.  We  are  restrict¬ 
ed  to  carrying  out  our  calculations  at  room  temperature  and 
for  these  particular  doping  levels  in  the  GaAs  since,  to  our 
knowledge,  the  only  data  available  for  the  behavior  of  a  in 
heavily  doped  GaAs  is  that  which  we  obtained  from  Ref.  10. 
In  the  calculations,  we  have  kept  the  AlAs  layer  thickness 
and  the  GaAs  top  layer  thickness  constant.  This  is  done  in 
order  to  isolate  the  effects  of  the  doping  level.  We  note  that 
the  peak  positions  shift  to  shorter  wavelengths  as  the  doping 
in  the  GaAs  is  increased.  The  difference  in  the  peak  position 
is  about  180  A  when  the  doping  level  is  changed  from 
3.3  X  1018  cm~3to6.7x  1018  cm-3  in  the  GaAs.  Since  these 
calculations  are  performed  at  room  temperature  and  the  ex¬ 
periment  at  80  K,  we  would  expect  the  absolute  peak  posi¬ 
tions  in  the  experiment  to  be  at  shorter  wavelengths  than  in 
the  calculations,  but  we  would  expect  the  relative  peak  posi¬ 
tions  for  similar  differences  in  doping  levels  to  be  the  same. 
We  see  that  in  Fig.  4  the  experimental  peaks  are  at  shorter 
wavelengths  than  the  calculated  peaks  (Fig.  5)  due  to  the 
temperature  difference  as  expected.  We  note  that  the  peaks 
in  the  spectra  for  samples  1  and  2  which  are  both  doped  at 
5x  I0'8  cm-3  are  at  the  same  position  and  are  at  shorter 
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FIG.  4.  Voltage  spectrum  for  samples  1  through  4  at  80  K,  illustrating  the 
relative  peak  positions.  The  band  gap  of  GaAs  at  80  K  is  aboui  151  eV 
(corresponding  to  about  8230  A),  samples  3  and  4  show  strong  signal  at 
wavelengths  corresponding  to  energies  below  this  value. 


wavelengths  than  the  peaks  in  the  voltage  spectrum  for  sam¬ 
ples  3  and  4.  This  is  consistent  with  the  fact  that  1  and  2  are 
more  heavily  doped  in  the  GaAs  region  than  samples  3  and 
4.  The  peaks  for  samples  3  and  4  are  also  at  the  same  posi¬ 
tion,  and  this,  too,  is  consistent  with  the  fact  that  they  have 
similar  doping  levels  in  the  GaAs  regions.  The  difference  in 
the  peak  positions  is  about  125  A,  and  this  is  consistent  with 
the  predicted  difference  in  peak  positions  when  one  consid¬ 
ers  the  fact  that  the  difference  in  doping  levels  used  in  the 
calculations  is  larger  than  the  difference  in  doping  levels  in 
the  actual  samples.  Changing  the  AlAs  layer  thickness  over 
the  range  found  in  these  samples,  while  maintaining  a  con¬ 
stant  doping  level,  does  not  result  in  such  a  shift  in  the  peak 
position. 

In  the  experimental  spectra  presented,  we  see  that  there 
are  undulations  which  are  not  as  marked  but  nonetheless 
visible  in  the  calculated  spectra.  There  are  two  possible  ori¬ 
gins  for  these  undulations.  One  explanation  is  that  they  are 
the  result  of  the  optical  properties  of  the  structures,  in  parti¬ 
cular,  the  detailed  behavior  of  a.  In  our  calculations,  for 
example,  we  have  assumed  a  constant  doping  profile,  while 
in  the  actual  samples  the  doping  in  the  GaAs  is  a  function  of 
position.  Another  possible  explanation  is  that  they  are  a  re¬ 
sult  of  resonant  effects  in  the  transmission  of  electrons  over 
the  AlAs  barrier.  The  peak  position  is  constant  as  we  move 
from  one  device  to  another  across  the  surface  of  a  given  wa¬ 
fer.  The  undulations  are  reproduced  from  device  to  device 
across  a  given  wafer  in  terms  of  their  amplitudes  but  not 
necessarily  in  terms  of  their  position.  Either  origin  for  the 
undulations  would  be  very  sensitive  to  small  variations  in  t  lie 
properties  of  the  structure  from  one  position  to  another  on  a 
wafer.  In  the  case  of  the  optical  origin,  any  variation  of  the 
doping  profile  or  layer  thicknesses  from  point  to  point  would 
affect  the  undulations,  while  any  resonant  effects  would  be 
very  sensitive  to  the  exact  thickness  of  the  AlAs  barrier. 

Larger  voltages  are  generally  observed  for  samples  with 
larger  zero-bias  impedances.  By  zero-bias  impedance  we 
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FIG.  5.  Calculated  difference  in  the  integrated  light  intensities  in  front  of 
and  behind  the  AlAs  barrier  as  a  function  of  wavelength.  The  model  in  all 
three  cases  was  for  a  top  layer  of  GaAs  which  was  4.4  /um  thick  and  for  an 
AlAs  layer  thickness  of 240  A  The  coefficient  of  absorption  a  used  in  calcu¬ 
lating  these  curves  was  for  a  GaAs  n-type  doping  of  2. Ox  10"  cm 
3.3  X  10"  cm  "  \  and  6.7x  10"  cm  ~ J,  as  indicated  in  the  figure. 
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using  a  model  for  the  structure.  We  have  presented  data 
showing  the  behavior  of  the  voltage  spectra  as  a  function  of 
temperature  which  we  have  explained  by  considenng  the 
change  in  the  band  gap  of  GaAs.  We  have  also  presented 
data  showing  the  change  in  the  spectra  as  a  function  of  dop¬ 
ing  level  in  the  GaAs  layers.  We  presented  the  results  of 
calculations  which  explain  the  difference  in  peak  positions 
for  the  spectra,  resulting  from  samples  with  different  doping 
levels  in  the  GaAs,  and  different  AlAs  layer  thicknesses,  as 
being  a  result  of  the  change  in  the  optical  absorption  proper¬ 
ties  of  the  structures  which  result  from  the  different  doping 
levels. 
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but  we  do  not  have  as  yet  a  satisfactory  explanation  for  this 

SUMMARY 

In  summary,  we  have  presented  results  of  an  expen- 
mental  technique  that  we  have  employed  to  study  the  trans¬ 
port  of  electrons  past  an  energy  barrier.  We  have  presented 
experimental  results  which  imply  that  the  voltage  we  are 
measuring  is  a  result  of  electrons  which  have  been  promoted 
in  the  conduction  band  of  GaAs  by  free-carrier  absorption  to 
energies  which  allow  them  to  flow  over  the  barrier  presented 
by  the  AlAs.  We  have  argued  that  the  driving  force  for  the 
flow  of  the  electrons  through  the  AlAs  barrier  is  the  differ¬ 
ence  in  the  concentration  of  excited  electrons  on  either  side 
of  the  barrier  which  is  a  result  of  the  difference  in  light  inten¬ 
sity  on  either  side  of  the  barrier.  Thus,  two  conditions  must 
be  satisfied  for  there  to  be  a  measured  voltage.  The  wave¬ 
length  of  the  incident  light  must  correspond  to  an  energy 
larger  than  the  barrier  provided  by  the  AlAs,  and  there  must 
be  a  difference  in  the  light  intensity  in  front  of  and  behind  the 
AlAs  layer.  The  difference  in  light  intensity  was  calculated, 
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We  present  a  model  of  the  metal-semiconductor  junction,  for  heavily  doped  GaAs,  so  that 
tunneling  dominates  the  current.  It  is  assumed  that  the  imaginary  part  of  the  wave  vector  in  the 
semiconductor  is  given  by  the  two-band  model.  Modifications  in  the  barrier  potential  due  to 
image  charge,  negative  charge  near  the  interface,  and  the  degenerate  doping  of  the  semiconductor 
are  included  The  role  of  the  L-point  minimum  in  the  GaAs  in  determining  the  position  of  the 
Fermi  level  in  the  semiconductor  is  included.  The  energy  distribution  of  the  conductance  as  a 
function  of  doping  and  barrier  height  is  given.  The  contact  resistance  as  a  function  of  doping  and 
barrier  height  is  also  presented.  The  results  suggest  that  previous  calculations  are  substantially  in 
error  due  to  the  simple  models  that  were  used  for  the  dependence  of  the  imaginary  part  of  the  wave 
vector  on  energy. 


I.  INTRODUCTION 


II.  THEORY 


The  development  of  GaAs  technology  has  led  to  a  continu¬ 
ing  search  for  methods  of  making  Ohmic  contacts  to  a -type 
GaAs.  One  of  the  very  interesting  recent  developments  has 
been  the  growth  of  heavily  doped  layers  using  molecular 
beam  epitaxy  (MBE].1  Furthermore,  the  most  widely  used 
theory  of  Ohmic  contacts  is  that  due  to  Cljang,  Fang,  and 
Sze  (CFS),’  in  which  a  metal-semiconductor  contact  is  treat¬ 
ed.  Attempts  to  compare  theory  with  experiment  have  been 
mainly  based  on  this  theory.3  4  5 

In  CFS’s  theory,  they  treat  the  tunneling  of  carriers 
through  a  Schottky  barrier  as  a  function  of  doping  in  the 
semiconductor.  The  theory  includes  the  correction  of  the 
image  potential  to  the  barrier.  It  assumes  that  the  tunneling 
through  the  barrier  is  characterized  by  an  imaginary  wave 
'ector  vs  energy  that  is  simple  parabolic  (one-band  I  relation¬ 
ship.  Given  the  current  interest  in  the  theory,  it  is  desirable 
to  improve  the  theory  so  that  it  treats  the  best  model  avail¬ 
able.  With  these  improved  results,  one  might  hope  for  a  more 
realistic  comparison  between  the  theory  and  experiment. 

In  this  paper,  we  present  the  results  of  a  more  realistic 
treatment  of  the  contact  resistance  for  electrons  tunneling 
through  the  Schottky  barrier  between  a  metal  and  GaAs, 
which  is  heavily  doped  n  type  The  theory  includes  the  ef¬ 
fects  of  image  charge,  negative  charge  at  the  interface,  and 
the  presence  of  the  4-point  minimum  in  the  GaAs.  The  big¬ 
gest  correction  is  due  to  the  improvement  in  the  model  of  the 
imaginary  wave  vector  vs  energy  relation.  We  assume  a  two- 
band  model  that  uses  the  correct  energy  gap  and  conduction 
band  effective  mass  for  GaAs.  We  find  results  for  the  contact 
resistance  which  differ  by  as  much  as  an  order  of  magnitude 
in  the  doping  ranges  of  interest  for  devices. 

This  paper  is  organized  in  the  following  way  Section  II 
presents  the  theoretical  model.  In  Sec.  Ill  we  give  the  results 
for  the  model,  and  compare  them  with  CFS.  In  Sec.  IV,  we 
show  an  example  of  fluctuations  in  the  barrier  potential  due 
to  the  discrete  nature  of  the  doping.  Finally,  in  Sec.  V  we 
present  our  conclusions. 


For  a  metal-n-semiconductor  junction  in  reverse  bias, 
with  the  semiconductor  being  degenerate,  the  energy  dia¬ 
gram  is  shown  in  Fig.  1. 

The  corrections  due  to  the  image  force,6  negative  charge 
at  the  interface,7  and  nonparabolic  corrections  to  the  poten¬ 
tial  due  to  the  presence  of  conduction  electrons8  in  pan  of 
the  depletion  region  are  all  included.  The  metal  is  modeled  as 
a  degenerate  electron  gas,  with  the  radius  of  the  Fermi 
sphere  (the  Fermi  energy)  being  ~  7-8  eV.  (This  is  not  shown 
in  Fig.  1,  as  it  would  be  offscale.) 


Potential  at  a  metal  — 
semiconductor  Interface 


conduction 


band 
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Fh.  I  Diagram  of  a  metal-n  semiconductor  junction  for  an  applied  fuas  > 
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In  Fig.  1,  the  energy  is  measured  upwards  from  the  con¬ 
duction  band  edge  far  from  the  depletion  region;  ,  £FM 
are  the  Fermi  levels  in  the  semiconductor  and  metal,  respec¬ 
tively,  and  V  is  the  applied  bias.  Then,  following  Ref.  9, 
m,k„Tq  /* 

J  =  —  —  --  dET[E) 

lirff  Jo 


X  In 


exp[  -l E  -  £j.s )/kHT  ]  +  1 


(1) 


,exp[  —  {E  —  EFM  )/kgT  \  +  1 

where  E  is  the  energy  of  electron  normal  to  interface,  mc  is 
the  conduction  band  effective  mass  at  the  /’point,  and  T(E )  is 
the  transmission  probability.  Equation  (1)  neglects  phonon- 
assisted  processes,  i.e.,  the  transverse  wave  vector  k  is  as¬ 
sumed  to  be  conserved  for  tunneling  particle.  The  large  Fer¬ 
mi  radius  of  the  metal  in  reciprocal  space,  relative  to  the 
distribution  of  free  carriers  in  reciprocal  space  for  doped 
GaAs,  or  indeed  most  semiconductors,  makes  this  a  good 
approximation. 10  In  the  derivation  of  Eq.  (1),  a  number  of 
other  approximations  have  been  made.  These  approxima¬ 
tions  are  detailed  in  Ref.  11.  The  approximations  are  very 
commonly  made  in  treating  a  Schottky  barrier  and  are  not 
likely  to  result  in  large  errors.  One  major  assumption  is  used 
to  guarantee  that  zero  current  will  result  when  zero  voltage 
is  applied,  by  assuming  that  the  boundary  conditions  on  the 
current  from  the  semiconductor  to  the  metal  are  appropriate 
for  the  current  in  the  other  direction. 

By  the  WKB  approximation,  take 


T[E)  = 


exp^  —  2/J  klxjdx'J,  ifO<,E<qVm 


(2) 


J.  E>^m 

where  x,  and  x2  are  the  turning  points  of  the  forbidden  re¬ 
gion  for  k.  Assuming  a  direct-gap  semiconductor,  a  two- 
band  model,  from  k-p  theory,  gives  for  k  everywhere  within 
the  forbidden  region,12 


k  =  -  i 


2mc 

ZE 


j_=a/=  r 
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Rc  is  the  specific  contact  resistance  of  the  interface,  mea¬ 
sured  in  fl  cm2,  and  G  {E  |  is  the  conductance  distribution 
function.  (Often,  in  the  literature,  Rc  is  defined  at  V  =  0,  but 
there  is  no  need  to  confine  ourselves  to  this.)  Then, 
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Rc  and  G(E)  are  found  from  Eq  1 5 )  by  numerical  methods 
[Note  that  at  zero  bias  the  second  term  in  Eq  |5|  is  0  ]  To 


obtain  EFS  in  Eq.  (5),  n,  the  free  carrier  concentration,  is 
compared  with W,-,  the  effective  density  of  states  in  the  con¬ 
duction  band.  We  assume  that  T  is  sufficiently  large  to  ne¬ 
glect  carrier  freezeout.  Then 


£Fs  — 


(kB  /In  (n/Wc) 
fr 


2m, 


-{Itth) 


n  <.VC, 
n>Nc. 


For  carrier  concentrations  in  the  GaAs  sufficiently  large 
such  that  the  /,-point  states  are  occupied,  the  semiconductor 
Fermi  level  is  given  by 


n  = 


1 

3  v2 


where  mL  is  the  effective  mass  at  the  L  point,  ana  Erl  is  the 
energy  separation  between  the  L  point  and  the  r  point.  The 
factor  of  4  in  the  second  term  on  right-hand  side  of  Eq.  ~ 
arises  because  each  Brillouin  zone  has  a  net  of  four  L  points 
The  potential  energy  of  the  electron  is  given  by 


d>(x}=  - 


r 


+  <t>Ax), 


1 8 ; 


167rex 

where,  for  degenerate  doping,  satisfies  the  implicit  rela¬ 
tionships 
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For  nondegenerate  doping,  where  E ps  <  0,  let  ne  =0.  Equa¬ 
tions  (9)  and  (10)  are  solved  for  <b ,  by  the  method  of  successive 
integration.13  In  doing  this,  the  boundary  condition  of 
^i(x)  =  0  at  x  =  1.2  X  the  depletion  length  is  chosen  to  en¬ 
sure  convergence.  The  complete  solution  for  0,  would  tend 
to  0  asymptotically  as  x— *oo.  Thus  the  boundary  condition 
lowers  <b ,  near  the  edge  of  the  depletion  region.  However, 
this  occurs  at  large  tunneling  distances,  relative  to  the  tun¬ 
neling  distances  at  higher  energies,  and  so  will  have  little 
effect  on  the  results. 

The  effect  of  negatively  charged  surface  states  can  be  in¬ 
cluded  by  adding  the  following  term  to  the  right-hand  side  of 

Eq.  (8): 

dqlN 


■  exp(  —  x/d  ), 


till 


where  .V  is  the  area  density  of  surface  states  and  d  is  the 
penetration  length  of  the  states. 

III.  RESULTS  AND  DISCUSSION 

We  have  used  the  following  values  for  the  parameters  in 
the  theory:  T  =  300  K.  Ec  =  1.42  eV,  ErL  =  0.284  fV, 
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mc  =  0.063,  mL  =  0.55,  *  =  12.85,  and  Nc  =4.21  X  10'7 
cm’’.  The  value  of  Ec  is  taken  from  Ref.  14,  while  the  other 
values  are  taken  from  Ref.  15. 

To  illustrate  the  contributions  to  the  conductance  per  unit 
area  per  unit  energy,  for  various  doping  concentrations,  we 
have  plotted  in  Fig.  2  the  barrier  shapes  and  values  of  G(E\ 


Potential  and  Conductance 
at  several  dopings 
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1  O20 
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Fig  2  Distributions  of  the  specific  contact  conductance  as  a  function  of 
energy  To  illustrate  the  energy  position  on  the  barrier,  the  barrier  shape  is 
plotted  using  the  same  energy  scale  The  log  of  the  conductance  is  shown, 
with  the  conductance  being  measured  in  fl  1  m  !  J  1 


for  four  different  dopings:  1017,  10'“,  1019,  and  102°  cm"  '. 
These  are  at  <t>Bn  =  0.8  eV,  T  =  300  K,  and  V  =  0  V.  Nega¬ 
tive  surface  charges  are  included,  with  the  choice  of 
.V  =  5  X  1014  cm-2,  d  =  5  A.  G  (E )  is  plotted  on  a  log  scale 
due  to  its  large  variation  when  the  energy  is  varied  from  0  to 
1  eV.  For  example,  at  a  doping  of  1019  cm-3,  G[E)  ranges 
over  10  orders  of  magnitude.  Notice  that  G  (E )  has  a  discon¬ 
tinuity  in  its  slope  at  E  =  qVm.  This  is  due  to  Eq.  (2),  where 
the  transmission  below  the  barrier  is  given  by  the  WKB  ap¬ 
proximation,  while  the  transmission  above  the  barrier  was 
taken  to  be  1. 

For  adoping  of  1017  cm  ~ 3,  the  maximum  in  G  (E )  occurs  at 
q  ,  and  the  total  conductance  has  roughly  equal  contribu¬ 
tions  from  carriers  going  over  the  barrier  in  thermionic  emis¬ 
sion  and  Fowler-Nordheim  tunneling.16  Observe  that  the 
tunneling  contribution  is  significant  only  for  energies  down 
to  —0.2  eV  below  qVm,  since  the  tunneling  length  increases 
strongly  for  decreasing  energy.  There  is  no  direct  tunneling, 
since  £FM  =  Ep^  <0.  At  a  doping  of  10'“  cm-3  the  maxi¬ 
mum  in  G  (E )  increases  and  it  occurs  at  E  =  0.42  e  V,  which  is 
below  q  Vm .  Due  to  the  narrowing  of  the  barrier,  we  get  sig¬ 
nificant  conductance  over  most  of  the  tunneling  energies. 
Fowler-Nordheim  tunneling  dominates  the  conductance, 
while  some  direct  tunneling  occurs,  for  E  <Ers  —  0.06  eV. 
When  the  doping  is  raised  to  1019  cm-3,  G{E)  increases 
strongly,  with  its  maximum  occurring  at  E  =  0. 27  eV ■^Eri . 
Now  direct  tunneling  is  comparable  to  Fowler-Nordheim 
tunneling.  Finally,  for  a  doping  of  1020  cm-3,  Vm  falls  con- 


Fic.  3  Specific  contact  resistance  vs  doping,  for  4,  =0.6  eV  and  V  =  0  V, 
with  and  without  negative  surface  charge,  with  the  results  from  CFS 
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siderably  after  being  approximately  constant  at  lower  dop¬ 
ings.  This  is  due  to  the  negative  surface  states  “annulling" 
most  of  the  thin  barrier.  Here  EFS  =0.33  eV,  and  so,  direct 
tunneling  is  seen  to  dominate  the  conductance. 

The  specific  contact  resistance  as  a  function  of  doping  is 
presented  for  three  different  intrinsic  barrier  heights:  0.6, 
0.8,  and  1.0  eV  in  Figs.  3-5,  respectively.  Our  results  include 
all  of  the  effects  described  above.  For  comparison  we  have 
plotted  the  results  assuming  no  negative  surface  charge,  and 
also  those  from  CFS.  In  Fig.  4,  where  <pB0  =  0.8  eV,  we  also 
show  the  parabolic  result.  This  occurs  if  in  Eq.  ( 10]  we  set  ne 
=0,  regardless  of  whether  or  not  the  semiconductor  is  de¬ 
generately  doped. 

Comparison  of  the  results  in  these  figures  allows  us  to 
draw  a  number  of  conclusions.  First,  let  us  describe  those 
which  are  common  to  botn  CFS's  and  our  models.  The  con¬ 
tact  resistance  is  seen  to  increase  as  the  barrier  height  in¬ 
creases,  at  fixed  doping.  For  example,  at  a  doping  of  10''' 
cm-3,  in  going  from  <t>BQ  =  0.6  to  1.0  eV,  the  contact  resis¬ 
tance  increases  by  about  two  orders  of  magnitude  in  both 
models.  For  dopings  between  1018  and  10 19  cm-3,  the  con¬ 
tact  resistance  falls  steeply  by  about  seven  orders  of  magni¬ 
tude.  Recall  from  Fig.  2  and  the  discussion  above,  that  over 
this  doping  range  for  d>B0  =  0.8  eV,  tunneling  increased  to 
dominate  the  total  conductance,  and  G[E)  increased  strongly 
over  all  the  tunneling  energies.  Similar  results  are  seen  in  our 
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Fig  4  Specific  contact  resistance  vs  doping,  for  <j>n  -  0  8  eV  and  V  —  0  V, 
with  and  without  negative  surface  charge,  with  the  results  from  CFS  For 
comparison,  the  results  assuming  a  parabolic  conduction  hand  in  f  he  deple¬ 
tion  region  are  shown 


fj  t - - - , - - - - - - - T _ _ 
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log  DOPING  (ern-1) 

Fig.  5.  Specific  contact  resistance  vs  doping,  for  d>a  =  10  e V  and  V  ~  1 
with  and  without  negative  surface  charge,  with  the  results  from  CFS. 

model  for  </>BO  =  0.6  and  1.0  eV.  Thus  we  may  consider  the 
doping  range  from  1018  to  1019  cm-3  as  a  transition  from  a 
Schottky  barrier  to  an  Ohmic  contact.  As  a  final  point  of 
similarity,  note  that  CFS’s  results  and  our  nonparabolic 
curves  (with  no  surface  charge)  tend  to  merge  for  dopings 
near  10l8cm-3.  This  is  to  be  expected,  since  the  two  models 
differ  in  their  expressions  for  the  tunneling  transmission, 
and  hence  should  yield  the  same  results  when  the  tunneling 
contribution  is  small. 

Next,  consider  the  difference  between  the  results  of  the 
two  models.  The  nonparabolic  curves  tend  to  lie  below 
CFS’s  curves  for  dopings  greater  than  1018  cm-3.  Between 
dopings  of  1018  and  1019  cm  “ 3,  our  nonparabolic  curves  fall 
with  steeper  slope,  and  continue  doing  so  up  to  1021  cm  "  \ 
whereas  by  10:o  cm-3  the  slopes  on  CFS’s  curves  are  much 
reduced  Thus  the  gap  between  the  nonoarabolic  and  CFS's 
curves  increases  with  increasing  doping.  That  the  differences 
between  the  curves  increases  with  doping  is  expected,  since 
the  tunneling  contribution  dominates  the  conductance,  and 
the  change  in  the  function  relating  the  imaginary  part  of  the 
wavevector  to  the  energy  plays  an  important  role,  and  ac¬ 
counts  for  most  of  the  changes.  Thus,  at  high  dopings  (  >  1 0 1 9 
cm~3)  it  may  be  possible  to  obtain  metal-semiconductor 
junctions  with  substantially  lower  contact  resistance  than 
previously  thought  feasible. 

To  illustrate  the  relative  importance  of  the  various  factors 
in  our  theory,  consider  Fig.  4.  For  dopings  below  1018  cm  ' 
or  greater  than  1020  cm'3,  there  is  negligible  difference 
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POTENTIAL  IN  DEPLETION  REGION 


Fig  6,  Potential  in  the  depletion  region,  for  doping  =■  8  v  10|4‘  cm  <&„ 
=  0  8  eV,  and  V  =  0  V  The  r  direction  is  normal  to  the  interface  and  the  c 
direction  is  along  the  interface,  wuh  resolutions  ofO  5  and  10  A.  respective¬ 
ly  Also,  the  direction  is  in  the  <100)  direction  of  life  lattice 

between  the  parabolic  and  nonparabolic  results.  Indeed, 
between  these  doping  limits,  the  difference  is  less  than  0.2  on 
the  log  scale.  [That  the  parabolic  curve  falls  below  the  non¬ 
parabolic  curve  is  obvious  from  Eq.  (9).]  By  our  choice  of 
parameters  for  the  surface  charges,  greater  difference  is  seen 
between  this  curve  and  the  nonparabolic  curve  than  between 
that  and  the  parabolic  curve.  Qualitatively,  including  sur¬ 
face  charges  causes  little  difference,  except  at  dopings 
greater  than  102Ocm-3,  where  the  curve  flattens.  The  reason 
for  this  can  be  seen  from  Fig.  2.  At  102u  cm  the  negative 
surface  charges  have  greatly  reduced  the  barrier.  Effectively, 
there  is  no  barrier,  so  further  increases  in  doping  will  have 
little  effect.  Quantitatively,  we  may  say  that  the  presence  of 
any  negative  surface  charge  will  accentuate  the  difference 
between  our  results  and  CFS’s  results.  The  effect  of  neglect¬ 
ing  image-force  lowering  is  slight;  it  causes  the  curves  to  shift 
upwards  slightly  by  —0.2. 

IV.  FLUCTUATIONS 

One  significant  effect  that  has  not  been  included  in  this 
theory  is  fluctuations  in  the  potential  due  to  the  discreteness 
in  the  doping.  A  simple  Monte  Carlo  model  has  been  used  to 
calculate  the  fluctuations,  wherein  dopants  have  been  ran¬ 
domly  placed  as  substitutional  impurities  within  a  GaAs  lat¬ 
tice,  in  the  depletion  region.  One  such  result  is  illustrated  in 


Fig.  6,  for  a  doping  of  8x  1019  cm-3  and  <t>B0  =  0.8  eV.  The 
potential  is  shown  in  a  45  A  x  1 95  A  rectangle  in  an  xz  plane, 
where  x  is  normal  to  the  interface  and  z  is  parallel  to  the 
interface.  The  potential  is  seen  to  exhibit  significant  varia¬ 
tions.  In  a  future  paper,17  we  shall  describe  our  model,  and 
calculate  the  effects  of  such  fluctuations  on  the  conductance 
as  a  function  of  doping,  and  compare  these  to  the  results 
presented  here. 


V.  CONCLUSION 

We  have  presented  a  model  for  the  calculation  of  contact 
resistance  and  conductance  as  a  function  of  energy  of  the 
tunneling  electron,  for  a  metal-n-GaAs  junction,  where  the 
GaAs  is  heavily  doped,  so  that  tunneling  is  significant.  This 
model  uses  the  WKB  approximation  and  the  two-band  mod¬ 
el  to  obtain  a  nonanalytic  expression  for  the  transmission 
probability.  The  contact  resitance  at  zero  bias  is  compared 
with  that  from  CFS’s  model.  We  suggest  that  our  model 
determines  the  transmission  probability  more  accurately 
than  CFS’s  model,  which  assumes  constant  transmission  for 
all  tunneling  energies. 
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Resistance  fluctuations  in  ohmic  contacts  due  to  discreteness  of  dopants 

W.  J.  Boudville  and  T.  C.  McGill 

California  Institute  of  Technology.  Pasadena.  California 't l Iff 
(  Received  30  September  1985,  accepted  for  publication  28  January  1986) 

The  role  of  fluctuations  in  the  potential  due  to  randomly  distributed  dopants  in  the  depletion  layer 
of  a  metal-semiconductor  junction  is  explored  To  be  specific,  the  case  of  n-GaAs  is  considered. 

Monte  Carlo  simulation  techniques  are  used  to  calculate  the  potential  in  thejunction  By  using  the 
WKB  approximation  and  the  two-band  model,  the  small-signal  resistance  at  zero  bias  is  found  to 
be  lowered  by  up  to  half  an  order  of  magnitude  from  the  result  assuming  a  continuum  distribution 
of  dopant  charge.  The  resistance  is  found  to  vary  by  an  order  of  magnitude  in  the  plane  of  the 
interface  resulting  in  a  very  nonuniform  distribution  of  the  current. 


One  of  the  standard  models  of  ohmic  contacts  is  a  metal- 
semiconductor  junction  in  which  the  semiconductor  is  suffi¬ 
ciently  heavily  doped  that  tunneling  is  the  major  transport 
mechanism  through  theSchottky  barrier.1  This  ideal  model 
of  the  ohmic  contact  has  been  realized  in  recent  years.  Oh¬ 
mic  contacts  to  n- GaAs  have  been  made  by  molecular  beam 
epitaxy,  where  the  layers  of  GaAs  have  been  heavily  doped 
with  Sn  (Ref.  2)  or  Si,’  with  dopings  in  the  range  of  lO1'*- 
lO^cm"3.  At  these  dopings,  the  current  through  a  metal-n- 
GaAs  junction  thus  formed  is  dominated  by  tunneling.  The 
model  of  ohmic  contacts  by  Chang,  Fang,  and  Sze  (CFS)1 
has  been  most  commonly  used  to  analyze  results  such  as 
those  in  Refs.  2  and  3. 

In  a  previous  paper4  we  have  presented  an  improved 
theory.  Here  we  extend  that  theory  to  include  the  effects  of 
fluctuations  in  the  potential  due  to  the  randomness  of  the 
doping  within  the  depletion  region.  These  fluctuations  could 
be  important  in  the  range  of  doping  where  ohmic  contacts 
are  formed.  For  example,  at  dopings  around  10'9  cm  ’, 
where  the  average  interdopant  separation  is  —  30  A,  the  de¬ 
pletion  length  is  —  1 20  A.  Thus  an  electron  traversing 
through  the  depletion  region  would  be  likely  to  encounter, 
not  a  potential  due  to  a  continuum  of  ionized  dopants,  but  a 
potential  influenced  strongly  by  a  few  dopants,  and  hence 
varying  considerably. 

The  WKB  approximation  is  used  to  give  the  transmis¬ 
sion, 

expf  —  2/  f  if  0 <,E<*qV„ 

T(E)=  \  J*.  J  ,  (1) 

„1  if  E>qVm 

where  Vm  is  the  height  of  the  potential  barrier,  with  respect 
to  the  conduction  band  outside  the  depletion  region,  and  x,, 
x2  are  the  turning  points  of  the  forbidden  region  for  k,  where 
the  x  direction  is  normal  to  the  interface.  Within  the  forbid¬ 
den  region,  from  the  two-band  model,5 


u„y'\ 


where  is  the  band  gap  energy  and  U(x)  is  the  potential 
energy  in  the  depletion  region  At  zero  bias,  the  small-signal 
resistance  is  given  by 

—  =  G,  =  1  dE - - . 

R  2rrft'  J„  exp{  (E  -  E ,w  )/kB  T  j  +■  1 


Previously,4  U  was  found  by  including  nonparabolic  correc¬ 
tions  and  the  L  point  states.  The  depletion  region  was  as¬ 
sumed  to  consist  of  a  uniform  continuum  of  ionized  dopants. 

To  include  the  effects  of  fluctuations  in  the  potential,  we 
have  taken  a  different  approach.  In  principle,  one  could  put 
down  the  dopants  in  a  random  fashion,  solve  for  the  electro¬ 
static  potential,  and  then  solve  for  the  transport  through  the 
potential  barrier.  However,  this  problem  is  much  too  diffi¬ 
cult  and  we  were  forced  to  make  some  approximations. 

Consider  a  rectangle  with  one  side  in  the  plane  of  the 
interface  and  the  other  side  normal  to  the  interface,  and  of 
length  xd.  For  the  purpose  of  calculation,  the  length  of  the 
rectangle  in  the  plane  of  the  interface  is  taken  to  be  200  A.  To 
compute  the  potential  in  this  rectangle,  we  have  made  a 
number  of  simplifying  approximations.  First,  we  assume 
that  the  width  of  the  depletion  region  does  not  vary  in  the 
plane  of  the  interface.  We  take  the  depletion  width  to  be 
given  by  the  standard  continuum  approximation, 

fie  V/J 

Xd  —  ^ — ~  <*.0  "F  Efs)J  >  (4) 

where  <f>b0  is  the  intrinsic  barrier  height,  n  is  the  density  of 
donors,  and  En  is  the  fermi  level  in  the  semiconductor. 
is  found  from  Eqs.  (6)  and  (7)  in  Ref.  4. 

Second,  we  divide  the  contnbutions  to  the  potential  seen 
by  the  electron  into  three  parts, 

U  =  -^-+Ux  +  U2.  (5) 

16 irex 

The  first  term  on  the  right-hand  side  is  the  image  term  of  the 
electron.  The  term  (/,  is  due  to  dopants  far  from  the  electron, 
and,  hence,  should  cause  small  fluctuations  in  the  potential. 
The  term  U2  is  the  potential  due  to  ionized  dopants  near  the 
electron,  and  is  responsible  for  most  of  the  fluctuations  in  the 
potential. 

To  be  specific,  the  regions  contributing  to  the  potential 
seen  by  an  electron  are  illustrated  in  Fig.  1.  We  divide  the 
depletion  region  into  two  parts:  a  cylinder  of  radius  R  and 
depth  xd,  and  the  region  outside  the  cylinder.  The  radius  of 
the  cylinder  is  taken  to  be  three  times  the  thickness  of  the 
depletion  layer.  [/?  =  lxd  in  Eq.  (6)  below.]  The  contribu¬ 
tion  to  the  potential  from  the  ionized  donors  outside  the 
cylinder  is  calculated  by  assuming  that  the  charge  density  is 
continuous.  The  image  in  the  metal  is  included;  however,  the 
image  potential  due  to  the  undepleted  GaAs  is  not  included 
because  of  the  small  electron  density  when  compared  to  a 
metal. This  yields 
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FIG.  1  Contnbution  to  the  potential  at  (x,  y.  2)  in  the  depletion  region, 
from  a  differential  charge  at  (jr„,  y,„  z„)  outside  the  cylinder,  is  shown  Also 
illustrated  is  the  contribution  from  the  differential  image  at  (  -  >>,,,  z„) 


FIG  2.  Random  potential  for  one  configuration  of  ionized  donors  in  a  rec¬ 
tangle  in  the  depletion  region.  The  average  doping  is  2.5  X  10'“  cm  \  The 
Schottky  barrier  height  is0.8eV  The  applied  bias  is  zero.  The  a  direction  is 
normal  to  the  interface,  and  the  z  direction  is  along  the  interface  The  x 
direction  is  the  ( 100)  direction  of  the  GaAs  lattice. 
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where  <t>^,  —  U{ 0)  in  the  absence  of  image  force  lowering. 

The  charged  donors  inside  the  cylinder  are  taken  to  be 
discrete,  and  their  contribution  to  the  potential  is  calculated 
in  detail  including  their  image  in  the  metal.  (The  image  con¬ 
tribution  from  the  undepleted  GaAs  is  neglected  again.)  Of 
course,  this  contribution  to  the  potential  will  depend  in  de¬ 
tail  on  the  spatial  arrangement  of  the  ionized  donors  found 
in  the  depletion  region.  The  spatial  configuration  of  the  ion¬ 
ized  donors  is  given  by  using  a  random  number  generator  to 
distribute  the  ionized  dopants  on  one  of  the  two  atomic  sites 
in  a  zincblende  lattice.  The  density  is  set  by  the  doping.  The 
atoms  are  inside  a  rectangular  slab  of  GaAs  whose  thickness 
equals  xd  The  dimensions  of  the  rectangular  slab  are  chosen 
to  be  large  enough  to  include  the  rectangle  within  which  we 
are  finding  U  and  cylinders  centered  on  any  point  within  the 
rectangle. 

The  contnbution  to  conductance  per  unit  area  for  a  giv¬ 
en  cross-sectional  area  is  calculated  using  Eqs.  ( l )-( 3).  The 
transport  is  calculated  as  if  the  potential  depended  only  on 
the  x  direction.  The  formulas  are  evaluated  on  an  equally 
spaced  grid  with  a  gnd  spacing  of  5  A 

To  examine  the  variations  in  potential  and  transport 
that  would  be  obtained  for  other  configurations,  we  have 
used  the  random  number  generator  to  generate  a  number  of 
different  spatial  arrangements  of  atoms  The  entire  transport 
calculation  was  earned  out  for  each  case 

Figure  2  shows  an  example  of  a  potential  found  by  this 
method  at  a  doping  of  2  5  •  1 0 ' ^  cm  1  and  <3,.,,  0  8  eV 

The  contour  lines  are  lines  of  constant  z,  spaced  5  A  apan. 


where  z  is  in  the  plane  of  the  interface.  It  is  seen  that  the 
potential  fluctuates  significantly.  Thus,  at  a  given  tunneling 
energy,  the  length  of  the  forbidden  region  can  vary  widely, 
especially  if  a  dopant  is  nearby,  causing  the  potential  to  fall. 
From  Eqs.  ( 1 ) ,  ( 2 ) ,  and  ( 3 ) ,  it  is  seen  that  this  variation  in 
length  can  cause  a  large  variation  in  Rc . 

In  Figs.  3  and  4  we  show  the  results  of  our  calculations 
for  bamer  energies  of  0.8  and  1.0  eV  along  with  the  results 
for  various  models.  The  a  direction  has  been  chosen  to  be  the 
< 100)  direction  of  the  lattice.  The  results  are  presented  as 
specific  contact  resistance  (the  reciprocal  of  the  conduc¬ 
tance  per  unit  area)  as  a  function  of  doping.  Each  dot  in  the 
figures  is  obtained  from  a  different  random  configuration  of 
dopants  and  represents  Gc ,  averaged  over  200  A  in  the  plane 
of  the  interface,  where  individual  Gt ’s  are  found  at  40  :  co¬ 
ordinates  spaced  5  A  apart.  Between  four  and  six  simula¬ 
tions  are  done  at  each  doping  value.  Having  obtained  an 
averaged  Gc  for  each  dopant  configuration,  these  averaged 
results  are  in  turn  averaged  at  each  doping  value.  (These  are 
not  displayed  in  Figs.  3  and  4. )  A  least-squares  polynomial  is 
then  fitted  to  the  resultant  [  (log)0  n,  log,„  rc )  ]  at  each  bar¬ 
rier  height  The  curves  are  shown  in  the  figures  under  the 
label  “Monte  Carlo  (100)."  For  comparison,  the  results 
f.om  CFS1  and  Ref.  4  (labeled  “nonparabolic" )  are  also 
displayed.  The  results  from  Ref.  4  are  obtained  by  assuming 
a  continuum  distribution  of  charged  dopants  within  the  de¬ 
pletion  region,  rather  than  the  discrete  case  considered  here. 
The  results  in  Figs.  3  and  4  took  —  30  days  CPU  time  on  a 
VAX  1 1/785 
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FIG  3  Specific  contact  resistance  vs  doping.  The  value  of  the  Schottky 
bamer  is  0  8  eV  The  calculations  are  carried  out  at  zero  bias.  The  points  are 
the  results  of  the  Monte  Carlo  simulations.  The  dashed  line  is  produced 
from  the  Monte  Carlo  results  by  the  procedure  described  in  the  text.  The 
solid  line  is  produced  using  the  theory  described  in  Ref.  4.  The  dash-dotted 
line  is  the  result  obtained  from  the  original  theory  by  Chang,  Fang,  and  Sze 
( see  Ref  I ) 

A  number  of  important  points  can  be  obtained  by  study¬ 
ing  the  results  presented  in  these  figures.  From  the  figures,  it 
is  seen  that  for  dopings  <  10'“  cm"1,  the  Monte  Carlo  re¬ 
sults  exhibit  small  fluctuations  about  the  continuum  curve. 
However,  as  the  doping  is  increased,  the  average  resistance 
falls  about  half  an  order  of  magnitude  below  the  continuum 
curve.  Furthermore,  the  fluctuations  increase  to  as  much  as 
one  and  a  half  orders  of  magnitude  between  results  from 
different  dopant  configurations  at  the  same  doping  and  bar¬ 
rier  height.  Thus  current  through  the  junction  may  flow 
preferentially  through  “hot  spots” — regions  of  high  conduc¬ 
tance 

In  summary,  we  have  found  that  the  discreteness  of  the 
ionized  dopants,  within  the  depletion  region  of  a  metal-n- 
GaAs  junction,  causes  a  significant  lowering  in  the  small- 
signal  resistance  at  zero  bias,  as  a  function  of  doping.  The 
small-signal  resistance  can  be  up  to  half  an  order  of  magni¬ 
tude  less  than  those  resulting  from  assuming  a  continuous 
charge  distribution,  and  up  to  an  order  and  a  half  less  than 
the  results  obtained  by  CFS  '  More  importantly,  large  fluc¬ 
tuations  in  the  conductance  are  found.  These  fluctuations 
may  have  serious  implications  for  the  transport  properties  of 
metal-semiconductor  junctions  The  results  indicate  that  for 
a  semiconductor  which  is  uniformly  doped  on  a  macroscopic 
scale,  the  conductance,  and  hence  the  current,  can  differ  by 


FIG.  4  Specific  contact  resistance  vs  doping  The  value  of  the  Schottky 
barrier  is  1 .0  eV.  The  calculations  are  earned  out  at  zero  bias.  The  points  are 
the  results  of  the  Monte  Carlo  simulations  The  dashed  line  is  produced 
from  the  Monte  Carlo  results  by  the  procedure  desenbed  in  the  text  The 
solid  line  is  produced  using  the  theory  desenbed  in  Ref  4  The  dash-dotted 
line  is  the  result  obtained  from  the  original  theory  by  Chang,  Fang,  and  Sze 
( see  Ref  1 ) . 

over  an  order  of  magnitude  in  different  spatial  regions  of  the 
interface.  This  large  spatial  variation  in  the  current  may  be  a 
significant  contribution  to  the  noise  of  a  metal-semiconduc¬ 
tor  junction,  and  hence  may  adversely  affect  the  perfor¬ 
mance  of  micron  or  submicron  size  devices  employing  such 
junctions. 

We  expect  similar  phenomena  to  occur  in  other  doped 
barrier  structures.  Studies  are  under  way  to  explore  the  role 
of  potential  fluctuations  in  other  device  structures  with 
small  characteristic  dimensions. 
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The  role  of  defects  in  heterojunctions  was  investigated.  The  density  of  such  defects  required  to  pin 
the  Fermi  level  or  to  affect  the  band  offset  was  estimated  using  simple  electrostatic  considerations. 

We  conclude  that  it  is  very  unlikely  that  defects  play  any  role  in  determining  the  band  offsets,  but 
they  might  affect  the  Fermi-level  position  at  the  interface. 

PACS  numbers:  73.40.Lq,  71.25.Tn,  61.70.  -  r 


I.  INTRODUCTION 

Experiments  performed  on  clean,  cleaved  surfaces  of  III-V 
semiconductors  have  demonstrated  that  the  Fermi  level  at 
the  surface  is  pinned  by  deposition  of  very  thin  coverages 
(submonolayer  to  several  atomic  layers)  of  either  metal,'-9 
semiconductor,10-12  or  oxygen1,2,4,8,13-17  atoms.  The  fact 
that  the  Fermi-level  position  at  the  surface  was  almost  the 
same  for  very  different  adatoms,  as  well  as  the  fact  that  clean, 
cleaved  III-V  semiconductors  do  not  have  surface  states  in 
the  gap,2,18-21  led  Spicer  et  a/.4,22-24  to  propose  the  unified 
defect  model.  According  to  that  model,  chargeable  native 
defects  are  formed  upon  deposition  of  minute  quantities  of 
adatoms  (being  either  metal,  semiconductor  or  oxygen 
atoms),  and  these  defects  may  determine  the  position  of  the 
Fermi  level  at  III-V  semiconductor  interfaces. 

The  unified  defect  model  of  Spicer  et  al.,  which  explains 
the  position  of  the  Fermi  level  at  semiconductor  interfaces  in 
the  limit  of  very  thin  coverages  and  is  successful  in  predict¬ 
ing  some  interfacial  properties,  is  currently  used  in  the  litera¬ 
ture  to  predict  the  behavior  of  semiconductor  interfaces 
when  the  coverage,  either  metallic  of  semiconducting,  is 
thick.  For  example,  it  was  predicted  that  in  a  metal-semi¬ 
conductor  interface,  the  Fermi  level  would  be  pinned  at  the 
various  defect  levels, 3,9-25  jnd  that  this  pinning  position 
would  depend  on  the  semiconductor  type26  but  not  on  the 
metal  work  function.  Other  researchers  have  applied  the  un¬ 
ified  defect  model  to  semiconductor  heterojunctions  and 
have  predicted  that  the  defect  levels  might  pin  the  Fermi 
le»el  at  the  interface.26 

In  this  paper  we  examine  the  role  defects  might  play  in 
determining  band  offsets  and  the  position  of  the  Fermi  level 
at  semiconductor  heterojunctions.  In  a  recent  paper,27  we 
have  argued  that  the  Fermi-level-pinning  mechanism  for 
very  thin  metallic  coverages  of  semiconductor  surfaces  is 
different  from  the  mechanism  of  Fermi-level  pinning  at  the 
interface  between  two  bulk  materials,  the  main  difference 
between  the  two  cases  being  the  source  of  charge  that  bal¬ 
ances  the  charge  captured  on  the  surface  defects.  In  the  case 
of  very  thin  coverage,  this  source  must  be  the  depletion 
charge  in  the  bulk  semiconductor.  Therefore,  only  surface 
acceptors  can  be  charged  on  bulk  n-type  semiconductors, 
and  only  surface  donors  can  be  charged  on  bulk  p-type  semi¬ 
conductors,  even  if  both  surface  acceptors  and  surface  do¬ 
nors  are  present  in  both  cases.  This  will  result  in  pinning  of 


the  Fermi  level  at  the  surface-acceptor  level  on  n-type  and  at 
the  surface-donor  level  on  p-type  semiconductors.  The  den¬ 
sity  of  defects  required  to  pin  the  Fermi  level  is  comparable 
to  the  total  depletion  charge,  which  is  of  the  order  of  10’ 2 
cm  - 2  for  a  doping  concentration  of  1 0 1 7  cm  ~ 3  in  most  wide¬ 
band-gap  semiconductors.  Such  a  low  defect  density  can  be 
formed  by  a  submonolayer  coverage  of  the  semiconductor.  It 
was  also  argued  that  if  the  defect  density  will  be  further  in¬ 
creased,  the  two  “pinned"  Fermi-level  positions  would  move 
and  eventually  merge.28  The  “pinning”  at  submonolayer 
coverage  is,  therefore,  not  complete. 

In  contrast  to  the  case  of  very  thin  coverages,  a  defect 
concentration  of  1012  cm-2  was  shown28  to  be  approximate¬ 
ly  two  orders  of  magnitude  too  low  to  pin  the  Fermi  level  at  a 
macroscopic  metal-semiconductor  junction.  In  this  paper 
we  shall  investigate  the  effect  of  defects  on  semiconductor 
heterojunctions  for  both  thin  and  thick  coverages  and  con¬ 
trast  it  with  our  previous  papers  concerning  a  metal-semi¬ 
conductor  interface.27,28  We  shall  show  that  the  effect  of  de¬ 
fects  in  a  bulk  heterojunction  is,  qualitatively,  very  similar  to 
the  case  of  a  submonolayer  coverage.  (This  was  implictly 
assumed  by  Kroemer29  in  his  treatment  of  interface  charges 
in  heterojunctions.)  Quantitatively,  we  shall  show  that  a  de¬ 
fect  concentration  of  1012  cm-2  is  also  too  low  to  pin  the 
Fermi-level  at  the  interface  between  two  semiconductors, 
although  it  does  have  an  effect  on  the  Fermi-level  position, 
which  is  more  pronounced  at  semiconductor  heterojunc¬ 
tions  than  at  metal-semiconductor  interfaces.  A  defect  den¬ 
sity  of  at  least  an  order  of  magnitude  more  will  be  shown  to 
be  required  to  pin  the  Fermi  level  in  a  heterojunction. 

We  would  like  at  this  point  to  clarify  the  meaning  of  “Fer¬ 
mi-level  pinning”.  This  term  means  that  the  position  of  the 
Fermi  level  at  the  interface,  relative  to  the  band  edges,  is 
independent  of  some  external  conditions.  The  degree  of  this 
independence  however,  is  important  to  our  discussion.  The 
experiments  leading  to  the  unified  defect  model  show  that 
the  Fermi  level  at  the  surfaces  of  III-IV  semiconductors  is 
"pinned”  by  submonolayer  coverages.  The  word  pinned  is 
used  in  this  context  to  mean  that  the  curve  of  Fermi-level 
position  versus  coverage,  flattens  out  and  reaches  some  “fi¬ 
nal”  position  at  submonolayer  coverages.  There  are,  how¬ 
ever,  differences  of  about  0. 1  eV  between  these  “final”  Fer¬ 
mi-level  positions  as  measured  on  surfaces  covered  with 
different  adatoms,  as  well  as  differences  of  about  0.2  eV 
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between  the  final  Fermi-level  positions  measured  on  n-  and 
p-type  semiconductors  covered  with  the  same  adatom.'6 
Moreover,  more  recent  experiments  using  a  thicker  coverage 
of  metallic  atoms30'3 1  demonstrate  that  this  final  Fermi-level' 
position  is  not  final  and  can  still  move  about  when  the  cover¬ 
age  becomes  thicker.  It  appears,  therefore,  that  enough  de¬ 
fects  are  created  at  the  surface  to  affect  significantly  the  Fer¬ 
mi  level  but  not  quite  to  pin  it.  In  our  paper  we  shall  reserve 
the  term  “pinning”  to  the  case  in  which  the  Fermi-level  posi¬ 
tion  at  the  interface  is  determined  by  the  defects  but  not  by 
the  bulk  doping. 

In  our  calculations,  we  assume  that  interface  defects  have 
no  effect  on  the  band  offsets.  To  justify  this  assumption,  we 
note  that  the  interface  defers  would  have  to  form  a  dipole 
layer  in  order  to  affect  the  band  offsets;  a  single  layer  of 
trapped  charge  cannot  abruptly  change  the  potential,  only 
its  slope.  Therefore,  three  conditions  have  to  be  satisfied  if 
the  defects  form  a  dipole  layer  which  affects  the  band  offsets. 
First,  both  donor-  and  acceptor-type  defects  should  be  pres¬ 
ent  at  the  interface.  Second,  the  layer  of  donors  should  be 
separated  from  the  layer  of  acceptors  by  some  distance.  This 
requirement  is  different  from  the  case  of  a  metal-semicon¬ 
ductor  interface,  in  which  the  metal  is  guaranteed  to  supply 
charges  of  the  opposite  sign  by  screening.  A  reasonable  value 
for  this  separation  is  several  angstroms.  It  cannot  be  smaller 
than  one  atomic  layer  separation,  and,  on  the  other  hand,  it 
cannot  be  larger  than  a  few  tens  of  angstroms,  since  the  band 
offset  is  abrupt.  The  third  condition  that  has  to  be  satisfied  is 
that  the  defect  density  would  be  high  enough.  We  can  esti¬ 
mate  the  charge  density  required  to  affect  the  band  offset  by 
0. 1  V  using  the  parallel-plate-capacitor  formula: 

ee° 

For  a  separation  of  d  =  5  A  between  the  defect  donors  and 
acceptors,  ar  d  assuming  a  dielectric  constant  on  the  order  of 
10,  one  obtains  a  defect  density  of  approximately  a  =  10’3 
cm-2.  This  defect  density  seems  to  be  very  high,  especially 
when  one  considers  the  fact  that  this  density  of  a  =  10'3 
cm  " 2  counts  only  those  defects  which  participate  in  this  nar¬ 
row  dipole  layer,  that  is,  only  those  defects  which  are  sepa¬ 
rated  by  several  angstroms  from  a  layer  of  defects  charged 
with  the  opposite  sign.  We  must  conclude,  therefore,  that  it 
is  very  unlikely  that  interface  defects  determine  the  band 
offsets  or  have  a  significant  influence  on  them.  Some  experi¬ 
mental  evidence  supports  this  assumption.32  33 

This  paper  is  organized  in  the  following  way.  In  Sec.  II,  we 
discuss  the  theoretical  methods.  In  Sec.  Ill,  we  present  the 
results  based  on  our  calculations.  The  conclusions  of  this 
paper  are  given  in  Sec.  IV. 


II.  THEORETICAL  METHODS 

The  schematic  of  our  model  is  given  in  Fig  1 .  Our  model 
consists  of  two  bulk  semiconductors  (infinitely  thick)  extend¬ 
ing  from  x  =  0  to  <  ither  side.  The  interface  at  x  =  0  is  as¬ 
sumed  to  be  abrupt.  The  band  discontinuities  at  the  interface 
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Fic.  1  The  geometry  *nd  energy  levels  of  our  heterojuncuon  model.  We 
assume  a  very  thin  defect  layer  between  the  two  bulk  semiconductors,  con¬ 
taining  both  defect  donors  and  defect  acceptors.  The  ionization  energies  of 
these  defects  relative  to  the  band  edges  of  the  two  semiconductors  are  de¬ 
noted  by  “£>"  and  "A  ",  respectively  We  assume  that  the  band  offset  ide- 
noted  by  |  is  independent  of  the  defects. 


are  prescribed,  and  the  interface  contains  a  layer  of  defects 
(both  donors  and  acceptors)  with  prescribed  densities  and 
ionization  energies. 

The  calculation  of  the  Fermi  level  at  the  interface  is  based 
on  conservation  of  charge.  It  uses  a  property  of  semi-infinite 
semiconductors,  namely,  that  the  Fermi-level  position  at  the 
surface  (x  =  0)  determines  the  total  charge  per  unit  area  in 
the  semi-infinite  semiconductor.  Thus,  if  tj{x)  denotes  the 
position  of  the  Fermi  level  relative  to  the  conduction  band 
edge, 

V  (x)=fi-Ec(x),  (2.1) 

and  Q[x)  denotes  the  total  number  of  charges  per  unit  area 
between  x  and  <®, 


Q\x) 


-J>- 


+  N /  -N;  ]dx, 


(2.2) 


then  Poisson’s  equation  relates  t/(x|  to  Q(x ): 

-  —  ew. 

ax  ee0 


dQ 

dx 


=  n{7j)  -  p{rj)  +  N  a  \rj)  —  /V  /  [rj). 


(2.3) 


Here,  q  is  the  electron  charge  (taken  positive),  e  is  the  dielec¬ 
tric  constant  of  the  semiconductor,  n  and  p  are  the  densities 
of  the  electrons  in  the  conduction  band  and  holes  in  the  va¬ 
lence  band,  and  N  J  and  N~  are  the  densities  of  ionized 
donors  and  acceptors,  respectively.  Expressions  for  n,  p, 
N  /  ,  and  N  ~  can  be  found  in  textbooks  on  the  physics  of 
semiconductors34  and  are  not  repeated  here. 

Dividing  the  two  equations  in  Eq.  (2.3)  by  each  other,  inte¬ 
grating  from  rf{  oo )  to  rj  and  taking  the  square  root,  we  obtain 
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Q(v)  =  ± 


[n{v‘)  -Piv')  -  *v;  l»7')  +  y  a  .[v')]drj' . 


(2.4) 


Equation  (2.4)  relates  the  Fermi-level  position  at  the  inter¬ 
face  ( 7 j)  to  the  total  number  of  charges  per  unit  area  Q  (77).  The 
sign  in  front  of  the  square  root  in  Eq.  (2.4)  is  positive  if 
77(0)  <  rj{  00 ),  negative  otherwise.  tj{  00 ),  the  bulk  Fermi-level 
position,  is  calculated  from  the  local-charge-neutrality  con¬ 
dition: 


n[v)  -  p\v)  +  N a  (77)  -  IV/ (77)  =  0.  (2.5) 


Now,  we  can  put  the  interface  together.  Given  77,(0),  the 
Fermi-level  position  at  the  interface,  relative  to  the  conduc¬ 
tion-band  edge  of  semiconductor  1,  we  obtain  772(0),  the  Fer¬ 
mi-level  position  at  the  interface  relative  to  the  conduction- 
band  edge  of  semiconductor  2,  from  the  band  discontinuity 
(which  is  an  input  to  this  model) 

772(0)  -  77,(0)  =  -AEC.  (2.6) 

From  77,(0)  and  t72(0),  we  obtain  the  total  number  of  charges 
per  unit  area  in  the  two  bulk  semiconductors,  Q{  and  Q2 
using  Eq.  (2.4).  Finally,  the  Fermi-level  position  at  the  inter¬ 
face  determines  the  fraction  of  ionized  interface  defects.  This 
fraction  multiplied  by  the  density  of  these  defects  gives  the 
number  of  ionized  defects  per  unit  area  Q  ~  and  Q  * . 


q :  = 


1  +  gd  exp 


(^) 


= 


1  +  ga  expi 


(W 


(2.7) 


Here,  ad  and  oa  are  the  densities,  gd  and  ga  are  the  degener¬ 
acy  factors,  and  E ^  and  £M  are  the  ionization  energies  of 
interface  donors  and  acceptors,  respectively. 

The  Fermi-level  position  at  the  interface  is  now  deter¬ 
mined  by  an  overall  charge-neutrality  condition: 

Q,  +Q2  +  Q*  -Qz  =0.  (2.8) 

In  the  case  shown  in  Fig.  1,  both  semiconductors  are  n  type. 
Therefore,  the  depletion  regions  on  both  sides  of  the  inter¬ 
face  (Q 1  and  Q2)  are  positively  charged.  The  Fermi  level  (74)  at 
the  interface  approximately  equals  the  defect  acceptor  level 
(“A  ”),  and,  therefore,  approximately  half  the  defect  accep¬ 
tors  will  be  negatively  charged  and  balance  the  depletion 
charge.  Since  the  Fermi  level  at  the  interface  is  above  the 
defect-donor  level  (“D  "),  most  of  the  defect  donors  will  re¬ 
main  neutral. 


III.  RESULTS 


A.  Submonolayer  coverage 

The  case  of  submonolayer  coverage  was  treated  be¬ 
fore.2627  The  model  for  treating  this  case  assumes  that  the 
adatoms  cannot  contribute  any  significant  amount  of 
charge,  and  their  only  role  is  the  creation  of  defects.  There¬ 
fore,  one  might  assume  Q2  =  0  in  Eq.  (2.8),  and  that  the 
charge  due  to  the  surface  states  balances  the  substrate  deple¬ 
tion  charge.  Figure  2  shows  the  resulting  Fermi-level  posi¬ 
tion  at  the  interface  as  a  function  of  surface-defect  density  for 
both  /I-  and  77-type  semiconductors.  We  assume  that  both 


defect  donors  and  defect  acceptors  are  present  at  the  inter¬ 
face  at  the, same  concentration.  The  shape  of  both  curves  is  in 
good  qualitative  agreement  with  experimental  results  (com¬ 
pare  to  Figs.  3  and  8  in  Mbnch26). 

We  should  observe  here  that  our  assumption  about  the 
relative  abundance  of  defect  acceptors  and  defect  donors  is 
not  restrictive  in  any  way.  If  we  assume  any  other  fixed  ratio 
of  defect  acceptors  to  defect  donors,  there  will  be  a  corre¬ 
sponding  change  in  the  horizontal  scale  of  the  two  curves  in 
Fig.  2.  By  comparing  the  curves  of  Fig.  2  with  experimental 
results  (for  example,  the  point  of  maximum  curvature),  we 
may  conclude  that  both  defect  donors  and  defect  acceptors 
are  present,  and  their  concentration  is  of  the  same  order  of 
magnitude  though  not  necessarily  equal. 

The  curves  in  Fig.  2  seem  to  stabilize  at  a  defect  density  on 
the  order  of  101J  cm-2.  (In  fact,  these  two  curves  have  some 
small  but  nonzero  slope  at  that  density,  and  they  slowly  ap¬ 
proach  each  other).  This  defect  density  corresponds  to  a  dop¬ 
ing  concentration  of  1017  cm-3  and  scales  like  the  square 
root  of  the  doping.  Moreover,  the  pinning  positions  for  n- 
and  /7-type  semiconductors  are  substantially  different  and 
determined  essentially  by  the  defect-donor  level  for  a  /7-type 
semiconductor  and  by  the  defect-acceptor  level  for  an  n-type 
semiconductor.  One  should  be  aware  of  the  fact  that  this 
behavior  is  typical  for  the  case  of  very  thin  coverage.  A  thick 
coverage  of  either  a  metal  or  a  semiconductor  can  contribute 
a  significant  amount  of  charge  to  the  system.  In  other  words, 
Q2  in  Eq.  (2.8)  cannot  be  neglected,  and  as  a  result,  the  behav¬ 
ior  of  the  Fermi-level  position  at  the  interface  might  be  dif¬ 
ferent.  This  difference  in  behavior  between  thin  and  thick 


Fig  2.  The  position  of  the  Fermi  level  at  the  surface  as  a  function  of  defect 
concentration  for  GaAs.  The  positions  of  the  defect  levels  (either  acceptors 
or  donors)  are  shown  by  arrows.  The  density  axis  refers  to  either  the  defect 
donors  or  the  defect  acceptors. 
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coverage  was  treated  in  another  paper  for  the  case  of  a  met¬ 
al-semiconductor  interface.28  In  the  case  of  a  metal-semi-, 
conductor  interface,  the  charge  contributed  by  the  metal  can 
be  substantially  higher  than  that  contributed  by  the  semi¬ 
conductor  (Q2>Q i)-  As  a  result,  the  behavior  of  thick  cover¬ 
ages  was  shown  to  be  qualitatively  different  than  that  of  a 
thin  coverage.  We  now  treat  the  case  of  a  thick  semiconduc¬ 
tor  coverage,  i.e.,  a  bulk  heterojunction.  Here,  the  charge 
contributed  by  both  bulk  materials  is  of  the  same  order  of 
magnitude.  Therefore,  one  expects  a  similar  behavior 
between  thin  and  thick  coverages,  with  some  quantitative 
changes  due  to  the  presence  of  another  bulk  material. 

B.  Thick  coverage 

In  our  calculations,  we  have  used  the  Ge/GaAs  junction 
as  a  typical  heterostructure.  The  band  offset  was  chosen  to 
be35-36 

AEC  =  0.54  eV. 

This  band  offset  is  not  the  only  one  appearing  in  the  litera¬ 
ture.  Other  possible  values  for  AEC  are  0.411'2  and  0.13- 
0.2832.  The  interface-defect  levels  in  this  model  correspond 
to  the  GaAs  defect  levels.22  Their  energies  are  given  by 

£c(GaAs)  -  £*,  =  0.93  eV, 

-  £„(GaAs)  =  0.68  eV. 

Our  choice  of  band  offsets  and  defect  levels,  though  some¬ 
what  arbitrary,  is  reasonable,  but  more  importantly,  our 
conclusions  based  on  these  calculations  do  not  depend  on 
this  particular  choice. 

In  Fig.  3,  we  show  the  band  bending  in  the  ideal  Ge-GaAs 
heterostructure  (ideal  meaning  that  there  are  no  interface 
defects  or  that  the  defect  density  is  smaller  than  10‘ 1  cm"2). 
We  show  four  different  cases  corresponding  to  either  n-  or  p- 
type  GaAs  on  either  n-  or  p-type  germanium.  The  doping 
concentration  (either  donors  or  acceptors)  is  1017cm-2  in  all 


four  cases.  The  important  feature  of  the  ideal  heterojunction 
is  that  the  Fermi-level  position  can  be  moved  all  the  way 
across  the  germanium  band  gap  by  changing  the  doping.  In 
fact,  with  the  right  doping  concentration,  the  Fermi  level  can 
be  anywhere  in  the  gap  without  any  band  bending  in  either 
the  germanium  or  the  GaAs. 

In  Fig.  4,  we  show  the  result  of  introducing  interface  de¬ 
fects.  The  defect  density  in  this  figure  is  2  013  cm-2,  half  of 
them  donors  and  the  other  half  acceptors.  One  can  see  that 
the  possible  range  of  the  Fermi-level  position  at  the  interface 
decreases  substantially  as  a  result  of  this  high  defect  density. 
The  Fermi-level  position  at  the  interface  as  a  function  of 
defect  density  is  shown  in  Fig.  5.  The  four  curves  correspond 
to  the  four  combinations  of  either  n-  orp-type  germanium  on 
either  n-  or p- type  GaAs.  The  doping  in  all  four  cases  is  10 17 
cm-3.  One  can  divide  Fig.  5  into  three  regions.  For  a  defect 
density  below  1012  cm-2,  the  Fermi-level  position  at  the  in¬ 
terface  is  determined  by  the  doping  of  the  bulk  semiconduc¬ 
tors;  the  defect  density  is  too  small  to  have  any  effect.  For  a 
defect  density  of  10M  cm-2  or  more,  the  Fermi-level  position 
is  determined  by  the  interface  defects.  Between  these  two 
regions,  that  is,  for  a  defect  density  in  the  range  10‘2-1014 
cm-2,  the  Fermi-level  position  is  determined  both  by  the 
bulk  doping  and  by  the  defects. 

It  is  interesting  to  compare  the  case  of  a  thick  coverage, 
described  above,  to  either  the  thin  coverage,  described  in  the 
previous  subsection,  or  to  the  case  of  a  thick  metallic  cover¬ 
age.  27-28  The  behavior  of  the  Fermi-level  position  at  the  inter¬ 
face  in  Fig.  5  (representing  a  thick  semiconductor  coverage) 
is  similar  to  the  one  shown  in  Fig.  2  (thin  coverage).  In  both 
cases,  the  Fermi-level  position  shifts  from  being  determined 
by  the  bulk  to  being  determined  by  the  interface  defects  at  a 
defect  density  of  1012-1013  cm-  2  (in  the  case  of  a  bulk  hetero¬ 
junction,  approximately  twice  as  many  defects  are  required 
to  have  the  same  effect).  At  that  defect  density,  the  bulk 


n-Ge  n-GaAs  p-Ge  n-GaAs  n-Ge  n-GoAs  p-Ge  n  GoAs 
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Fig.  3.  The  energy  levels  near  the  interface  of  a  Ge-GaAs  heterojunction 
with  very  few  interface  defects  (less  than  10' '  cm " J).  The  four  band  struc¬ 
tures  shown  are  the  four  combinations  of  either  n-  or  p-  type  GaAs  on  either 
n-  or  p-  type  germanium.  The  doping  concentration  is  101’  cm  'on  both 
sides  of  the  interface.  The  band  offsets  correspond  to  Ref  36  The  Fermi- 
level  position  can  be  moved  throughout  the  germanium  gap  by  changing  the 
doping. 


n-Ge  p-GoAs  p-Ge  p-GaAs 


800  A  800  A 

Fig.  4.  The  energy  levels  near  the  interface  of  a  Ge-GaAs  heterojunction 
with  a  high  concentration  of  interface  defects  (10IJ  cm"1,  half  ot  them 
donors  and  half  acceptors!.  All  the  semiconductor  parameters  are  the  same 
as  in  Fig.  3.  The  interface-defect  levels  were  taken  from  Ref.  22.  The  Fermi 
level  is  not  quite  pinned  at  the  interface,  i.e.,  it  can  be  slightly  moved,  de¬ 
pending  on  the  doping,  but  the  range  of  this  movement  is  very  limited  due  to 
the  high  concentration  of  interface  defects. 
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Fig.  5.  The  Fermi-level  position  at  the  interface  of  a  Ge-GaAs  heterojunc- 
tion  as  a  function  of  interface-defect  density  (50%  acceptors  and  50%  do¬ 
nors).  The  semiconductor  parameters  (doping,  bands  offset,  and  defect  lev¬ 
els)  are  the  same  as  in  Figs.  3  and  4. 

doping  still  has  some  influence  on  the  position  of  the  Fermi 
level,  namely,  it  determines  on  which  defect  level  it  will  be 
pinned.  Only  when  the  defect  density  becomes  substantially 
higher,  does  the  influence  of  the  bulk  doping  on  the  Fermi- 
level  position  vanish.  This  behavior  should  be  contrasted 
with  that  of  a  bulk  metal-semiconductor  interface,  in  which 
the  effect  of  the  defects  is  very  small  below  1014  cm-2,  and 
the  Fermi-level  position  is  independent  of  the  semiconduc¬ 
tor  doping. 


IV.  SUMMARY  AND  CONCLUSIONS 

The  effect  of  interface  defects  on  semiconductor  hetero¬ 
junctions  was  investigated.  It  was  shown  that  a  defect  den¬ 
sity  of  approximately  1012  cm-2  (which  is  formed  at  sub- 
monolayer  coverages)  can  affect  the  Fermi-level  position  at 
the  interface  of  bulk  heterojunctions.  This  defect  density  is 
comparable  to  the  depletion  charge  in  the  bulk  semiconduc¬ 
tors.  It  was  also  shown  that  a  much  higher  density  10l3-1014 
cm-2  is  required  to  pin  the  Fermi  level  in  a  bulk  semicon¬ 
ductor  heterojunction  (compared  to  approximately  10M 
cm-2  in  the  case  of  a  metal-semiconductor  interface28). 
When  such  a  high  density  of  defects  is  present,  the  defect 
charge  is  much  higher  than  the  bulk  depletion  charge.  The 
experimental  evidence  supporting  the  unified  defect  model 
does  not  show  the  existence  of  such  a  high  defect  density. 

An  even  stronger  statement  can  be  made  regarding  the 
effect  (or  lack  thereof)  of  interface  defects  on  band  offsets. 
We  have  argued  that  if  the  defects  affect  the  band  offsets, 
then  they  have  to  satisfy  three  conditions.  Firstly,  there 
should  be  both  defect  donors  and  defect  acceptors  near  the 
interface.  Secondly,  the  layer  of  defect  donors  should  be  phy¬ 
sically  separated  from  the  layer  of  defect  acceptors  by  several 
angstroms.  Thirdly,  it  requires  a  large  density  of  such  defects 
(so  large,  in  fact,  that  it  will  affect  the  Fermi-level  position 
even  more,  thus  making  the  interface  undesirable  from  a 
device  point  of  view).  In  view  of  these  requirements,  we  can 


say  that  it  is  unlikely  that  interface  defects  determine  the 
band  offsets,  or  have  a  significant  effect  on  them. 
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Transition-metal  silicides  lattice-matched  to  silicon 
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We  have  used  a  systematic  search  to  determine  all  the  possible  transition-metal  silicides  that  are 
geometrically  lattice-matched  to  either  the  (100),  ( 1 10),  or  the  (111)  face  of  silicon.  A  short  table 
with  the  best  possible  matches  is  presented  here,  and  a  more  comprehensive  table  including 
slightly  worse  matches  is  deposited  with  the  editor. 


Recently,  there  has  been  a  considerable  interest  in  thin 
films  of  transition-metal  silicides  grown  epitaxially  on  sili¬ 
con.1'5  The  traditional  application  of  transition-metal  sili¬ 
cides  is  as  contacts  and  interconnects;  this  application  uses 
their  metallic  character  and  their  ability  to  be  oxidized  or  to 
withstand  high  temperatures.  Transition-metal  silicides  also 
have  possible  applications  in  novel  devices. 

The  problem  of  epitaxial  growth  is  complicated.  One  of 
the  factors  that  can  be  of  relevance  to  epitaxial  growth  is  the 
lattice  mismatch.  In  a  previous  paper10  we  have  described  a 
general  method  that  enables  one  to  find  systematically  all  the 
possible  matching  faces  of  any  pair  of  crystalline  materials. 
In  this  paper  we  apply  this  method  to  the  interfaces  of  transi¬ 
tion-metal  silicides  on  silicon,  and  thereby  obtain  the  most 
comprehensive  tabulation  of  such  lattice  matches  for  this 
system.  Partial  lists  were  published  before.’  "  The  silicides 
that  were  considered  are  those  of  the  three  transition  periods 
in  the  periodic  table  (columns  IIIA  through  IB)  including 
lanthanum,  but  not  the  other  lanthanides. 


The  precise  definition  of  lattice  match,  as  well  as  a  sys¬ 
tematic  procedure  to  find  them,  was  given  in  a  previous  pub¬ 
lication10;  we  shall  repeat  here  only  the  essential  definitions 
and  concepts.  By  “lattice  match”  we  mean  a  periodic  trans¬ 
lational  symmetry  of  the  interface  that  is  compatible  with 
the  crystal  structures  on  both  sides  of  that  interface.  To  un¬ 
derstand  this  definition,  refer  to  Fig.  1,  in  which  we  show  a 
hypothetical  lattice  match  betweenV3Si(  111)  and  Si(l  11).  To 
see  the  positions  of  the  atoms  relative  to  the  cubic  axes,  the 
atoms  are  shown  within  a  cube  of  20  A  on  a  side,  from  which 
one  corner  was  cut;  this  cut  comer  reveals  a  ( 1 1 1 )  face.  The 
lattice  translations  in  the  (111)  face  of  silicon  (face-centered 
cubic,  a  =  5.43 1  A)  form  a  two-dimensional  lattice.  This  lat¬ 
tice  is  shown  in  Fig.  1  as  a  fine-line  grid.  It  is  shown  both  at 
the  bottom  of  the  figure  to  demonstrate  the  translations  and 
the  atoms,  and  at  the  top  of  the  figure  for  clarity.  The  sides  of 
each  of  the  little  rhombes  (the  “unit  cells”  of  that  grid)  are 
5.431  A/>/2  =  3.840  A,  and  they  are  oriented  parallel  to  the 
[lTO],  [Oil],  or  the  [101]  directions.  Similarly,  the  lattice 


FIG  1 .  Lattice  match  of  V,Si|  111)  on  Si(  111).  The 
atomic  positions  for  both  materials  are  shown  in 
the  lower  part  of  the  figure  (VjSi  on  the  left,  sili¬ 
con  on  the  right).  The  atoms  are  shown  as  spheres 
packed  in  a  20  x  20x20- A  cube  oriented  parallel 
to  the  cubic  axes.  One  comer  of  each  cube  is  cut  m 
a  ( 1 1 1 )  plane,  and  the  viewer  is  looking  down  the 
( 1 1 1 1  direction.  The  translational  symmetry  par¬ 
allel  to  the  (111)  face  is  shown  as  a  grid  of  60" 
rhombes.  Larger  unit  cells  can  be  formed  on  those 
grids  to  obtain  a  good  match;  in  this  case  a  larger 
unit  cell  containing  three  silicon  unit  cells,  as 
shown  in  the  upper  right  comer,  matches  almost 
precisely  (0.4%)  the  V,Si  original  unit  cell. 
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translations  perpendicular  to  the  (111)  face  of  V,Si  [simple 
cubic,  0-W structure,  a  =  4.722  A)  form  a  similar  rhombic 
grid  whose  unit  cell  dimensions  are  4.722  A  ■  ,,2  =  6.678  A. 
One  can,  therefore,  form  a  superlattice  on  the  silicon  surface, 
whose  unit  cell  is  composed  to  three  silicon  unit  cells  as 
shown  in  Fig.  1 .  This  larger  unit  cell  is  oriented  such  that  its 
sides  are  parallel  to  the  [121],  [1 12],  or  [21 1]  directions,  and 
their  lengths  are  5.431  v3/2  =  6.652  A,  or  less  than  0.4%  oir 
the  VjSi  unit  cell.  These  two  unit  cells  are  shaded  in  the 
upper  part  of  Fig.  1. 

One  can  see  that  the  lattice  match,  as  defined  here,  is 
characterized  by  two  parameters:  the  mismatch  and  the 
common-unit-cell  dimensions.  In  the  example  given  in  Fig. 
1,  the  mismatch  is  0.4%  and  the  common  unit-cell  area  is 
approximately  38  A2  (38.32  for  the  silicon,  38.62  for  the 
V,Si).  In  the  most  general  case,  the  common  unit  cell  is  not 
necessarily  a  rhombus,  like  in  Fig.  1,  but  rather  a  parallelo¬ 
gram,  whose  sides  area  and  b,  with  an  acute  angle  a  between 
them.  These  three  parameters,  namely,  a,,  6,  and  «,  for  the 
substrate,  may  be  slightly  different  for  the  corresponding  a2, 
b2,  and  a2  of  the  film.  This  difference  determines  the  mis¬ 
match  (e)  of  the  match: 


The  area  of  the  common  parallelogram,  i.e., 

A  —  ab  sin  a  * 

is  the  other  parameter  relevant  to  lattice  match  (there  are 
actually  two  areas,  A ,  and  A  2,  corresponding  to  the  substrate 
and  the  film,  but  they  are  very  close  to  each  other|.  For  every 


choice  of  maximal  area  and  mismatch,  the  number  of  possi¬ 
ble  matches  is  finite,  and  all  of  them  can  be  found  in  a  w  ay 
that  was  described  in  our  previous  paper. 111  We  would  like  to 
emphasize  that  the  procedure  described  above  does  not  en¬ 
sure  that  the  silicon  atoms  in  the  substrate  are  aligned  with 
the  silicon  atoms  in  the  silicide.  It  does  ensure,  however,  that 
any  alignment  of  the  silicide  atoms  with  respect  to  the  silicon 
atoms  can  repeat  itself  periodically  along  the  interface. 

A  sample  of  our  results  is  given  in  Table  I.  In  this  table 
we  list  only  the  few  best  matches,  i.e.,  those  with  mismatch  of 
0.5%  or  less  in  all  three  parameters  a,  b,  and  a  of  the  com¬ 
mon  unit  cell,  and  unit-cell  area  of  50  A2  or  less.  Under  these 
conditions,  we  have  scanned  every  known  transition-metal 
silicide  and  tried  to  match  it  to  the  (100),  ( 1 10),  or  ( 1 1 1 )  faces 
of  silicon.  Since  the  restrictions  on  the  match  were  very  tight 
(0.5%  mismatch  and  50-A2  cell  area),  only  five  silicides  were 
found  to  match.  These  are  the  V3Si  [simple  cubic,  a  ~  4.722 
(Ref.  12)],  Ni3Si2  [trigonal,  a  =  5.617,  a  =  72°50'  (Ref.  1 3 i], 
NiSi;  [face-centered  cubic,  a  —  5.407  (Ref.  14)],  >'-Y,Su 
[hexagonal,  a  =  3.842,  c  =  4.144  (Ref.  15)],  and  RuSi  [sim¬ 
ple  cubic,  a  =  4.703  (Ref.  16)]. 

Much  more  extens:  e  tables,  similarly  structured  but 
with  less  tight  matches  (up  to  3%  mismatch  and  1 50-A2  unit¬ 
cell  area)  are  deposited  at  the  Physics  Auxiliary  Publication 
Service  (PAPS),17  and  can  be  requested  from  the  editor  of 
this  journal.  Matches  of  up  to  2%  mismatch  and  80-A2  unit¬ 
cell  area  are  shown  graphically  in  Fig.  2.  These  matches  are 
shown  on  the  mismatch  versus  unit-cell  area  plane,  and  they 
cluster  (as  expected)  near  the  upper  right  comer,  corre 
sponding  to  poor  match  (both  large  mismatch  and  large  unit¬ 
ed!  area).  For  every  pair  of  matching  faces  only  one  point 


TABLE  I  The  best  lattice  matches  of  transition-metal  silicides  on  silicon.  The  primitive  common  unit  cells  in  this  table  do  not  exceed  50  A:.  and  the 
mismatch  is  less  than  0.5%.  Under  these  conditions,  all  the  possible  matches  of  the  silicides  on  Si  1 1001. 1 1 10),  and  1 1 1 1 )  are  given  in  this  lable  For  each  possible 
match,  we  give  here  the  epitaxial  condition,  as  well  as  the  common  unit  cell  dimensions  on  each  side  of  the  interface.  The  epitaxial  condition  is  a  pair  of  crssijl 
directions,  one  on  each  side  of  the  interface,  that  w  ill  be  parallel  to  each  other  Many  such  pairs  arc  nossiblc,  and  only  one  of  them  is  given  here  The  cells’ 
dimensions  in  angstroms  and  degrees  are  given  here  for  comparison  The  mismatch  percentage  in  all  the  three  dimensions  of  the  common  unit  cell  is  giscn  in 
the  last  three  columns. 
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38. 

3.82 

10.12 

79.11 

3.84 

10.16 

79.11 

0.4 

0.4 

0.0 

NiSij  (100)/Si(  100) 

(0U|||(0ll| 

15. 

3.82 

3.82 

90.00 

3.84 

3.84 

90.00 

0.4 

0.4 

0.0 

NiSia  (221)/Si(  100) 

|1I01||[01I] 

44. 

3.82 

11.47 

90.00 

3.84 

11.52 

90.00 

0.4 

0.4 

0.0 

NiSij  ( 1 10)/Si(  1 10) 

|il°!|IUlo! 

21, 

3.82 

5.41 

90.00 

3.84 

5.43 

90.00 

0.4 

0.4 

0.0 

YjSij  (0001)/Si( III) 

[2110111(1101 

13. 

3.84 

3.84 

60.00 

3.84 

3.84 

60.00 

0.1 

0.1 

0.0 

RuSi(Ul)/Si(lll) 

[ioI]||(t2i! 

38. 

6.65 

6.65 

60.00 

6.65 

6.65 

60.00 

0.0 

0.0 

0.0 
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MISMATCH  (75) 


FIG  2.  Latnce  matches  of  transition-metal  silicides  on  silicon  1 100),  ( 1 10). 

1 1 II I.  All  (he  possible  matches  such  that  the  mismatch  is  smaller  than  2%, 
and  the  common  unit-cell  area  is  smaller  than  80  A:  are  considered.  These 
matches  are  characterized  by  mismatch  and  unit-cell  area,  and  are  shown, 
here  as  diamonds  for  matches  on  silicon  1 1001.  rectangles  for  matches  on 
silicon  1 1 101,  and  mangles  for  matches  on  silicon  (111).  Only  the  better 
matches  are  labeled.  All  the  matches  as  well  as  the  relative  orientations  can 
be  found  in  a  set  of  tables  deposited  at  the  PAPS  " 

(corresponding  to  the  smallest  cell  area)  is  shown  in  Fig.  2, 
for  lack  of  space  in  that  figure.  Different  faces  of  either  the 
silicide  or  the  silicon,  however,  are  shown  separately  on  Fig. 
2  For  example,  the  five  NiSi,  points  in  Fig.  2,  all  at  0.44%, 
correspond  to  NiSi;(  1 1 1)/Si(l  1 1),  NiSi:(100)/Si(100), 

NiSi2|110)/Si(110|,  NiSi,(51 1)/Si(l  1 1),  and'  NiSi2(221)/ 
Si(100),  respectively,  from  the  bottom  of  Fig.  2.  The  three 
faces  of  the  substrate  considered,  namely,  (100),  (110),  and 
(111)  are  shown  in  Fig.  2  as  diamonds,  rectangles,  and  trian¬ 
gles.  The  corresponding  faces  of  the  silicides  are  not  shown 
in  Fig.  2,  and  the  reader  is  referred  to  our  tables17  for  that 
information.  The  best  matches  at  the  lower  left  corner  of  Fig. 
2  are  the  same  matches  appearing  in  Table  I. 

It  is  interesting  to  compare  the  experimentally  known 
cases  of  epitaxial  growth  with  the  lattice-match  predictions 


presented  here.  Seven  transition-metal  silicides  are  currently 
known  to  grow  epitaxially  on  silicon.  These  are  NiSi-,IK"'s 
CoST,'’-23-27-9-1-’  PdjSi.1921-23'25'33-16  PtSi,21-2447  ‘ 

Pt2Si,52  51  FeSi,,55  and  CrSi2Sh  Out  of  these  seven,  the  three 
that  correspond  to  good  lattice  match,  namely,  NiSi,  (0.4% 
mismatch),  CoSi2  (1.2%),  and  Pd2Si  (1.8%-2.3%),  grow  as 
single  crystals.57  The  other  two,  PtSi  and  Pt2Si,  also  grow 
epitaxially  (in  the  sense  that  the  crystallites  in  the  film  have 
preferred  orientations  with  respect  to  the  substrate)  despite  a 
large  mismatch,  but  the  resulting  PtSi  and  Pt2Si  films  are 
polycrystalline.55  The  mismatch  between  the  Pt2Si[100]  and 
Si[110]  is  2.4%,  and  that  between  PtSi[001]  and  Si[  1 10]  is 
3.0%.  Ben  Ghozlene  et  al 55  list  several  other  preferred  ori¬ 
entations  of  the  platinum  silicides  with  respect  to  the  silicon 
substrate;  the  mismatch  in  all  of  them  is  at  least  3%.  No 
information  is  yet  available  for  the  crystalline  quality  of  epi¬ 
taxial  FeSi2  and  CrSi2,  but  they  would  correspond  to  good 
lattice  match  (0. 1  %  for  CrSi2,  1 .0%  for  FeSi2,  see  Fig.  2).  On 
the  other  hand,  epitaxial  y-Y3Si5  or  RuSi,  for  example,  have 
not  been  reported,  despite  a  very  good  lattice  match.  In  view 
of  these  experimental  results,  one  can  conclude  that  lattice- 
match  considerations  alone  are  insufficient  to  determine  all 
silicide  films  that  grow  epitaxially.  A  correlation  between 
lattice  match  and  epitaxy  does  seem  to  exist,  though.  Some 
other  systems,  however,  most  notably  silicon  (100)  on  sap¬ 
phire  (1102),  demonstrate  that  good  lattice  match  is  not  al¬ 
ways  necessary  for  a  single-crystalline  film58;  the  mismatch 
in  this  case  is  12.5%.  One  can  see,  therefore,  that  no  defini¬ 
tive  statements  regarding  heteroepitaxy  can  be  made  based 
on  the  criterion  of  lattice  match  alone.58  The  reason  is  that 
this  criterion  takes  into  account  only  the  geometry  of  the 
substrate  and  the  film,  but  not  the  chemistry  between  them. 
Hence,  only  a  combination  of  chemical  and  geometrical  con¬ 
siderations  may  yield  accurate  predictions  regarding  heter- 
oepitaxial  growth.  Despite  this  fact,  a  comprehensive  listing 
of  good  matches  is  useful  to  designate  some  silicide  films  as 
better  candidates  than  others  for  epitaxial  growth.  More¬ 
over,  the  lattice-match  criterion  is  straightforward  to  apply, 
and  one  can  scan  all  the  possible  phases  of  transition-metal 
silicides  (well  over  one  hundred  such  phases)  and  compare 
systmaticaily  ail  the  possible  crystal  faces  and  relative  orien¬ 
tations,  as  was  done  in  this  work. 
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